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Abstract 

The classic Cayley identity states that 

det(9)(detA:)* = s{s + 1) ■ ■ ■ {s + n - I) {dei Xf'^ 

where X = (xij) is an n x n matrix of indeterminates and d = (d/dxij) is the 
corresponding matrix of partial derivatives. In this paper we present straight- 
forward algebraic/combinatorial proofs of a variety of Cayley-type identities, 
both old and new. The most powerful of these proofs employ Grassmann al- 
gebra (= exterior algebra) and Grassmann-Berezin integration. Among the 
new identities proven here are a pair of "diagonal-parametrized" Cayley iden- 
tities, a pair of "Laplacian-parametrized" Cayley identities, and the "product- 
parametrized" and "border-parametrized" rectangular Cayley identities. 
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1 Introduction 



Let X = (xij) be an n X n matrix of indeterminates, and let d = (d/dxij) be 
the corresponding matrix of partial derivatives. The following beautiful identity is 
conventionalljll] attributed to Arthur Cayley (1821-1895): 

det{d) (det Xy = s{s + l)---{s + n- 1) (det X)'-^ . (1.1) 

[When n = 1 this is of course the elementary formula {d/dx)x^ = sx**~^.] A gen- 
eralization of fll.ip to arbitrary minors also holds, and is sometime^ attributed to 
Alfredo Capelh (1855-1910): if /, J C {1, . . . , n} with |/| = \ J\ = k, then 

det (9/ j) (det X)^ = s(s + 1) ■ ■ ■ (s + A; - 1) (det X)^^^ e(/, J) (det X/cjc) (1.2) 

where e{I,J) = (— l)^»ei*+I^jej-?. Analogous identities for symmetric and antisym- 
metric matrices were proved by Garding [JT] in 1948 and Shimura [3D] in 1984, re- 
spectively. 

Although these identities are essentially algebraic or combinatorial in nature, the 
simplest proofs currently available in the literature are analytic, exploiting Laplace- 
type integral representations for (detX)^ [351110] • Indeed, most of the existing alge- 



braic/combinatorial proofs [T^I5H I ^ H5 



partly because of old-fashioned notation 



I ^I1U3[[TTB] are somewhat difficult to follow. 



L 1 In this paper we would like to give straight- 
forward algebraic/combinatorial proofs of a variety of Cayley- type identities, some of 
which are known and others of which are new. The most powerful of these proofs 
employ Grassmann algebra (= exterior algebra) and Grassmann-Berezin integration]^ 
Nowadays, identities like f 1 1.1 1) are best understood as calculations of Bernstein- 
Sato type Pl32U64] for special polynomials. To see what is at issue, let P{xi, . . . , Xn) ^ 
be a polynomial in n variables with coefficients in a field K of characteristic 0. 
Then Bernstein [7] proved in 1972 that there exist a polynomial-coefficient partial 
differential operator Q{s,x,d/dx) and a polynomial b{s) ^ (both with coefficients 
in K) satisfying 

Q{s,x,d/dx)P{xy = b{s)P{xy-K (1.3) 

We call any pair {Q, b) satisfying f 1 1.3 1) a Bernstein-Sato pair for P. 

The set of all b for which there exists a Q satisfying (11.31) is easily seen to be an 
ideal in the polynomial ring K[s]. By Bernstein's theorem this ideal is nontrivial, so it 



^ But erroneously: see Section \2M\ below. 
^ Also erroneously: see again Section \2M 

^ Among the exceptions are IHJ Theorem 1.3 and Lemma 2.12] and Lemma 2.1]. 

^ We stress that our proofs are not "combinatorial" in the narrow sense of exhibiting bijections or 
exploiting double-counting. Rather, our proofs are based on straightforward algebraic manipulation 
combined with some elementary arguments of enumerative combinatorics (enumeration of permuta- 
tions by number of cycles, identities involving binomial coefficients, etc.). For this reason we have 
opted to describe our methods using the rather awkward adjective "algebraic/combinatorial". 
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is generated by a unique monic polynomial b{s), called the Bernstein-Sato polynomial 
(or b-function) of P|j Cayley-type identities thus provide Bernstein-Sato pairs for 
certain polynomials P arising from determinants. 

Bernstein-Sato pairs are especially useful in treating the problem of analytically 
continuing the distribution Pq, which can be posed as follows [5|[71-[9i[32 | H2 | [i3 | l52 [ l64] : 
Let P{xi, . . . ,Xn) ^ be a polynomial with real coefficients, and let fl C R"' be an 
open set such that P > on Q and P = on dQ. Then, for any complex number s 
satisfying Re s > 0, the function P^ is well-defined on fl and polynomially bounded, 
and thus defines a tempered distribution P^ G 5'(M'*) by the formula 



for any test function ip G iS(]R"). Furthermore, the function s ^ {P^,if) is analytic 
on the half-plane Res > 0, with complex derivative given by 



Thus Pq is a distribution-valued analytic function of s on the right half-plane. We 
want to know whether P^ can be analytically continued to the whole complex plane 
as a meromorphic function of s. This problem was first posed by I.M. Gel'fand ^2] 
at the 1954 International Congress of Mathematicians. It was answered affirmatively 
in 1969 independently by Bernstein and S.I. Gel'fand [S] and Atiyah [S], using deep 
results from algebraic geometry (Hironaka's resolution of singularities |1H])- A few 
years later, Bernstein [7] produced a much simpler proof based on using the Bernstein- 
Sato equation (11. 3p to extend P^ successively to half-planes Res > —1, Res > —2, 
etc. See e.g. P sections 7.1 and 7.3] or ^52l section 5.3] for details. 

The special case in which P is a determinant of a symmetric or hermitian matrix 
(and Q is e.g. the cone of positive-definite matrices) has been studied by several au- 
thors [T0 | [TT | [36 |[72]47it l8 m [8T] ; it plays a central role in the theory of Riesz distributions 
on Euclidean Jordan algebras (or equivalently on symmetric cones) [36', Chapter VII]. 
This case is also useful in quantum field theory in studying the analytic continuation 
of Feynman integrals to "complex space-time dimension" pilll55| [^ . In an analogous 
way, the parametrized symmetric Cayley identity (Theorem 12.111 below) will play a 
key role in studying the analytic continuation of integrals over products of spheres 
gN-i ^ "complex dimension N" [23], with the aim of giving a rigorous non- 

perturbative formulation of the correspondence found in [201 121] between spanning 

^ In the literature on Bernstein-Sato equations it is customary to shift our s by 1, i.e. write 
Q{s, X, d/dx)P{x)^~^'^ = 6(s)P(x)^, so that the usual Bernstein-Sato polynomial is our b{s + l). We 
choose here the slightly unconventional notation (|1.3I) because it seems better adapted to the Cayley 
identity pTTI) . 




(1.4) 



n 




(1.5) 



n 
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forests and the iV- vector model in statistical mechanics at = — 1. This latter appli- 
cation was, in fact, our original motivation for studying Cayley-type identities. The 
original Cayley identity f ll.ip was also rediscovered by Creutz [53] and used by him 
to compute certain invariant integrals over SU (n) that arise in lattice gauge theory. 

Many of the polynomials P treated here can also be understood as relative in- 
variants of prehomogeneous vector spaces [52||60]. When applicable, this connection 
allows the immediate identification of a suitable operator Q{d/dx) — namely, the 
dual of P itself — and provides a general proof that the corresponding b{s) satisfies 
degb = degP and is indeed (up to a constant factor) the Bernstein-Sato polynomial 
of pj^ Furthermore, this approach sometimes allows the explicit calculation of b{s) by 
means of microlocal calculus [59l[771[78l[85l[Ml99l [T00l[T06lfTU9] or other methods 

The purpose of the present paper is to give straightforward (and we hope elegant) 
algebraic/combinatorial proofs of a variety of Cayley-type identities, both old and 
new. Since our main aim is to illustrate proof techniques that may be useful in other 
contexts, we shall give, wherever possible, several alternate proofs of each result. One 
purpose of this paper is, in fact, to make propaganda among mathematicians for the 
power of Grassmann-Berezin integration as a tool for proving algebraic or combinato- 
rial identities. Among the new results in this paper are the "diagonal-parametrized" 
Cayley identities (Theorems 12.101 and 12. lip , the "Laplacian-parametrized" Cayley 
identities (Theorems 12.121 and I2.14| . and the "product-parametrized" and "border- 
parametrized" rectangular Cayley identities (Theorems 12. 161 and 12. 171) . We also give 
an elementary (though rather intricate) proof of the multi-matrix rectangular Cayley 
identity (Theorem 12. 9p that was proven recently by Sugiyama [100] . 

The plan of this paper is as follows: In Section [2] we state the identities to be 
proven and briefly discuss their interpretation. In Section [3] we give elementary alge- 
braic/combinatorial proofs of the three basic Cayley-type identities (ordinary, sym- 
metric and antisymmetric). In Section H] we give very simple proofs of these same 
identities, based on representing (detX)'* as a fermionic or bosonic Gaussian inte- 
gral]^ In Section \5\ we give alternate (and arguably even simpler) proofs, based on 
representing det(9) as a fermionic Gaussian integral; this method is very powerful 
and allows us to prove also the (considerably more difficult) "rectangular Cayley 
identities" . In Section [6] we prove the diagonal-parametrized Cayley identities, in 
Section [7] we prove the Laplacian-parametrized Cayley identities, and in Section [8] we 
prove the product-parametrized and border-parametrized rectangular Cayley identi- 

^ See [SU Corollary 6.1.1 and Theorem 6.1.1] [501 Proposition 2.22] for the first two points, 
and [52l Theorem 6.3.2] for the third. 

We are grateful to Nero Budur for explaining to us the connection between our results and the 
theory of prehomogeneous vector spaces. 

^ The proofs based on representing (det Xy as a bosonic Gaussian integral are closely related to 
the existing analytic proofs [361190) . The proofs using fermionic Gaussian integrals, by contrast, are 
really algebraic/combinatorial proofs, as fermionic "integration" is a purely algebraic/combinatorial 
construction (see Appendix IA.4I below) . 
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ties. Finally, in Section |9] we formulate some conjectures concerning the minimality 
of our Bernstein-Sato pairs. In Appendix |X] we provide a brief introduction to Grass- 
mann algebra and Grassmann-Berezin (fermionic) integration; we hope that this 
appendix will prove useful to mathematicians seeking a mathematically rigorous (but 
conceptually elementary) presentation of this powerful algebraic/combinatorial tool|^ 
In Appendix [B] we collect some formulae that will be needed in the proofs. 

We have tried hard to write this paper in a "modular" fashion, so that the reader 
can skip around according to his/her interests without having to read the whole thing. 
Indeed, after a brief perusal of Section |21 the reader can proceed directly to Section [31 
m or [5] as desired, consulting Appendices |A] and [B] as needed. 

2 Statement of main results 

Notation: We write [n] = {1, . . . ,n}. Give a matrix A, we denote its transpose 
by A^. For an invertible square matrix A, we use the shorthand A'"^ for (A"^)""" 
(which is also equal to (A^^)'^). If A = (ajj)"j=i is an n x n matrix and I, J [n], 
we denote by Ajj the submatrix of A corresponding to the rows / and the columns 
J, all kept in their original order. We write P to denote the complement of / in [n]. 
We define e(/) = (— l)l-'^l(l^l~i)/2(— i)Eig/«- jg gjgn of the permutation that takes 
the sequence 1 ■ ■ ■ n into IP when the sets / and P are each written in increasing 
order. We also define e(/, J) = e(J)e(J); it is the sign of the permutation that takes 
IF into J J''. In particular, if \I\ = \ J\, we have e(J, J) = (-l)Eie-f *+Eje,/J'. 

2.1 Ordinary, symmetric, antisymmetric and hermitian Cay- 
ley identities 

The basic Cayley-type identity is the following: 

Theorem 2.1 (ordinary Cayley identity) Let X = (xij) be annxn matrix of in- 
determinates, and let d = (d/dxij) be the corresponding matrix of partial derivatives. 
Then 

det(a) (detX)^ = s{s + I) ■■■ {s + n - 1) {det Xy-^ . (2.1) 
More generally, if I,J [n] with \I\ = \ J\ = k, then 

det{dij){detxy = s{s + 1) ■■■ {s + k - 1) {det Xy-^ e{I , J) {det Xjcjc) . (2.2) 



^ Unfortunately, most of the existing presentations of Grassmann-Berezin integration are aimed 
at physicists (see e.g. |1191 Chapter 1] for an excellent treatment) and may not meet mathematicians' 
standards of precision and rigor, or else are aimed at applications to differential geometry and hence 
involve heavier theoretical machinery than is needed for combinatorial applications. One exception 
is the brief summary given by Abdesselam 2, Section 2]; our presentation can be viewed as an 
enlargement of his. 
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Remark. Since det{d) and det{dij) are constant-coefficient differential operators, 
the matrix X can everywhere be replaced hj X + A for any fixed matrix A, and the 
identities (12.1 p / (12. 2 p remain valid. 

We consider next a version of the Cayley identity for symmetric matrices X''^™ = 
(xij). What this means is that only the variables {xij)i<i<j<n are taken as independent 
indeterminates; then Xij for i > j is regarded as a synonym for xji. 

Theorem 2.2 (symmetric Cayley identity) Let X*^^™ = (xij) be an n x n sym- 
metric matrix of indeterminates, and let (9*^^™ be the matrix whose elements are 

{d/dxii ifi=j 
\d/dx,, ifi<3 (2.3) 
\d/dxji ifi>j 

Then 

det(a^y"^) (det X'^^'y = s(s + i) ■ ■ • (^s + ^^^^ (det X^'^"')'-^ . (2.4) 
More generally, if I,J [n] with \I\ = \ J\ = k, then 

det(a57)(detX«y-)^ = s{s + l)---(^s + ^^^ (det X^H'"' e(/, J) (det X|i;;i) . 

(2.5) 

Remarks. 1. The matrix X*^^™ can everywhere be replaced by X^^^™ + A for any 
fixed symmetric matrix A. 

2. If we prefer to work over the integers rather than the rationals, it suffices to 
multiply d^^"^ by 2 and correspondingly multiply the right-hand side of (12. 4p [resp. 
fl23|) ] by 2" (resp. 2''). 

Next let us state a version of the Cayley identity for antisymmetric matrices 
j^antisym _ Hcrc ouly the variables (a;jj)i<i<j<„ are taken as independent 

indeterminates; then Xij for i > j is regarded as a synonym for —Xji, and Xa is 
regarded as a synonym for 0. As beffis antisymmetric matrices, the corresponding 
identity involves pfaffians in place of determinants: 



More precisely, alternating matrices, i.e. matrices satisfying = 0. The 

latter identity is a consequence of the former whenever the underlying ring of coefficients is an 
integral domain of characteristic 7^ 2 (so that 2x = implies x = 0), but not in general otherwise. 
See e.g. [65l sections XIII. 6 and XV. 9]. In this paper we use the term "antisymmetric" to denote 
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Theorem 2.3 (antisymmetric Cayley identity) Let x^^^^^^'^ = (xij) be a 2m x 
2m antisymmetric matrix of indeterminates, and let Q^^^^^^^y^ be the corresponding 
matrix of partial derivatives, i.e. 

ro tft = j 

^^antisym^.. _ J i < j (2.6) 

y-d/dxji ifi>j 

Then 

pf^^antisym^ ^antisym^. _ s(s + 2) • ■ ■ (s + 2m - 2) (pf X^'^^^y'")^-^ . (2.7) 

More generally, if I [2m] wzt/i |/| = 2A;, t/ien 

p^^^antisym^ ^pj^antisym^^ = S (s + 2) ■ • • (s + 2A; - 2) (pf X"'^*'^^"')"-^ e(/) (pf X^.'^f y'") . 

(2.8) 

As an immediate corollary we get a result for antisymmetric determinants: 

Corollary 2.4 (antisymmetric Cayley identity for determinants) Let X^'^*'^^™ = 
(xij) be a 2m x 2m antisymmetric matrix of indeterminates, and let g^'^^^^y^ 5e the 
corresponding matrix of partial derivatives. Then 

^g^^^antisym^ ^^^^ ^antisym^js _ (2s- l)(2s) ■ ■ ■ (2s + 2m - 2) (det X'^^"^^'^) . (2.9) 

More generally, if I [2m] with \I\ = 2k, then 

det(a^f"^"') (detX'^'^*'^^'")^ = (2s-l)(2s) ■ ■ ■ {2s+2k-2) (det X^'^^^y'")'^-^ (detX.T;^'") 

(2.10) 

Please note that in the antisymmetric case we are able at present to handle only princi- 
pal minors, i.e. we have been unable to find a general formula for det{dj^^^^^^) (det X^^^^^^"^ 
when I ^ J. 



[Zij 



Next we state a version of the Cayley identity for "hermitian" matrices Z^^™ = 
)0 By this we mean the following: We introduce indeterminates {xij)i<i<j<n and 
{yij)i<i<j<n, and define matrices X*^^™ and y^-i^tisym before; we then set 2'^'^™ = 

X^ym ^ ^yantisym ^^^^ ^herm ^ ^sym _ / 2) dy''"^'"' ■ lu tcrmS of the USUal COmplcX 

derivatives 

d 1/(9 d \ d 1/(9 d \ , , 



dzij 2 \dxij dyij J ' dzij 2 \dxij dy^j 



We put "hermitian" in quotation marks because the variables (xij) and [yij) in this identity 
are neither real nor complex numbers, but are simply indeterminates. 
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for i < j, this can be written as 

{d/ dxii if i = j 
d/dz,, iit<j (2.12) 
d/dzji ifi>j 

We then have: 

Theorem 2.5 ("hermitian" Cayley identity) Let Z^'^^"' and d^'^^^ be defined as 
above. Then 

det((9^'=™) (det Z^"™)^ = s(s + 1) • ■ ■ (s + n - 1) (det Z^<=™)^-i . (2.13) 

More generally, if I,J [n] with \I\ = \ J\ = k, then 

detid}}"^) (det Z^^™)^ = s(s + 1) ■ ■ ■ (s + A; - 1) (det Z^'"''^)'-^ e{I, J) (det Zf^'Jf) . 

(2.14) 

The resemblance of this theorem to the ordinary Cayley identity (12. ip / (12. 2p is no 
accident; indeed, the two identities are immediately interderivable. To see this, it 
suffices to notice that the action of derivatives on indeterminates is identical in the 
two cases: 

dijXiiji = Si^i'6jji (2.15a) 

^herm^herm ^ (2.15b) 

and that this relation completely determines the expressions (12. ip / (12. 2 p and fl2.13p / ( l2.14p . 
For this reason, we need not consider further the "hermitian" Cayley identityo 



2.2 Rectangular Cayley identities 

Let us now formulate some Cayley-type identities for rectangular matrices of size 
m X n with m < n. These identities are somewhat more complicated than the 
preceding ones, because the matrices appearing in the determinants (or pfaffians) are 
quadratic (or of higher order) rather than linear in the indeterminates. 

Theorem 2.6 (two-matrix rectangular Cayley identity) Let X = (xij) and 
Y = (yij) be m X n matrices of indeterminates with m < n, and let dx = {d/dxij) 
and dy = [d/dyij) be the corresponding matrices of partial derivatives. Then 

n + •^')(* + ^ - ^ + det(Xy^)^-^ . (2.16) 

Let us remark that, by contrast, in the analytic proofs [36l[90] it is convenient to consider the 
hermitian case instead of the ordinary case, as (det Zy for a positive-definite complex hermitian 
matrix Z has a simple Laplace-type integral representation. 
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More generally, if I,J [m] with \I\ = \ J\ = k, then 



det[((9x5y)7j] det(XF 



W{s + j){s + n -m + j) j det(Xr 
eiI,J) det[{XY^)j.jc 



T\s-1 



(2.17) 



If m = 1, Theorem 12.61 reduces to the easily-derived formula (Vx • Vy)(x ■ y)'* = 
s{s + n — l)(x-y)'^~"'^ for X = (a^i, . . . , x„) and y = {yi, . . . , ?/„). If m = ra, Theorem l2.6l 
can be derived by separate applications of Theorem 12. II to X and Y. In other cases 
it appears to be new. 

Theorem 2.7 (one- matrix rectangular symmetric Cayley identity) 

Let X = (xij) be an m X n matrix of indeterminates with m < n, and let d = (d/dxij) 
be the corresponding matrix of partial derivatives. Then 



' m— 1 



det(99^) det(XX' 



JJ(2s+j)(2s + n-m-l+j) j det(XX 

0=0 



T\s-1 



More generally, if I,J ^ [m] with \I\ = \ J\ = k, then 



^k-l 



det[{dd^)ij] det(XX^ 



Y[{2s + j){2s + n - m - 1 + j] 

Kj=0 



det(XX^ 



e(/, J) det[(XX^ 



(2.18) 



X 



(2.19) 



If m = 1, Theorem 12.71 reduces to the well-known formula A(x^)* = 2s (2s + n — 
2)(x^)''~^ for X = [xi, . . . ,Xn). If m = n, Theorem 12.71 can be derived from two 
applications of Theorem 12. 1[ Thegeneral case has been proven recently by several 
authors, using analytic methodso We call Theorem 12.71 a "symmetric" identity 
because the matrices XX"^ and dd"^ that appear in it are symmetric by construction. 

Finally, here is an analogue of Theorem 12.71 that involves matrices that are anti- 
symmetric rather than symmetric by construction; in place of the identity matrix we 
use the standard 2n x 2n symplectic form 



J 



1 




-1 






1 




-1 



V 



(2.20) 



Not surprisingly, this formula involves pfaffians rather than determinants: 
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See Faraut-Koranyi [36l section XVI.4], Khekalo [52l[58] and Rubin [83]. 



11 



Theorem 2.8 (one- matrix rectangular antisymmetric Cayley identity) 

Let X = (xij) be a 2mx2n matrix of indeterminates with m < n, and let d = [d/dxij] 
he the corresponding matrix of partial derivatives. Then 



' m—l 



-T\s-1 



pf((9J9^) pf(XJX^)" = |^JJ(s + 2j)(s + 2n-2m + l + 2j)j pf(XJX^ 

(2.21) 

More generally, if I [2m] with \I\ = 2k, then 



X 



(k-i 
Y[{s + 2j){s + 2n-2m + l + 2j) 
j=0 

pi{XJX^y-^e{I)pi[{XJX^)jcjc] . (2.22) 

li m = n, fl2.2ip reduces to the ordinary Cayley identity (12.1 p of size 2m since 
pi{AJA^) = det A. 

These three rectangular Cayley identities are roughly analogous to the ordinary, 
symmetric and antisymmetric Cayley identities, respectively, but their proofs are 
more intricate. 

Finally, here is a generalization of Theorems 12 . 1 1 and 12 . 61 to an arbitrary number i 
of rectangular matrices of arbitrary compatible sizes: 

Theorem 2.9 (multi-matrix rectangular Cayley identity) Fix integers i > 1 
and ni,...,ni > and write n^+i = Ui. For 1 < a < i, let X*^") be an Ua x 
Ua+i matrix of indeterminates, and let d^"'^ be the corresponding matrix of partial 
derivatives. Then 

( t m-i \ 
det(a«---9W) det(X«---XW)^ = J] J](s + n„-ni+j) det(X« ■ • ■ XW)^-^ 

\a=l j=0 / 

(2.23) 

More generally, if I,J ^ [rii] with \I\ = \ J\ = k, then 



X 



(I k-l 

det(X(i) ■ ■ ■ X^'^))^-^ e(J, J) det[(X(^) ■ ■ ■ X^'^)jcjc] . (2.24) 

We expect that there will exist "symmetric" and "antisymmetric" variants of The- 
orem 12.91 for matrix products of the form X*^^) ■ ■ ■ X'^^^X*^^^"'" ■ ■ ■ X*^^)""" or 
X(i) . . . XWX(^+1)XWT . . . X(i)T (with x(^+i) square and symmetric) for the symmet- 
ric case, and X(i) ■ ■ -X^JX^^ ■ ■ .Xd)^ or X(i) ■ . .xWx(^+i)xWt ■ ■ .Xd)^ (with 
square and antisymmetric) for the antisymmetric case (and analogous matri- 
ces of differential operators). But all good things must come to an end, and for lack 
of time we have chosen not to pursue this direction. 
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2.3 Diagonal-parametrized Cayley identities 

We would now like to formulate analogues of the ordinary and symmetric Cayley 
identities in which the diagonal elements of the matrix X are treated as parameters 
(i.e., not differentiated with respect to) and only the off-diagonal elements are treated 
as variables. 

Theorem 2.10 (diagonal-parametrized ordinary Cayley identity) LetX = (x 
be an n X n matrix of indeterminates, let a = (ai, . . . , a„) and (3 = {/3i, . . . , /3„) be 
arbitrary numbers (or symbols), and let Da^/s^s be the matrix of differential operators 
defined by 

s-ai^Xik-^-{l-ai)Y.XH-£- ifi = j 



Then 



det(D«,^,,) (det Xy = s(s + 1) ■ ■ ■ (s + n - 1) (f[ (det Xy~^ 



(2.26) 



More generally, if I,J [n] with \I\ = \ J\ = k, we have 
det((D«,^,,),j)(detX)^ = 

s(s + l)...{s + k-l) idetxy-' (n^^O f n^J^^O '(^' '^^ (^^*^^^^= 



(2.27) 



Please note that although the elements of the matrix Dct,i3,s belong to a non- 
commutative algebra of differential operators, this matrix has the special property 
that each of its elements commutes with all the elements not in its own row or column; 
therefore, the determinant is well-defined without any special ordering prescriptions. 

Remarks. 1. Expressions involving x^^ and x^jj can be understood as follows: 
we work in the Weyl algebra generated by {xij}i^j and {dij}i^j, augmented by s 
in the usual way (see Section 12. 7p as well as by the central elements yi = x^- and 
Zi = x]~^' , it being understood that X j^j^ IS 9j shorthand for yiZi. 

2. It is easy to see why f l2.26p holds with a right-hand side that is independent of 
the choice of /3. Indeed, the operators Da,/3,s and Da,/3\s are related by the similarity 
transformation 

^«,/3',« = diag(x5"''')D,,,^,,diag(xfr^0 (2-28) 

where the quantities xf^^^' commute with all entries in all matrices, so that det{Da^f3,s 
det{Dc^/3'^s)- Similar reasoning explains why the right-hand side of f l2.27p depends on 
f3 in the way it does. 
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3. It should be stressed that (12.261) provides a non-minimal Bernstein-Sato pair. 
In fact, a lower-order Bernstein-Sato pair can be obtained from f l2.27p by taking 
I = J = [n] \ {io} for any fixed io G [n]: 

det((D«,;3,,)|,„|c{i„}0(detX)^ = s{s + l) ■ ■ ■ {s + n-2) (flx.i) {detxy-' . (2.29) 

Theorem 2.11 (diagonal-parametrized symmetric Cayley identity) 

Let X^^'^ = {xij) be annxn symmetric matrix of indeterminates, let f3 = (/3i, . . . , /3„) 
be arbitrary numbers ( or symbols ), and let -D^^™ be the matrix of differential operators 
defined by 



d 



(2.30) 



ifi>j 



Then 



det(D|n(detX^n' = + + (^Hx,,^ (detX^n^"' • (2.31) 



More generally, if I,J [n] with \I\ = \ J\ = k, we have 



det((D^7),j) (detX^^-)^ 
k-1 



s(s + 



s + 



■sym ^ 



(2.32) 



Please note that in the symmetric case we are forced to take = | for all i. Note 
also that (I2.3ip provides a non-minimal Bernstein-Sato pair, and that a lower-order 
pair can be obtained from (I2.32p by taking I = J = [n] \ {zq} for any fixed io G [n]: 

det((D^7)|,,|c{,„}c)(detX^n^ = ^(^+|)---(^ + ^) (n^..)(detX^n^-'- 

(2.33) 



2.4 Laplacian-parametrized Cayley identities 

In the preceding subsection we treated the off-diagonal elements {xij}i^j as in- 
determinates and the diagonal elements Xu as parameters. Here we again treat the 
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off-diagonal elements as indeterminates, but now we use the row sums ti = ^ij 

as the parameters. (This way of writing a matrix arises in the matrix-tree theo- 
rem |2l[27ll69].) In other words, we define the row-Laplacian matrix with off-diagonal 
elements {xij}i^j and row sums 0, 

{Xij iiij^ j 

-Y.x.k if.=J (2.34) 

and the diagonal matrix T = diag(tj), and we then study det(T + x™'^"^^p) where 
the ti are treated as parameters. We shall need the matrix of differential operators 

( / dx ■ ■ if z 7^ 7 

^^row-Lap).. ^ } ' ^ ■> (2.35) 

I li i = j 

Theorem 2.12 (Laplacian-parametrized ordinary Cayley identity) 

Let x™™"'"'^^, T and be n x n matrices defined as above. Let U be the n x n 

matrix with all entries equal to 1. Then 

[det([/ + 5™"-^^P) - det(9™"-^'^P)] (det(T + x^^'^^^p))^ 

= (^ti) s(s + l)---(s + n-2)(det(T + X""-L^P))"-^ (2.36) 

i 

Let us recall that, by the matrix-tree theorem [21E71En], det(T + X^^-^^p) is the 
generating polynomial of rooted directed spanning forests on the vertex set [n], with 
weight riieiJ^* 'Wi'j^E{F)i~^ii) ^ (rooted directed) forest with roots at the vertices 
i & R and edges ij G E{F), directed towards the roots. In particular, the term linear 
in ti — whose coefficient is the principal minor of order n — 1, det(X™""'"'^P){j}c{j}.c — 
enumerates the directed spanning trees rooted at i. Taking this limit in fl2.36p . we 
obtain: 

Corollary 2.13 (Cayley identity for the directed-spanning-tree polynomial) 

For each i E [n], we have 

[det{U + 9™"-L"P) - det(9™"-^"P)] (det(X""-L"P){i}c{,|c)* 

= s(s + 1) • ■ ■ (s + n - 2) (det(X''°™-^^P)|i}c{i|c)^-i . (2.37) 

We also have an analogous identity for symmetric Laplacian-parametrized matri- 
ces. What this means is that we introduce indeterminates {xjj}i<j<j<„ and regard 
Xij for z > j as a synonym for xji] we then define the symmetric Laplacian matrix 



(j^sym-Lap-j 



Xij Hi < j 

Xji if « > j (2.38) 

k^i 
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and the corresponding matrix of partial derivatives 

{d/dxij if i < j 
d/dx,, iii>j (2.39) 
if ^ = J 

We then have: 

Theorem 2.14 (Laplacian-parametrized symmetric Cayley identity) 

Let X'^y™"^^^, T and (}sym-Lap ^ x n matrices defined as above. Let U he the n x n 
matrix with all entries equal to 1. Then 

[det(?7 + 9'y°^-L^P) - det(9'y'^-^^P)] (det(T + x^^-^-^^p))^ 

= (^ti)2s(2s + l)---(2s + n-2)(det(r + X^y'"-L"P))^-^ (2.40) 

i 

Similarly, det(T + X^^'^'^^'^) is the generating polynomial of rooted (undirected) 
spanning forests on the vertex set [n], with weight HiGi?^* nije£{F)(~-^ii) ^^^^ ^ (rooted 
undirected) forest with roots at the vertices i & R and edges ij G E{F). In particular, 
the term linear in ti — whose coefficient is det(X'^^™"^^P){j}c{j}c — is independent of 
i and enumerates the spanning trees. We thus obtain: 

Corollary 2.15 (Cayley identity for the spanning-tree polynomial) For each 
i E [n], we have 

[det{U + a^y'^-L'^P) - det(a^y'"-L'^P)] (det(X^y"^-L'^P){i|c{i|e)^ 

= 2s(2s + 1) ■ ■ ■ (2s + n - 2) (det(X^y"^-L^P){i|e|i}c)^-i . (2.41) 

See pi] for many interesting additional properties of the spanning-tree polynomial. 

2.5 Product-parametrized and border-parametrized rectan- 
gular Cayley identities 

Here we will present two curious Cayley identities for rectangular matrices that are 
much simpler than the identities presented in Section [2^2} because the indeterminates 
occur linearly rather than quadratically in the argument of the determinant. 

So let X = (xij) be an m X n matrix of indeterminates with m < n, and let 
d = (d/dxij) be the corresponding matrix of partial derivatives. One easy way to 
obtain square (m x m) matrices from X and d is to right-multiply them by n x m 
matrices A and B, respectively. We then have the following Cayley-type identity, in 
which A and B occur only as parameters: 
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Theorem 2.16 (product-parametrized rectangular Cayley identity) 

Let X = (xij) be an m X n matrix of indeterminates with m < n, and let d = {d/dxij) 
he the corresponding matrix of partial derivatives. Let A = (aij) and B = {pij) he 
n X m matrices of constants. Then 

det{dB) (det XA)' = det{A^B) s(s + 1) ■ ■ ■ (s + m - 1) (det XA)'-^ . (2.42) 

More generally, if I,J [m] with \I\ = \ J\ = k, then 

det{{dB)jj) (det XA)' = det(MA) s{s + 1) ■ ■ ■ {s + k - 1) (det XA)'-^ . (2.43) 

where M is anmxn matrix defined as follows: if I = {zi, . . . , i^} and J = {ji, . . . , jk} 
in increasing order, then 

\ Xap if a ^ I 
M^p = I / . (2.44) 

Note that if = m, then f l2.42p reduces to the ordinary Cayley identity (12. ip 
multiplied on both sides by (det 5) (det A) while fl2.43p specialized to A = B = 1^ 
gives the all- minors identity (12. 2p . In the general case m < n, fl2.43p reduces to (I2.42p 
when I = J = [m] (since we then have M = B'^), while (I2.43P reduces to the trivial 
identity (detXA)'^ = (detXA)'^ when I = .J = (since we then have M = X). 

A second easy way to complete X and d to square {n x n) matrices is to adjoin 
n — m rows of constants at the bottom: 

where A and B are [n — m) x n matrices of constants. We can think of X and d 
as "bordered" matrices obtained by filling out X and d. We then have the following 
Cayley-type identity, in which A and B again occur only as parameters: 

Theorem 2.17 (border-parametrized rectangular Cayley identity) Let X = 

(xij) he an m X n matrix of indeterminates with m < n, let d = (d/dxij) he the 
corresponding matrix of partial derivatives, and let A and B he [n — m) x n matrices 
of constants. Define X and d as in 1^2.45^ - Then 

detid) (det Xy = det{AB^) s(s + 1) ■ • ■ (s + m - 1) (det Xy~^ . (2.46) 

When m = n this formula reduces to the ordinary Cayley identity (12.11) . if we make 
the convention that the determinant of an empty matrix is 1. 

Please note that, by Laplace expansion, detX is a linear combination of m x m 
minors of X (and likewise for det d). Indeed, when n — m = l one can obtain all such 
linear combinations in this way, by suitable choice of the row vector A; for n — m>2 
one obtains in general a subset of such linear combinations. It is striking that the 
form of the identity — and in particular the polynomial h{s) occurring in it — does 
not depend on the choice of the matrices A and B. 

In Section [8] we will prove the "product-parametrized" and "border-parametrized" 
identities and then explain the close relationship between them. 
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2.6 Historical remarks 



As noted in the Introduction, the identity (11.11) is conventionally attributed to 
Arthur Cayley (1821-1895); the generalization (ll.2p to arbitrary minors is sometimes 
attributed to Alfredo CapeUi (1855-1910). The trouble is, neither ([HI]) nor ([L2D 
occurs anywhere — as far as we can tell — in the Collected Papers of Cayley [26j . 
Nor are we able to find these formulae in any of the relevant works of CapeUi [T514T9] . 
The operator Q = det{d) was indeed introduced by Cayley on the second page of his 
famous 1846 paper on invariants [2S]; it became known as Cayley's f2-process and went 
on to play an important role in classical invariant theory (see e.g. [3 H [39 | [75 | [89 ]1 113] ). 
But we strongly doubt that Cayley ever knew (II. ip . 

A detailed history of (II. ip and (II. 2p will be presented elsewhere |^. Suffice it 
to say that the special case for 2x2 matrices appears already in the 1872 book of 
Alfred Clebsch (1833-1872) on the invariant theory of binary forms p. 20]. But 
even for n = 3, the first unambiguous statement of which we are aware appears in 
an 1890 paper of Giulio Vivanti (1859-1949) [T08]FI And amazingly, in this very 
first paper, Vivanti proves not only the basic "Cayley" identity (II. ip but also the 
generalization (II. 2p for minors, for completely general n and k albeit only in the case 
I = J (i.e., principal minors). In fact, his inductive method of proof works only 
because he is handling the "all-principal-minors" version; it would not work for the 
"simple" identity (II. ip alone. 

Proofs (by direct computation) of the "Cayley" identity (II. ip for n = 3 can be 
found in the early-twentieth-century books of Grace and Young (1903) |16], Weitzen- 
bock (1923) [ng and Turnbull (1928) pTOH]. Weitzenbock also states dH]) without 
proof for general n, saying that it is obtained "by completely analogous calcula- 
tion" |112| p. 16]; similarly, Turnbull states both (II. ip and (II. 2p for general n and 
leaves them as exercises for the reader |103| pp. 114-116]. (Unfortunately, TurnbuU's 
old-fashioned notation is very difficult to follow.) We are not convinced that the 
extension from = 3 to general n is quite so trivial as these authors imply. We 
will, in any case, provide an elementary algebraic/combinatorial proof of (ll.ip / f lTT^ 
in Section [3?T] below. 

The symmetric analogues (12. 4p / ( 123|) are due to Carding in 1948 [H]; see also [361 
Proposition VII. 1.4]. The antisymmetric analogue (12. 7p is due to Shimura [90] in 1984; 
see again Proposition VII. 1.4], where the quaternionic hermitian determinant 
(Moore determinant [4j) is equivalent to a pfaffian, and see also [HH Corollary 3.13]. 

Shimura [9^ also gives generalizations of all these formulae in which det{d) is 
replaced by other homogeneous differential operators. Similarly, Rubenthaler and 
Schiffmann [82] especially Section 5] and Faraut and Koranyi [36] Proposition VII. 1.6] 
give generalizations in which both det(9) and detX are replaced by suitable prod- 



This paper is also cited in Muir's massive annotated bibliography of work on the theory of 
determinants [70l vol. 4, p. 479]. Indeed, it was thanks to Muir that we discovered Vivanti's paper. 
Malek Abdesselam has also drawn our attention to the papers of Clebsch (1861) [28i PP- 7-14] and 
Gordan (1872) [45^ pp. 107-116], where formulae closely related to (|l.ip can be found. 
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ucts of leading principal minors. All these proofs are analytic in nature. Elegant 
algebraic/combinatorial proofs of identities in which detX is replaced by a product 
of minors have been given by Canfield, Williamson and Evans [HI Theorem 1.3 and 
Lemma 2.12] [IISl Theorem 4.1]0 

The one-matrix rectangular symmetric Cayley identity fl2.18p has been proven 
recently by Faraut and Koranyi [361 section XVI. 4], Khekalo [57l|58j and Rubin [83] , 
using analytic methods. 

Many of the polynomials P treated here can also be understood as relative invari- 
ants of prehomogeneous vector spaces [5^1^: and the corresponding 6-functions have 
been computed in that context, mostly by means of microlocal calculus [59 | l86 |[T00] . 
See [601 Appendix] for a table of all irreducible reduced prehomogeneous vector spaces 
and their relative invariants and 6-functions. All but one of the "generic" cases in 
this table (i.e., those involving matrices with one or more arbitrary dimensions) corre- 
spond to identities treated here: namely, the ordinary, symmetric and antisymmetric 
Cayley identities [f l2.ip . (12. 4p and (12. 7p ] correspond to cases (1), (2) and (3), respec- 
tively, while the one-matrix rectangular symmetric and antisymmetric Cayley iden- 
tities [ (I2.18P and (I2.2ip ] correspond to cases (15) and (13). The 6-functions hsted in 
that table of course agree with ours0 (The final "generic" case (30) involves matrices 
X of order 3 x 2?7, with P{X) = tr(X JX""")^, which falls outside our methods since it is 
neither a determinant nor a pfaffian.) Finally, the most difficult result obtained in the 
present paper — namely, the multi-matrix rectangular Cayley identity identity (I2.23P 
— has very recently been proven independently by Sugiyama }100 -. Theorem 0.1] in the 
context of prehomogeneous vector spaces associated to equioriented quivers of type A; 
his proof uses a decomposition formula found earlier by himself and F. Sato [85] . In- 
deed, Sugiyama proved an even more general result jlOOl Theorem 3.4], applying to 
quivers of type A with arbitrary orientation!^ 

It is worth stressing that the Cayley identity (II. ip — though not, as far as we 
can tell, the all-minors version (II. 2p — is an immediate consequence of a deeper 
identity due to Capelli [TBHTS] . in which the operator H = (detX)(det9) is repre- 
sented as a noncommutative determinant involving the Ql{n) generators X^d: see 
e.g. |1071 p. 53], [nH pp. 569-570] or [221 Appendix] for the easy deduction of Cayley 
from Capelli. Likewise, the symmetric Cayley identity (12. 4 p follows from a symmetric 
Capelli- type identity due to TurnbuU |104] (see also [IIOJ ). and the antisymmetric 
Cayley identity (12. 7p follows from an antisymmetric Capelli-type identity due inde- 
pendently to Kostant and Sahi [53] and to Howe and Umeda [HI] (see also [EI])- 

The one-matrix rectangular symmetric Cayley identity (I2.18P for m = 1 follows 
from a Capelli-type identity given in |1131 pp. 291-293] and [107^ p. 61]. For m = n 



See also TurnbuU [105 for a similar result, but expressed in diflieult-to-follow notation. 

After making the translation of conventions s ^ s + 1: see footnote |S] above. See also F. Sato 
and Sugiyama ^85ji Section 3.1 and Lemma 4.2] for an alternate approach to (12. ip . (|2.7p and (|2.18p . 

We are grateful to Nero Budur for explaining to us the connection between our results and the 
theory of prehomogeneous vector spaces, and for drawing our attention to the work of Sugiyama |100j. 
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it of course follows from the ordinary Capelli identity. For 2<m<n — Iwedo not 
know any Capelli-type identity. 

Proofs of the Capelli-type identities based on group-representation theory have 
been given by Howe and Umeda (HDIEI]- Combinatorial proofs of the Capelli and 
Turnbull identities have been given by Foata and Zeilberger [37j. We have recently 
given very simple algebraic proofs of these same identities as well as some generaliza- 
tions [22]. See also Weyl [HSl pp. 39-42] and Fulton-Harris |39l Appendix F.3] for 
more traditional proofs. Further information on Capelli-type identities can be found 

in mmnz]. 

2.7 Some algebraic preliminaries 

A few words are needed about how the identities ( 12. ip ff. — or more generally, 
Bernstein-type identities of the form 

Q{s,x,d/dx)P{xY = b{s)P{xy-^ (2.47) 

where to be interpreted. On the one hand, they can be 

interpreted as analytic identities for functions of real or complex variables xi, . . . ,Xn, 
where s is a real or complex number; here P(a;)* denotes any fixed branch on any 
open subset of or where it is well-defined. Alternatively, these formulae can 
be regarded as purely algebraic identities, in several different ways: 

1) For integer s > 1, as an identity in the ring R[xi, . . . ,Xn] of polynomials in 
the indeterminates xi, . . . , x„ with coefficients in some commutative ring R (for 
instance, R could be Z, Q, M or C)@ 

2) For any integer s (positive or negative), as an identity in the field K{xi, . . . , Xn) 
of rational fractions in the indeterminates Xu.. . ,Xn with coefficients in some 
field K (for instance, K could be Q, M or C)l3 

3) For s interpreted symbolically, as an identity in a module defined as follows [221 
pp. 93-94] [Ml PP- 96 ff.]: Let be a field of characteristic 0, let Xi , . . . , Xn 
and s be indeterminates, and let y4„(i^')[s] be the i^'-algebra generated by 
Xi, . . . , Xn, d/ dxi, . . . ,d/ dxn and s with the usual commutation relations. (That 
is, it is the algebra of differential operators with respect to Xi, . . . , in which 
the coefficients are polynomials in xi, . . . , x„, s with coefficients in K.) Now fix 
a nonzero polynomial P G K[xi, . . . ,x„], and let ii"[x, s,P~^] denote the ring 
of rational fractions in the indeterminates Xi, . . . , x^, s whose denominators are 

When b{s) contains fractions [e.g. ((2^ / ((2?5l) and ((23T1) / ((232)) ] we should assume that the 
coefficient ring R contains those fractions [i.e. in this case i]. 

1^ When b{s) contains fractions [e.g. (|2^ / (|2?5l) and (|2?3T1) / (|232l) ] we should assume that the 
coefficient field K contains those fractions [i.e. in this case Usually we will take if to be a field 
of characteristic 0, so that K contains the rationals Q as a subfield. 
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powers of P. (It is a subring of the field K{xi, . . . , a:„, s) of all rational fractions 
in xi, . . . , Xn, s.) Let K[x, s, P~^]P'^ be the free K[x, s, P~^]-module consisting 
of objects of the form /P** where / G K[x, s, P~^]; here is treated as a formal 
symbol. We can define formal differentiation by 



where df /dxi is the standard formal derivative of a rational fraction. This 
differentiation is easily extended to an action of A„(K)[s] on s, P^^jP'*, 
making the latter into a left Ar, (K) [g]-moduleP°l 

Let us now show that all these interpretations are equivalent. 

We begin by recalling some elementary facts. Let p{xi, . . . polynomial 
with coefficients in some commutative ring R, and let di be the degree of p with respect 
to the variable Xi. Suppose that there exist sets Xi, . . . ,X„ C R with > di for 
all i, such that p{xi, . . . ,Xn) = whenever {xi, . . . ,Xn) G Xi x . . . x Xn- Then 
p must be the zero polynomial, i.e. all its coefficients are zero. Note in particular 
that if the sets Xi, . . . , Xn are infinite, then this reasoning applies to polynomials of 
arbitrary degree. As a special case of this, if P = M or C and p{xi, . . . , Xn) = for 
all ) lying in some nonempty open set U of M" or C", then p must be 

the zero polynomial. 

Now let K he a. field of characteristic 0, and let P(x), Q{s,x,d/dx) and b{s) be 
polynomials with coefficients in K [as always we use the shorthand x = (xi, . . . , x„)]. 
Then elementary algebraic manipulations allow us to write 

Q{s,x,d/dx)P{xy - b{s)P{xy-^ = R{s,x)P{xy^'^ (2.49) 

for some polynomial R{s, x) and some integer m > 0. Combining this fact with the 
preceding observations, we obtain immediately the following two propositions: 



More generally, we can proceed as follows: Let R be an integral domain and let A be an abelian 
group. We then define to be the commutative ring with identity generated by the symbols x° 
(x e i?, a e A) subject to the relations x'^x^ = x°-^^, x^y" = (xy)"' and a;° — 1. In particular, if A 
contains the integers as a subgroup, then we can consider i? as a subring of R^ by identifying x R 
with x^ e R^. 

Now suppose that i? is a polynomial ring S[xi, . . . , a;„] where S is an integral domain of charac- 
teristic 0, and that A is a subgroup of the additive group of S (where 1 e Z C A is identified with 
1 & S). Then we can define an action of the differential operators d/dxi on R^ by 

— (p^) = a — P''-^ 
dxi dxi 

together with the usual product rule. This makes i?'^ into a left A„(S')-module [where An{S) is the 
Weyl algebra in n variables over S] . 

In order to handle Bernstein-type identities, we will introduce an indeterminate s and take A = 
7L + si and S = K\s\ for some field K. Then we will work within the submodule of R^ consisting of 
elements of the form jP^^"- for / e K\x\., . . . , a;„, s], a G Z and some fixed P g K\x\^ . . . , x„]. This 
submodule is isomorphic to K\x^s^P~^\P^ . 
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Proposition 2.18 (Equivalence theorem for symbolic s) Let K be a field of char- 
acteristic 0, and let P, Q, h and R he as before. Then the following are equivalent: 

(a) holds as an algebraic identity in the module K[x , s , P~^]P'^ , symbolically 
in s. 

(b) R{s,x) = in the polynomial ring K[s,x]. 
Furthermore, if K is infinite, then (a)-(b) are equivalent to: 

(c) For infinitely many s & K, there exist infinite sets Xi, . . . ,X„ C K (possibly 
depending on s), such that R{s, Xi, . . . , Xn) = for (xi, . . . , x„) G Xi x . . . x X„ 
and the given value of s. 

Finally, if K = 'R or C, then (a)-(c) are equivalent to: 

(d) For infinitely many s & K, there exist a nonempty open set U C K"^ (possibly 
depending on s) and a branch of P{x)^ defined on U such that ^2.4'T\) holds for 
all X & U. 

(e) For every s E K, every nonempty open set U C K"- and every branch of P{x)^ 
defined on U , ( [^.^Tp holds for all x E U. 

Proposition 2.19 (Equivalence theorem for fixed s) Let K be a field of char- 
acteristic 0, let P, Q, b and R be as before, and fix some element s E K . Then the 
following are equivalent: 

(b) R{s,x) = in the polynomial ring K[x] (for the given value of s). 
Furthermore, if K is infinite, then (a)-(b) are equivalent to: 

(c) There exist infinite sets Xi,...,X„ C K such that R{s,xi, . . . ,Xn) = for 
{xi, . . . , Xn) G Xi X . . . X Xn (for the given value of s). 

Finally, if K = M. or C, then (a)-(c) are equivalent to: 

(d) There exist a nonempty open set U C K"- and a branch of P{x)^ defined on U 
such that ^2-4'^ holds for all x E U (for the given value of s). 

(e) For every nonempty open set U C and every branch of P{xY defined on U , 

holds for all X eU (for the given value of s). 

In particular, Proposition 12.181 shows that it suffices to prove f l2.47p for infinitely 
many positive or negative integers s, using the elementary interpretation (1) or (2) 
above; it then holds automatically for arbitrary s as an identity in the module 
K[x, s, P~^]P^ and as an analytic identity. Likewise, it suffices to specialize Xi, . . . , a;„ 
to real or complex variables and to prove f l2.47p for some nonempty open set in 
or C"". In what follows, we shall repeatedly take advantage of these simplifications. 
(Many previous authors — especially the earlier ones — have done so as well, but 
without making Proposition 12.181 explicit.) 
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3 Elementary proofs of Cayley-type identities 



In this section we give proofs of the three main Cayley-type identities for square 
matrices (Theorems 12. 11 - 12. 3p that use nothing but elementary properties of determi- 
nants (notably, Jacobi's identity for cofactors) along with the elementary formulae 
for the derivative of a product or a power. 

The general situation we have to handle in all three cases is as follows: Let F be 
a finite index set, let (i?^)-ygr be given n x n matrices with elements in some field K, 
and let (x^)-ygr be indeterminates. Now define the matrix X = Yl-yer^i-^i- If ^ = 
{Ay)^^r is a K-valued vector, we write Oa = J^yer ^-y^/^^i = J^-yer ^^^-y^ 

so that OaX = Ea- We need a formula for successive derivatives of (detX)**: 

Lemma 3.1 Let (i?^)^gr benxn matrices, let (x^)-^gr be indeterminates, and define 
X = J2-yeT 3^7-^7 ■ ^(^"1^ (^"iT-y sequence Ai, . . . , of K -valued vectors, we have 

k \ 

IIOaA (detxy 

i=l J 

= (-l)'=(detXr5^(-s)#(^^^i-°^-) n tr {X''Ea^^...X-'Ea J i3.1e.) 

'''^Sk cycles of r 

C = (ai, . . .ai) 

= (detX)^^sgn(r)s#(^y^i^^°^-) J] ti {X'^Ea^^ ... X~'Ea J . (3.1b) 

TG^fc cycles of t 

C = {ai,...ae) 

Proof. By induction on k. The case k = 1 follows from Cramer's rule or alternatively 
from the relation 

(detX)' = exp(s trlogX) . (3.2) 
For the inductive step we shall need the identity 

dA{X~') = -X-\dAX)X-' , (3.3) 

which follows from dA{X~^X) = 0. Assume now that the theorem is valid for k and 
let us apply dAk+i- When this derivative hits (detX)*, it creates a new cycle {k + 1), 
with prefactor s = (— 1)(— s); these terms correspond to permutations r G Sk+i in 
which the element A; + 1 is fixed. Alternatively, the derivative can hit one of the X~^ 
factors in one of the traces; by fl3.3p this inserts the element k + 1 into one of the 
existing cycles at an arbitrary position, and produces an extra factor —1; these terms 
correspond to permutations r G Sk+i in which the element -|- 1 is not fixed. □ 

For 1 < i, j < n, let E^^ be the matrix with a 1 in position ij and zeros elsewhere, 

i.e. 

{E'^),,j, = S,,i,S,j, . (3.4) 
We will express the matrices E^ in each of our three cases in terms of the E^K 
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3.1 Ordinary Cayley identity 



Proof of Theorem [2TTJ In this case the index set F is simply [n] x [n], and we write 
^ = TH,j=i XijE'K Now let I = {ii, ...,ik} with ii < . . . < and J = {ji, . . .,jk} 
with ji < ■ ■ ■ < jk, so that 

detidjj) = J2 sgn(a) n ■ (3-5) 

For each fixed a G Sk, we apply Lemma [3?T] with Oa^ = In the traces we 

have Ea = E^°"'^'^^°'r) and hence 

tT(X"'E'"^^^f-iK..X-'E'-^^^^^-e)) = Xr] Xr] ...X;] ■ (3.6a) 

I 

nX-^ , , (3.6b) 

r=l 

where = (X^^)'^ and we have used the fact that, for r as in Lemma [3. ![ r(aj) = 
ttj+i for i = 1, 1 and r(a^) = ai. We can now combine all the different traces 

into a single product. We obtain 

k 

dei{d,j) (det Xy = (det Xf ^ sgn(a) sgn(r) s*^^^^'^^ °^ X{{X-^\,r,3.ir, ■ 

(T€Sk reSk r=l 

(3.7) 

Let us now define the permutation n = a o t ^ and change variables from (a, r) to 
(r, vr), using sgn((T) sgn(r) = sgn(7r). The product over r can be written equivalently 
as a product over t = r(r). We have 

k 

Y,^g^{'rr)l[{X'^)uj.i,, = det{{X-^)u) = [detxy' e{I, J) {detXi.j.) (3.8) 

neSk t=i 

by Jacobi's identity [Lemma lA.l( e)]. while 

^ ^#(cyclcsofr) ^ + 1) ■ ■ ■ {s + k - 1) (3.9) 

(see e.g. jlTi p. 263, eq. (6.11)] or [521 Proposition 1.3.4] for this well-known equality). 
□ 
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3.2 Symmetric Cayley identity 



Proof of Theorem 12.21 In this case the index set F consists of ordered pairs 
E [n] X [n] with i < j, and we write X^^^^ = J2i<jXijiE'^ + E^') + Y^i^aE''- 

Then 

^^ym^sym _ ^{^E'^ + E^') (3.10) 

in all three cases {i < j, i > j and i = j). Now let us apply det(9jj™) to (detX^^™)^ 
and compare to what we had in the previous proof. On the one hand we have a factor 
2~'' coming from the k derivatives. On the other hand, each E^i is now a sum of two 
terms E^^ and E^^; in each cycle (ai, . . . , ai) of r, the argument of the trace becomes 

j^-i _g'i<T(c«i)««i j (_£;*"2i<T(c2) + _E'i<T(c«2)««2) . . . (3-11) 

(to lighten the notation we have written X instead of X*^^™). We therefore need to 
sum over all the 2' ways of choosing E^^ or E^'^ in each factor within the given trace 
(hence 2'' choices overall). Performing the trace, we will obtain terms of the form 

Let us now fix one of the 2^ choices and sum over the permutation a. If one or more 
of the factors is of the form ■ , then the sum over a will vanish because the 

exchange between a{ar) and a{ar+i) takes a —1 from sgn((T). Therefore, in each cycle 
of r there are only two nonvanishing contributions, corresponding to the two ways 
of coherently orienting the cycle. One of these has X~"^ (as in the previous proof) 
and the other has X~^, but these are in fact equal since X is symmetric. We thus 
have, compared to the previous proof, an extra factor 2'^^^^^^'^'^ °^ We now change 
variables, as before, in the sum over permutations. The sum over vr gives fl3.8p exactly 
as before, while the sum over r now gives 



Y^^2s)*^cyclesolr) = + 1) . . . + (3.13) 



by O. □ 



3.3 Antisymmetric Cayley identity 

Let us recall the definition of the pfaffian of a 2n x 2n antisymmetric matrix: 

Pf^ = 5Z Sgn(cr) a^(l)<x{2) ■ ■ ■ aa(2n-l)a(2n) • (3.14) 
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Proof of Theorem 12.31 In this case the index set F consists of ordered pairs 
e [n] X [n] with i < j, and we write X^"*''^^'" = Xij{E'^ - E^'). Then 

^antisymji^antisym _ j^ij _ j^ji j^g^ 

Now let / = {«!,..., «2fc} with ii < ... < i2k, and let us apply pf{d^f'^'^^^) to 
^pf;ss:antisym)s ^ (^^^ ^antisym)s/2. ^gj^g ^j^g representation (Km for piidjf"^""), we 

obtain 



pf(^antisym) ^antisym)s = ^ ^ Sgn(a) (fj ^— ^- ) {det X^^^'^n'^' . 



d_ 



(3.16) 

Now apply Lemma [3TT] as before. In each cycle (ai, . . . , ai) of r, the argument of the 
trace becomes 

(_£;V(2«i-l)i<T(2c«i) _ _g'V(2«i)*<7(2c«i-l)) Ji^-l (_g'V(2«2-l)*<T(2Q2) _ _f;«<T(2c2) V(2c«2 -1) ) .... 

(3.17) 

Once again we have 2^ choices in the E factors; but here these choices correspond 
simply to pre-multiplying a by one of the 2^ permutations that leave fixed the pairs 
{1, 2}, . . . , {2k — 1, 2k}] therefore, after summing over a we simply get a factor 2^. 
Let us now introduce the permutation a'^ defined by 

a^{2r - 1) = (T(2r(r) - 1) (3.18a) 
(T^(2r) = a{2r) (3.18b) 

For each r G Sk, the map cr i— t- o""^ is an automorphism of and satisfies sgn(cr'^) = 
sgn((T) sgn(r). We have 

k 

o-GS2fc T&Sk '■=1 

(3.19) 

We now define tt = cr"^ and change variables from (cr, r) to (r, tt) as before. The sum 
over TT gives 

^ E sg^(^)ri(^~^)M..-.,M..) = Pf((^'^)//) = (pfX)-^e(/)(pfX,.,.) 

(3.20) 

by the pfaffian version of Jacobi's identity [cf. ( 1A.21I) ]. while the sum over r gives 

2'= ^(s/2)#(=y^'<''^°f") = s(s + 2)---(s + 2fc-2) (3.21) 

by dSH). □ 
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4 Proofs of Cay ley- type identities by Grassmann/ 
Gaussian representation of (detX 



Let us now give simple proofs of Theorems 12 . H - [2]3| based on representing (det X)^ 
as a fermionic or bosonic Gaussian integral. A brief introduction to fermionic and 
bosonic Gaussian integration can be found in Appendix lAl 

4.1 Ordinary Cayley identity 

Proof of Theorem 12.11 Assume that s is a positive integer, and let us introduce 
Grassmann variables ■^/'j"'' , ■^j-"'* for i = 1, . . . , n and a = 1, . . . , s. We can then write 

(detX)^ = jv{ip,i,) e"^^^ , (4.1) 

where we have used the shorthand 

= E E ^'t^^'-.^f^ ■ (4-2) 

Now let / = {ii, . . . , with ii < . . . < ik and J = {ji, . . .,jk} with ji < . . . < jk, 
so that 

det(a,,) = J2 ^g^(^) n ■ (4.3) 
Applying this to f l4.ip . we obtain 

det{djj){detxr = ^V{^,^|J) E (ll^l'Vit)) ^'^'"'^ • (4-4) 

o-eSfc Q!i,...,Ofe=l ^r=l ^ 

When X is an invertible real or complex matrix, Wick's theorem for "complex" 
fermions (Theorem IA.16P gives 



/■p(^,^) (nV^frVt.O = (detX)^ J]sgn(r)n(X-T),.,^^^^^„5„,„^^^, . 

(4.5) 

Let us now define the permutation n = a o t and change variables from (a, r) to 
(r, tt), using sgn(cr) sgn(r) = sgn(7r). We then have 

k 

5^sgn(7r)J](X-T),^^.^^^^ ^ det((X-T),^) = (det X)"^ e(/, J) (det X,c^.) (4.6) 
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by Jacobi's identity, while 

s k 

E E = E^*^^'^''^°'^^ = s{s + l)---{s + k-l) (4.7) 

-rgcSfc «iv,"fc = l '" = 1 Tg5fc 

by (13. 9p . This proves (12.21) for integer s > 1, when X is an invertible real or complex 
matrix. By Proposition I2.18i this is sufficient to prove the identity. □ 

Alternate proof of Theorem 12.11 Instead of using "complex" fermions, we 
can use complex bosons. So assume that s = —m where m is a positive integer, 
and that X is a complex matrix whose hermitian part is positive-definite; and let 
us introduce bosonic variables <^j-"'' for i = 1, . . . , n and a = 1, . . . , m. We then 
have 

(detX)-'" = jv{ip,ip) e-^^^ . (4.8) 
By the same method as before, we obtain 

/m /■ k 

p(^,^) Esg^(^) E (n 

(4.9) 

Wick's theorem for complex bosons (Theorem IA.4|) now gives 



— 1 / ^^r c. — 1 



k 

*rjCT(T(r))"araT(r) 



(4.10) 

Once again we define the permutation vr = cr o r and change variables, using now 
sgn(cr) = sgn(r) sgn(7r). We again have 



k 

5^sgn(7r)n(^"^)v..M = det{{X'^)ij) = (det X)-V(/, J) (det X,.,., 

7rG5fe r=l 

(4.11) 

while now 

m k 

T-eSk ai,...,afc = l r=l reSfe 

_ (^_]^^#(odd cycles of r) ^_^^#(cycles of r) 

reSfe 

= i-l)Ks{s + l)---is + k-l) (4.12) 

by (13.91) since s = —m. This proves (12.21) for integer s < — 1, when X is a complex 
matrix whose hermitian part is positive-definite; we conclude by Proposition 12.181 as 
before. □ 
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4.2 Symmetric Cayley identity 



Proof of Theorem I2.2[ In this case we use real bosons. So assume that s = — m/2 
where m is a positive integer, and that X^^^ is a real symmetric positive-definite 
matrix; and let us introduce bosonic variables (f'f'^ for i = 1, . . . ,n and a = 1, . . . ,m. 
We then have 

(detX^y"^)-"^/^ _ Jvipe--^^^'""^^ . (4.13) 

The operator S^^J™ = ^{l + 6ij)d/dxij is exactly what is needed to bring down a factor 
— when acting on e~3'^^"^™'^. (To lighten the notation, let us henceforth write 

X instead of X^^^.) We therefore have 



/m ^ k 

J2 sgn(a) J2 (n 

. (4-14) 

Wick's theorem for real bosons (Theorem IA.3P applied to JVip Y[r=i fi^''^ f^j^^}) 
gives rise to two types of contractions: those that only pair i's with j's, and those 
that pair at least one i with another i (and hence also a j with another j). The 
pairings of the first class are given by a sum over permutations r, and yield 

k 

(det X)-/2 Yl n(^"')-^-.(.«)^-".« • (4.15) 

TG5fc r=l 



Changing variables as in the alternate proof of Theorem 12. H the sum over tt again 
yields (detX)~^e(J, J) (detX/cjc), while the sum over r yields 

(-l)#(°'i<i ^y^i^^ °f (_^)#{cycies of r) _ (_i)'=(2s)(2s + 1) • ■ ■ (2s + A; - 1) 

(4.16) 

by (13. 9p since s = —m/2. Inserting these results into (I4.14p . we obtain (12. 5p . 

Let us now show that for each pairing of the second class, the sum over a yields 
zero. By hypothesis, at least one j is paired with another j, say jp with jg. Now let 
TTpq G Sk be the permutation that interchanges p and q while leaving all other elements 
fixed. Then cr i— ?■ o" o jipg is a sign-reversing involution for the sum in question. 

Once again we argue that because the equality (12. 5p holds for infinitely many s 
and for X in a nonempty open set, it must hold symbolically in s. □ 



4.3 Antisymmetric Cayley identity 

Proof of Theorem 12.31 In this case we use "real" fermions. So assume that s is a 
positive integer, and let us introduce "real" Grassmann variables 6^"^ for i = 1, . . . , 2m 
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and a = 1, . . . , s. We can then write 

^pf ^antisym^s = J VO e^^^^-^^^""^ . (4.17) 

Now let J = {zi, . . . , i2k} with ii < . . . < i2k, so that 

0-e52j. r = l V(2r-l)«<T(2r) 



Applying this to f l4.17p . we obtain 

k 



pf(a,t--)(pfx— )^ = fve^ J2'sn{a) Y: (nC-.C))^^'''^"""^' 

(4.19) 

(To lighten the notation, let us henceforth write X instead of X'^'^*'^^™.) When X is an 
invertible real or complex matrix, Wick's theorem for "real" fermions (Theorem I A. 151) 
gives 



I J. J. *CT(2r-l) *(T(2r) / 



r=l 



1 

= *^P^"^)''2^ ^ sgn(r)]^(X )v(,(2,_i))i„(,(2r))<^arr(2r-i)i"r-(2r)i • 

■ rg52fc r=l 

(4.20) 

Once again we define the permutation n = a o t and change variables from (a, r) to 
(r, tt), using sgn((T) sgn(r) = sgn(7r). The sum over vr gives 

k 

^ J] sgn(7r)n(X"T),^^,^_^,,^(,^, = pf((X-T),,) = (pfX)-^(/)(pfX,.,.) 

' 7rG52fc r-=l 

(4.21) 

by the pfaffian version of Jacobi's identity [cf. (IA.2ip ]. The remaining sums give 



^ s k 

_L \^ \^ TTr^ - # (cycles of A/UMo) U r^c.^ 

2kk\ ll"ar-(2r-i)iar-{2r)l ~ ' l^-^^J 

where the latter sum runs over all perfect matchings M of 2k elements, and Mq is 
some fixed perfect matching [in our case (12)(34) ■ ■ ■ {2k — 1 2k)]; we observe that if 
M and Mq are thought of as edge sets on [2k], then M U Mq is the edge set of a 
graph in which each vertex has degree 2, and so is a union of cycles, showing that 
the right-hand side is well-defined. Let us now show that this latter sum equals 
s(s + 2)-- -(5 + 2^-2): 
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Let ei,...,efc be the edges of Mq. Then each matching M G M.2k induces a 
decomposition of the set {ei, . . . ,6^} into cycles, according to how those edges are 
traversed in M U Mq; in other words, M induces a permutation of [k]. Moreover, 
for each cycle C in this decomposition, there are 21'^'^^ ways of connecting up the 
vertices. Thus each permutation tt of [k] arises from HcGTr^''^'^^ ~ 2'^~#('=y'=^'^'^ °^ 
different matchings M. We therefore have 

^#(cycles of MUMo) _ 2^ (g/2) °^ '"^ (4 23a) 

= s(s + 2)---(s + 2A;-2) . (4.23b) 
Once again we invoke Proposition 12.181 to conclude^ □ 



5 Proofs of Cayley-type identities by Grassmann 
representation of det(^) 

In this section we give an alternate Grassmann-based approach to proving Cayley- 
type identities: now it is the differential operator det(9) that is represented as a 
fermionic Gaussian integral. This technique is in our opinion very powerful: not only 
does it give the slickest proofs of the three main Cayley-type identities for square 
matrices; it also gives the only direct algebraic/combinatorial proofs (thus far) of the 
(considerably more difficult) rectangular Cayley identities. 

The basic fact we will need is that an operator exp(a ■ d) generates translation by 
a. More precisely, let i? be a commutative ring containing the rationals; then for any 
polynomial P (2:1, ... , Zn) with coefficients in R and any constants G R, we 

have the translation formula 

exp (^a,-^^ P{z) = P{z + a) . (5.1) 

(Here exp is defined by Taylor series; note that all but finitely many terms will 
annihilate P.) Indeed, the identity (15.11) is nothing other than Taylor's theorem for 
polynomials P. In particular we will use (15. ip when the commutative ring R consists 
of the even elements of some Grassmann algebra. Moreover, in our applications the 
elements a, will be nilpotent, so that the Taylor series for exp(a ■ d) is in fact finite. 

We will also need a formula for the change of a determinant under a low-rank 
perturbation: see Appendix IB.2I 

Let us begin by explaining the general structure of all these proofs. We introduce 
a Grassmann integral representation for the differential operator det(9) and let it 

We thank Alex Scott for help in cleaning up our proof of (j4.23|) . 



31 



act on (detX)'^. After a change of variables in the Grassmann integral, we obtain 
the desired quantity (detX)*~^ [or its generalization for minors J, J] multiplied by a 
purely combinatorial factor that is independent of the matrix X. We then proceed 
to calculate this combinatorial factor, which turns out to be an explicit polynomial 
in s. 

Unfortunately, the "all-minors" versions of these proofs (i.e. those for |/| = | J| = 
k < n) are slightly more complicated than the "basic" versions (i.e. those for I = J = 
[n]). We have therefore structured our presentation so as to give the "basic" proof 
first, and then indicate the modifications needed to handle the "all- minors" case. 



5.1 Ordinary Cay ley identity 



Proof of Theorem 12. 1[ We introduce Grassmann variables r]i,r]i {I < i < n) 
and use the representation 

detid) = jVn{r^,n)^''^\ (5.2) 

where the subscript on V serves to remind us of the length of the vectors in question, 
and we have employed the shorthand notation 

^ d d 

[Indeed, (15.21) is simply a special case of the fermionic Gaussian integral flA.74p where 
the coefficient ring R is the ring Q[d] of polynomials in the differential operators 
d/dxij.] Applying (15. 2p to (detX)** where s is a positive integer and using the trans- 
lation formula fl5.1D. we obtaiirj 



det(a) (det Xy = j D„(r/, r/) det(X + r/r/^)^ (5.4a) 
= {deiXy IVn{r],f]) det{I + X-^f]7i^y . (5.4b) 



(We assume here that X is an invertible real or complex matrix.) Let us now change 
variables from (r/, f/) to [rj'^f]') = {r],X~^f]); we pick up a Jacobian det(X~^) = 
(detX)~^ and thus have 



det(9) (detX)^ = (detX)^-^ Vn{v',v') det(/ + r/V^)^ (5.5) 



We have assumed here that s is a positive integer, because we have proven (|5.1|) only for 
polynomials P. Alternatively, we could avoid this assumption by proving (|5.1|) also for more general 
functions (e.g. powers of polynomials) when all the at are nilpotent. 
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Formula flS.Sp expresses det{d) (detX)'' as the desired quantity (detX)* ^ times a 
purely combinatorial factor 



P(s,n) = JVniv.v)detiI + r]r]'y, (5.6) 

which we now proceed to calculate. The matrix I + fji]'^ is a rank-1 perturbation of 
the identity matrix; by Lemma [B. Ill we have 

det{I + f]r]^y = {l-f]'^7])-' (5.7a) 
= E(-1)'(7)(^V (5.7b) 



where 



Since 



n 

ffrj = ^ f]irii . (5. 

i=l 



V^{r^,f,){fr^Y = n\5,^n. (5.9) 



it follows that 

P(s,n) = (-l)"(^~^^n! (5.10a) 

= s{s + l)---{s + n-l) . (5.10b) 

This proves (12. ip when X is an invertible real or complex matrix and s is a positive 
integer; the general validity of the identity then follows from Proposition 

Let us now indicate the modifications needed to prove (12. 2p for a minor /, J, where 
/, J C \n\ with |/| = \ J\ = k. We begin by introducing a Grassmann representation 
for (\et{dij): 

where 

(n^^)j.jc = VhVh---r]i^-,Vjr.-k (5.12) 



^■^ Alternatively, one can observe (see the Remark after Theorem 12. ip that the Cayley identity 
p.ip is equivalent to 

det(9) det(X + Ay = s{s + 1) • • • (s + n - 1) det(X + A)"-^ 

for any fixed matrix A. This latter identity can be proven exactly as above, with replaced by 
{X + A)^^ , and it works whenever X + A is invertible. Taking, for instance, A = —cl, we prove the 
identity for all real or complex matrices X whose spectrum does not contain the point c. Putting 
together these identities for all c € C, we cover all matrices X. 
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and ii, . . . ,in-k (resp. ji, . . . , jn-k) are the elements of P (resp. J'^) in increasing 
order. Applying the translation formula (15. ip as before, we obtain 

det(a,j)(detX)^ = e{I,J) JVn{7],r]) (^l[m)^^j^det{X + r]r]''y (5.13a) 

= e(J,J)(detX)'^ [V^{V,V) {Hvv) ^ ^ det{I + X-'rjv'^y . 

(5.13b) 

Once again we change variables from (//,//) to {f]' ^fj') = {i],X^^fi), picking up a 
Jacobian det(X~^) = (detX)~^; and we use the identity (15. 7p . obtaining 

(det9,,) (detX)^ = e(/, J) (det Xy' J V,,iv,v) (n(^^)^),.,. E(-l)' (T) ^^^^^ 



i=0 

(5.14a) 

X 



^ _ p=i / e=o ^ ^ 

„ /n—k \ 



The rules of Grassmann integration constrain the integral to be zero unless i = k 
and (ri, . . . , Vn-k) = ■ ■ ■ ,ja(n-k)) for some permutation a e Sn-k- We therefore 

have 

(det9,,)(detX)^ = e(J,J)(detX)-i J] sgn(a) ( Hx.^.^^^, ] (-1)^ x 

y"^n(^, ^) f n %^iv ) ^' ( n ^^^^ ) (5.15a) 



vP=i / VieJ 



e(/, J) (det X)^-i s(s + 1) ■ ■ ■ (s + - 1) ^ sgn(a) J] X^^ 



J<t(p) 

cr€S„-k P=l 



(5.15b) 

e(/, J) (det Xy-^ s{s + l)---{s + k-l) det(X/cjc) . (5.15c) 



This proves (Q- □ 



5.2 Two useful lemmas for the all-minors case 

Let us now pause to abstract the type of reasoning that was just used in proving the 
all- minors identity (12.20 . as similar reasoning will be needed to prove the all-minors 
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versions of the symmetric and rectangular Cayley identities and the all-principal- 
minors versions of the antisymmetric identities. The reader who is interested mainly 
in the case I = J = [n] can skip this subsection on a first reading. 

The key results of this subsection will be a pair of general formulae: Lemma 15.11 
for "complex" fermions and Lemma [5.31 for "real" fermions. Important special cases 
of these formulae (which also have easier direct proofs) will be stated in Corollar- 
ies [5]2] and |5]11 respectively. These corollaries are, in fact, all we need to handle the 
all-minors symmetric Cayley identity (Section 15. 3p and the all-principal-minors anti- 
symmetric Cayley identity (Section l5.4p . However, the rectangular Cayley identities 
(Sections l5.5H5.9p will need the full strength of Lemmas 15.11 and 15.31 to handle the 
all-minors case. 

Let /, J C [n] with \I\ = \ J\ = and let f^, be the complementary subsets. We 
denote by ^i, . . . ,in-k (resp. ji, . . . ,jn-k) the elements of (resp. J'^) in increasing 
order. Let r]i,f]i (1 < i < n) be Grassmann variables as before. Then, for any n x n 
matrices A,B, we define 



In particular, {{Ari){Bri))^ ^ = 1. 

Now suppose that we have further sets of (real) Grassmann variables 6^"^ (1 < 
i < n, 1 < a < N) — the case = is also allowed — and suppose that /(//, f/, 9) 
is a polynomial in the scalar products f]^i], ff-0^°'\ r]'^9^°'^ and 6'^°''^'^ 6^l^\ Then a 
Grassmann integral of the form 



obviously takes the same value for all sets L C [n] of the same cardinality. 

Let us now show how a Grassmann integral involving {J\{Aff){Bri)) jc jc and /(?7, f/, 9) 
can be written as a determinant containing A and B multiplied by a purely combi- 
natorial factor: 

Lemma 5.1 Let I,J(^ [n] with \I\ = \ J\ = k, let A,B be n x n matrices, and let 
/(?7,f/,^) be a polynomial in the scalar products as specified above. Suppose further 
that the number N of additional fermion species is even. Then 



(n(^^)(5r/))^^^^ = {Ar^UBv)n■■■m^.-.iBr^) 



Jn—k 



(5.16) 




(5.17) 




(5.18) 



where L C [n] is any set of cardinality n — k. 
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Proof. The expansion of the product fl5.16p produces 



n—k 

(n(^^)(^^)),.,,. = E X{K..B,^,s,f^^^r^s, . (5.19) 

ri, . . . ,r„_fc G [n] P=l 
SI,... ,s„_k e [n] 

Note that, by nilpotency of Grassmann variables, the product on the right-hand 
side of (15.191) is nonvanishing only if the indices ri, . . . ,r„_fc are all distinct and the 
indices Si, . . . ,Sn-k are also all distinct. Let us now integrate (I5.19P against one of 
the monomials arising in the expansion of f{r],f],6). For each "site" i G [n], this 
monomial contains an even number of factors rji, fji or O^""^; and it must contain each 
one of the factors ol"^ (1 < a < A^) precisely once if the Grassmann integral over 6 is 
to be nonvanishing. Since is even, this means that at each site i we must either have 
both of the factors rji and f/j, or neither. In order to have a nonvanishing Grassmann 
integral over t] and f/, the former situation must occur at k sites and the latter at 
n — k sites; moreover, the latter n — k sites must correspond precisely to the factors 
fjrpVsp in (I5.19p . We can therefore assume that (ri, . . . , Vn-k) = (^^{i), • • • , ■So-(n-fc)) for 
some permutation a G Sn-k] the contributing terms in (I5.19P are then 

n—k 

ri, . . . , r„_j. G [n] o-e5„_fe p=l 
all distinct 

The factors -Bjp,r .i^ ^ can be reordered freely; reordering of the Grassmann factors 
'7r^_i(pj yields a prefactor sgn((T). We therefore obtain 

n—k 

E E ^S'^i'^) Il^ip^rpBj^^^^^rpVrpVrp ■ (5.21) 

ri,...,r„_fc G [n] '^^■Sr^-k P=l 
all distinct 



By the remarks around (I5.17p . the integral 

„ /n—k \ 

J P„(r/, r/, 6) VrpVrpj fiv, V, 6) (5.22) 

is independent of the choice of ri, . . . ,r„„fc (provided that they are all distinct) and 
hence can be pulled out. We are left with the factor 



n—k 

E E ^g^(^) n^^..^.%(.).'-.- (5.23) 

ri,...,r„_fe G [n] o-65„_fc p=l 
all distinct 
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We can now remove the restriction that the ri, . . . ,rn-k be all distinct, because the 
terms with two or more rj equal cancel out when we sum over permutations with the 
factor sgn((T). So fl5.23p equals 

n—k n—k 
ri,...,r„_fcG[n] o-e5„_fc p=l aeSn-k P=l 

= det[{AB^)icjc] . (5.24b) 

□ 

Remark. If is odd, the situation is different: at each site i the monomial 
coming from / must now provide exactly one of the factors rji and fji. This means 
that we can get a nonzero contribution only when n is even and k = n/2. For instance, 
suppose that = 1, n = 2, J"^ = {i} and J'^ = {j}, and that f{ri,f],9) = {ri^6){ff-6). 
Then 

jv2{r],f],9){j[{An){Br^))^^^^{ife){f9) = -AaB,, + A.^Bj, . (5.25) 
We shall not consider this situation further, as we shall not need it in the sequel. □ 



The following identity is what was used in the proof of the all-minors ordinary 
Cayley identity (12. 2p and will also be used in the proof of the all-minors symmetric 
Cayley identity fl23|l : 

Corollary 5.2 Let I,J^ [n] with \I\ = \ J\ = k, let A, B be n x n matrices, let M 
be an invertible n x n matrix, and let i be a nonnegative integer. Then 

jv^{r],r]) (^l[{Ar]){B7])') JfMr]Y = k\6e,,{detM) det[{AM-^ B^)icj.] . 

(5.26) 

First proof. The case M = / is an easy consequence of Lemma 15.11 The case 
of a general invertible matrix M can be reduced to the case M = / by the change 
of variables [fj^f]') = [Mr],-?]), which picks up a Jacobian det M and replaces B by 
BM-\ □ 

Here is a simple direct proof of Corollary 15.21 that avoids the combinatorial com- 
plexity of the full Lemma 15.11 and instead relies on standard facts from the theory of 
Grassmann-Berezin integration (Appendix Rlj^ : 

Said another way: The proof of Wick's theorem for "complex" fermions fTheorem lA.16|) requires 
combinatorial work similar in difficulty to that occurring in the proof of Lemma l5.ll But since Wick's 
theorem is a "standard" result, we can employ it without reproving it ab initio. 
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Second proof. It is easy to see that the Grassmann integral is nonvanishing 
only when £ = k. So in this case we can replace {ff- Mrf)^ in the integrand by 
k\ e'x^{ff- Mrf). The claim is then an immediate consequence of Wick's theorem for 
"complex" fermions (Theorem IA.16p . □ 



In order to handle the all-principal-minors versions of the antisymmetric identities, 
we shall need formulae analogous to Lemma 15.11 and Corollary 15.21 but for "real" 
rather than "complex" fermions. So let I C [2m] with |/| = 2fc, and let be the 
complementary subset. We denote by ii, . . . ,i2m-2k the elements of P in increasing 
order. Let 6i, . . . , 62m be "real" Grassmann variables. Then, for any 2m x 2m matrix 
C, we define 

(5.27) 



«2m-2fc 



(note that the order of factors 6 is crucial here because they anticommute). In par- 
ticular, iUico))^ = 1. 

Now suppose that we have N further sets of (real) Grassmann variables x!"'' 
(1 < z < 2m, 1 < a < A^) — the case = is also allowed — and suppose 
that f{0,x) is a polynomial in the symplectic products ^9J9, 9Jx^°'\ ^X^^^Jx^"^ ^nd 
X^^'^Jx^^K where the 2m x 2m matrix J is defined by 



J 



1 




-1 






1 




-1 



(5.28) 



Now consider a Grassmann integral of the form 



'D2m{0,X) 




f{0,x) 



(5.29) 



for some set L C [2m], where the product YlieL understood to be written from left 
to right in increasing order of the indices. Clearly \L\ must be even for the integral 
fl5.29p to be nonvanishing. But more is true: L must in fact be a union of pairs 
{2j — l,2j}, otherwise fl5.29p will again vanish@ To see this, it suffices to observe 
that all the variables 9 or with index 2j — 1 or 2j appear in pairs in f{9,x)j 
therefore, in order to have an integrand that is even in these variables, the set L must 
contain either both of 2j — 1 and 2j or neither. Let us call a set L C [2m] well-paired 
if it is a union of pairs {2j — l,2j}. Obviously the integral fl5.29p takes the same 
value for all well-paired sets L C [2m] of the same cardinality. 



This statement of course pertains to our specific choice of J, and would be modified in the 
obvious way if we had chosen a different convention for J. 
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Let us now show how a Grassmann integral involving (n(^^))j''= fi^yX) can 
be written as a pfaffian containing C multiplied by a purely combinatorial factor: 



Lemma 5.3 Let I C [2m] with \I\ = 2k, let C be a 2m x 2m matrix, and let f{6, x) 
be a polynomial in the symplectic products as specified above. Then 

T^2m{e,x){l[(ce)')^j{e,x) 

= pf[(CJCT),.,c] jv2U0,x) (jlo}j f{9,x) (5.30) 
where L C [2m] is any well-paired set of cardinality 2m — 2k. 
Proof. The expansion of the product f l5.27p produces 

(2m-2fc \ 
n C..,rA ■ ■ 9r,_,, . (5.31) 

By nilpotency of Grassmann variables, the product on the right-hand side of fl5.3ip 
is nonvanishing only if the indices ri, . . . ,r2m-2k are all distinct. Moreover, as noted 
above, the set R = {ri, . . . , r2m-2k} must be well-paired if the integral is to be nonvan- 
ishing. Let us denote by r^, . . . , r'2,^_2f^ the elements of R in increasing order. Then we 
can write (ri, . . . , r2m-2k) = (r^(i), • • • , r'^{2ni-2k)) for some permutation a G S2m-2k- 
The contributing terms in fl5.3ip are then 



^2m-2fc 

V \ T\ Ci r' \dr' ■■■Or> ■ (5.32) 

/ ^ / ^ I 11 P' t(p) / V(l) V(2m-2fe) ^ ^ 

R C [2m] o-eS2m_2fc \ P=l / 
|_R| = 2m - 2fc 
_R well-paired 

We can now reorder the factors 6 into increasing order, yielding a factor sgn(cT), so 
that (15.32P becomes 

(2m— 2k \ / \ 

n a.;(J n^o • (5-33) 
^ _ - _ . p=i / Vie/? / 

\R\ = 2m - 2fc 
H well-paired 

We now multiply by f{6, x) and integrate P2m(^5 x)- As noted previously, the integral 

j'D2m{e,x) {X{e}jf{e,x) (5.34) 
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is independent of the choice of R (provided that it is well-paired and of cardinality 
2m — 2k) and hence can be pulled out. In order to calculate the prefactor, let us 
substitute the defining structure of well-paired sets R by writing r2f^_i = 2jh — 1 and 
^2h — for suitable indices ji < . . . < jm-k in [fn] ■ We are left with 

2m-2k 

E E ^g^(^) n ■ (5-35) 

1 < il < . . . < < m creS2m-2k P=l 

Alternatively, we can sum over all distinct ji, . . . jm-k ^ [fn] by inserting a factor 
l/(m — A;)!o Remarkably, we can now remove the restriction that the ji, ■ ■ ■ ,jm-k 
be all distinct, because the terms with two or more ja equal cancel out when we sum 
over permutations with the factor sgnfa)!^ We can also reorder the factors C freely. 
Therefore f l5.35p equals 

2m-2fc 



{m — k)\ 

il, ■ ■ ■ , jm-fc e H cr&S2m-2k P-1 

^ m—k 



E> sgn(cr) TT Ci , r-i Ci , 



(rn — k)\ ^ ^ '<T-l(29-l)'^2g-l V-l(29)''^2q 

il, ■ ■ ■ , im-fc 6 M '^6<S2m-2fc 9=1 



m—k 



2-^-^{m-k)\ 

si, . . . , e [2m] a-e52m-2fe <?— 1 

tl, . . -^t^-k G [2m] 

= pf[(CJC^)/cjc] . (5.36) 

□ 



The following identity is what will be used in the proof of the all-principal-minors 
antisymmetric Cayley identity (12. 8p : 

Corollary 5.4 Let I C [2m] with \I\ = 2k, let C be a 2m x 2m matrix, let M he an 
invertible antisymmetric 2m x 2m matrix, and let i be a nonnegative integer. Then 

I (n(^^)),. i^^^Mey = k\ h^k (pf M) pf [(CM-TcT),ve] . (5.37) 



Such a factor might not exist in a general commutative ring R. But we can argue as foUows: 
First we prove the identity when i? = ]R; then we observe that since both sides of the identity 
are polynomials with integer coefficients in the matrix elements of C, the identity must hold as 
polynomials; therefore the identity holds when the matrix elements of C are specialized to arbitrary 
values in an arbitrary commutative ring R. 

If ja = jp with a 7^ /3, it suffices to use the involution exchanging r'2a-i with 7*2^ (or 
alternatively the involution exchanging with r2^). 
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First proof. The case M = J is an easy consequence of Lemma 15.31 For a 
general real antisymmetric matrix M, we can use the decomposition M = AJA^ 
(Lemma IB.17P and the change of variables 6' = A^O (which picks up a Jacobian 
det A = pf M) to reduce to the case M = J. 

Finally, if M is an antisymmetric matrix with coefficients in an arbitrary com- 
mutative ring R, we argue that both sides of (15.371) are polynomials with integer 
coefficients in the matrix elements of M; since they agree for all real values of those 
matrix elements (or even for real values in a nonempty open set), they must agree 
as polynomials] and this implies that they agree for arbitrary values of the matrix 
elements in an arbitrary commutative ring R. □ 

Here is a simple direct proof of Corollary 15.41 that avoids the combinatorial com- 
plexity of the full Lemma 15.31 and instead relies on standard facts from the theory of 
Grassmann-Berezin integration: 

Second proof. It is easy to see that the Grassmann integral is nonvanishing only for 
i = k. So in this case we can replace {\0^ MOY in the integrand by k\ exp(^^"'^M^). 
The claim is then an immediate consequence of Wick's theorem for "real" fermions 
(Theorem lAHD. □ 

5.3 Symmetric Cayley identity 

With Corollary 15.21 in hand, we can this time proceed directly to the proof of the 
all-minors identity. 

Proof of Theorem 12.21 Recall that the matrix 5^^™ is given by 



As before, we introduce Grassmann variables rii,r]i {1 < i < n) and use the represen- 
tation 




(5.38) 




') 



"^ki^iVj + VjVi) 



d 



(5.39) 



-i<j 



By the translation formula (15. ip . we have 



det(af7)/({x,,},<,) 



(5.40) 
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for an arbitrary polynomial /. We shall use this formula in the case f{{xij}i<j) 
det(X'^y™)'^ where s is a positive integer, so that 



(5.41) 

It is convenient to introduce the shorthand 



^^sym^) trans ^ ^sym ^ 1(^^T_^-1 



(5.42) 



for the argument of det. 

Let us now change variables from (77, f/) to [t]', fj') = (77, (X^^^™)"^?]), with Jacobian 
(det X*^^™)"^. Dropping primes from the new variables, we observe that the expression 
for the translated matrix can be written as 



SO that 

det(X^y")*=^'^'^^ = (det X^y°^) det \l + |(r/r/^ - {X'^'^)-^r]f]'^X'y'^) 
Applying Corollary IB. 121 to the rightmost determinant yields 



det 



(5.43) 



(5.44) 



(5.45) 



so that we are left with the Grassmann-integral expression 
det(9^7) (detX^y'")^ 



(detX^y-)-^(/,J) /l)„(77,77) (n(^^'"^)(^)),. 



^7] 7]) 



Now insert the expansion 



in'"'!)-'- 



E(4)'(-f) 



(5.46) 



(5.47) 



into ( I5.46P and use Corollary 15.21 we obtain 

det(a?7) (det X'^'^y = (det X'^'^y-^ e(/, J) (det X^l%) k\ (- 



l\k 
2) 



-2s 

k 



(5.48) 



This proves (12. 5p when X**^™ is an invertible real or complex symmetric n x n matrix 
and s is a positive integer; the general validity of the identity then follows from 
Proposition [2TT81 □ 
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Remark. A slight variant of this proof employs the factorization X^^"^ = AA^ 
when X**^™ is a real symmetric positive-definite n x n matrix (Lemma IB. 160 : we 
then use the change of variables {rj'jfj') = {A~^ri, A~^fi), with Jacobian (det = 
(det X'^y™)^^. This slightly shortens the calculations, but has the disadvantage that 
the proof is no longer purely algebraic / combinatorial, because it invokes a decompo- 
sition that is valid for real symmetric positive-definite matrices in order to define the 
needed change of variables (which involves A and not just X^^™). We therefore prefer 
to avoid matrix decompositions wherever we can (which is unfortunately not always). 
See also the Remark at the end of Section 15. 4| and the discussion in Appendix IB. 31 



5.4 Antisymmetric Cayley identity 

With Corollary 15.41 in hand, we can again proceed directly to the proof of the 
all-(principal-)minors identity. 

Proof of Theorem 12.31 The 2m x 2m matrix ^^ntisym jg given by 

if ^ = J 



(5.49) 



d/dxij if z < j 
— d/dxji if i > j 

We introduce "real" Grassmann variables 6i {1 < i < 2m) and use the representation 

d 



pf{d^f-y^) = e{i) jvud) (n^),.^^P 

By the translation formula (15. ip . 



i<j 



(5.50) 



(5.51) 



for an arbitrary polynomial /. We shall use this formula in the case f{{xij}i<ij) = 
(pf X^°*''^^™)'^ where s is a positive even integer. It is convenient to introduce the 
shorthand 

^ji^antisym^ trans j^antisym _j_ QQT^ 52) 

for the argument of pf . 

Suppose now that X'^"*"^^™ is a 2m x 2m real antisymmetric matrix of rank 2m. 
Then Lemma IB. 171 guarantees that we can find a matrix A G GL{2m) such that 
j^antisym ^ AJA^, where 



J 



1 




-1 






1 




-1 



(5.53) 
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is the standard 2m x 2m symplectic form (note that with this convention pf J = +1 
for any m). We have 

(^antisym^trans _ A{J + 99^ A-^) (5.54) 

and 

pf^-5j.antisym^trans _ (pf X^'^"'^^'") pf ( J + ^-^^^^A"^) . (5.55) 

Now change variables from 9 to 9' = A^^9, with Jacobian (detA)^^ = (pfX)~^. 
Dropping primes, we are left with 

pf^^antisym^ ^antisym) s ^ ^^j^ J ^^^^g^ (n^^^)),^ P^i^+^^^T ' 

(5.56) 

We can now write 

pf(J + ^^^)* = det{J + 99^y/^ = det{I - J99"^y/^ (5.57) 

since pf J = det J = +1 and = — J@ The matrix / — J99^ is a rank-1 perturba- 
tion of the identity matrix; applying Lemma [B. Ill with vectors u = —J9 and v = 9, 
we obtain 

det{I-J99^) = {l-9^J9)-\ (5.58) 
Using Corollary 15. 4[ we have 

= (-2)^=^1 (^"^/^^ (pf J) pi[iAJ-^A^M 

= s(s + 2) ■ ■ • (s + 2m - 2) pf (X^^^f^^^) (5.59b) 

since pf J = +1 and J^"^ = J. This proves ( 12. Sp when x^'^*''^^™ is a real antisymmetric 
matrix of rank 2m and s is a positive even integer; the general validity of the identity 
then follows from Proposition 12.181 □ 

Remark. Although the statement of Theorem l2.3l is purely algebraic/combinatorial, 
the foregoing proof is unfortunately not purely algebraic/combinatorial, as it invokes 
the decomposition X'^"*''^^™ = AJA"^ that is valid for real full-rank antisymmetric ma- 
trices in order to define the needed change of variables (which involves A and not just 

28 The equality pf(J + 69^) = dct(J + e9^y/'^ follows from the general fact (pf M)^ = detM 
together with the observation that both pf(J+6'0''") and det{J + 96'^) are elements of the Grassmann 
algebra with constant term 1, and that an element of the Grassmann algebra with constant term 1 
has a unique square root in the Grassmann algebra with this property. 
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j^antisym-j_ ^ similar decomposition will be invoked in the proofs of the rectangular 
Cayley identities (Sections 15. 5[ 15. 6[ 15.71 and 15. 9p . This contrasts with the proofs of 
the ordinary and symmetric Cayley identities (Sections 15.11 and 15. 3p . where we were 
able to define the needed change of variables in terms of the original matrix X or 
^sym rjj^g matrix factorization lemmas needed for the former proofs are collected 
and proven in Appendix IB. 31 



5.5 Two- matrix rectangular Cayley identity 



In the remaining subsections of this section we shall prove the various rectangular 
Cayley identities that were stated in Section 12.21 It is convenient to begin with the 
two-matrix rectangular Cayley identity (Theorem 12. 6p . whose proof is somewhat less 
intricate than that of the corresponding one-matrix identities ( Theorems I2.7l and l2.8p . 
Indeed, the one-matrix rectangular symmetric and antisymmetric Cayley identities 
are related to the two-matrix identity in roughly the same way as the symmetric and 
antisymmetric Cayley identities are related to the ordinary one. 

Here we will need the full strength of Lemma 15.11 to handle the all-minors case. 

Proof of Theorem 12.61 We begin once again by representing the differential 
operator as a Grassmann integral: exploiting Corollary IB. 61 we have 



det{dxd' 



det 



Orn 






In 



Here (1 < i < m) and rjj,7]j (1 < j < "n.) are Grassmann variables, and the 

subscripts on V serve to remind us of the length of each vector; shorthand notations 



for index summations are understood, e.g. ijj dx'H = YlT=iYTj=i^iVjd / dxij. For a 



dei[{dxdl)ij] = det 





J\ = 


0^ 


dx 




In 



fc, we have (writing L = {m+1. 



IVJL.JVJL 



m+n}) 
(5.61a) 



= e(/,J) yP„(V^,^)I)„(r7,r7)(n^V^),.^/"'^^"''"'^^"'^'- 

(5.61b) 

Applying the translation formula (15. ip to the whole set of variables {xij^yij} pro- 
duces 

det[{dxdl)jj]f{X,Y) = e(/,J) ^©^(tA, V^) I)„(r7, r^) ( J] V^^)^^ ^^Z"" /(X+^r^^, 

(5.62) 

for an arbitrary polynomial /. We shall use this formula in the case /(X, Y) = 
det{XY^y where s is a positive integer. It is convenient to introduce the shorthands 



^trans 
y^trans 



= x + iJT]^ 

= Y + ijjf]^ 



(5.63a) 
(5.63b) 
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for the arguments of /. 

Suppose now that X and Y are real mxn matrices of rank m that are sufficiently 
close to the matrix /^^ defined by 

r 1 if i = ?■ 

(Imnh = r . ^ . (5-64) 

[Note that = {Im,Omx{n-m)) whcu m < n, and = otherwise. Henceforth 
we shall drop the subscripts mn on /^n to lighten the notation.] Then by Lemma [B. 181 
we can find matrices P, P' G GL{m) and Q G GL{n) such that X = PIQ and 
Y = P'Tq-^. We have 

det{XY^) = det{P7QQ-^PP'^) = det{PP'^) = det(P) det(P') (5.65) 

and 

(5.66) 

Let us now change variables from {ip, ip, rj, ff) to (■?/'', ■?/'', 77', fj') = {P'^^ip, P~^ip, Q^'^rj, Qfj), 
with Jacobian (det P)~^(det P')~^ = det(Xy"'")~^. In the new variables we have 
(dropping now the primes from the notation) 

^trans^ytrans^T _ P {f + ^jjT]'^ ) (P - f]^/j^ ) P'"^ ^ (5.67) 

and the translated determinant is given by 

det[(X*''^"^)(F*''^"^)^] = det(Xr'^) det[(/ + ^/'?7^)(/^ - #^)] , (5.68) 

so that 

det[{dxd^)ij] det{XY^y 

= e{I,J) detiXY^-' jvUi^,i,)n,{r],r])(^l[iPi>){P'ij)')^^^^ x 

e*'^" det[(J + iJv'^){P - #^)]^ . (5.69) 
Let us now split the vectors t] and fj as 

) (5.70a) 

{f]l,...,f]n) = (Ai,...,Am,Xl,---,Xn-m) (5.70b) 

so that 

(/ + V^r7T)(J^-#T) = + ^A^ - AV^T + c^V^T (5.71) 

with c = A"'"A + x^X- This matrix has the form of a low-rank perturbation + 
Z]a=i""^I' with vectors {t;^} given by 



a 






1 




A + c?/) 


2 


A 





(5.72) 
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By Lemma IB.llI we can write the needed determinant as the determinant of a 2 x 2 
matrix: 



, , 1 + A^V^ c + A^A , 
det-M ^ \ 5.73b 



= [(1 - ^Ij'm + A V) - {x'x){m-' , (5.73c) 

where the second equahty is a row operation (row 1 — t- row 1 + c row 2). We therefore 
have 

dei[{dxdl)ij] dei{XY^y = e{I,J) dei{XY^y-^ j V„{x,x) e^""" x 



(5.74) 



Note in particular that the integrand depends on if), if), A, A only via scalar products. 
This allows us to apply Lemma ET] to the integral over ■?/', ■?/', A, A; using also the fact 
that P(P')^ = XY"^, we obtain 

<lei[{dxdl)ij] det(XF^)" = e(/, J) dei^XY^Y'^ dei[{XY^)icjc\ x 
j V„{x, x) e*"^ j ^) P^(A, A) e^"^ \{ ^-aV-a^ 

[(1 - ^^A)(l + F^) - {x^xWn-' ■ (5.75) 



X 



This formula expresses det[(9xc?y)/j] det(XF'^)* as the desired quantity det(XF^)'* ^ : 
e(/, J) (lQi[[XY'^)icjc\ multiplied by the purely combinatorial factor 



P"{s,m,n,k) = yP„(^,^)P„(A,A)P„.„(x,x) X 



^a=A:+l 



which we now proceed to calculate. 

First note that the factor H^fe+i i^ai^a forces the Taylor expansion of the square 
bracket to contain no variables ipa, '4' a with k + 1 < a < m. We can therefore drop 
the factor IHLfc+i V^aV'a, forget about the variables V'a, "^a with + 1 < a < m, and 
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consider ip^ip henceforth as vectors of length k. Let us also rename the vectors A, A 
as 

(Ai,...,Am) = (Ai, . . . , Afc, /ii, . . . , /im-fc) (5.77a) 
(Ai,...,Am) = (Ai, . . . , Afc, /ii, . . . , /im-fc) (5.77b) 
Then the quantity (15.761) is equivalent to 

where scalar products involving ip, ijj, A, A are understood as referring only to the first 
k variables. Integration over the variables fi, ft is trivial and produces just a factor 1. 
So we are left with 



Note that pt^rect^^^ ^-j (^^gpg^ds on n and m only via the combination n — m. 
The binomial expansion of the integrand in f l5.79p yields 

[(i-^^A)(i+A^^)-(x^x)(v^^v^r^ = E 'm-¥m+>^^n-'-'[-{x'x){¥% 



h=0 



(5.80) 

For each fixed value of integration over variables Xi X gives 



I)„_„(X, X) e^^^ {-X^xf = h J ■ ^^-^^^ 

Now the factor [ijp^il:)^ will contain h pairs of variables, which can be chosen in h\ 
ways (counting reorderings) ; then in the term /(?/', Tp^ A, A) = [{l — 'ip^\){l + A^ip)]~'^~^ 
we will have to use only the other k — h pairs of variables for both ip, ip and A, A, and 
with a reasoning as in Lemma 15.11 we can choose them to be the first k — h indices, 
so that 

"v.ii^, ij) V,i\, A) e^"^ (V^^^) V(^, ^, A, A) 

^"j h\ I Vk-h{^P, ij) Vk-h{\ A) e^''^ /(V^, ^, A, A) , (5.82) 



and hence 



'Pfc_,(V', ^) Pfc-/^(A, A) e^"" [(1 - ^ta)(i + A^V')]-^-'^ . (5.83) 
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Let us now make the change of variables — )■ Aj, Aj — ?• —ipi (whose Jacobian is 
cancelled by the reordering in the measure of integration); we have 



Vk-hi^, ij) Vk-H{\ A) e^"^ [(1 - ^'X){1 + A 



-s—h 



^) Vk-H{\ A) e-^ [(1 + + A'^'A)] 



k-h 

Collecting all the factors, we have 

oo 
h=0 



Vk_h{^, ij) Vk-niX, A) [(1 + V)(l + A)] 



T\\i-s-h 



(5.84a) 
(5.84b) 
(5.84c) 

(5.84d) 

(5.84e) 



n — m 



h 



I n m \ 2 / 1 \k-h 



ik\) 

{kin-i 



h=0 

oo 

(fc!)^ 
(fc!)^ 

k-l 

Ylis + j){s + n - m + j) , 



-s — h 
k-h 



—s — h 
k-h 



1 2 



(A;-/i)! 



s + A; - 1 
k-h 



s + k-1 
k 



h=0 



k-h 



s + k — 1\ /s + k + n — m — 1 



k J \ k 



(5.85a) 
(5.85b) 

(5.85c) 

(5.85d) 

(5.85e) 

(5.85f) 



j=0 



where the sum over h was performed using the Chu-Vandermonde convolution (Lemma lB.2p 

This proves Theorem 12.61 when X and Y are real m x n matrices of rank m lying 
in a sufficiently small neighborhood of Imn, and s is a positive integer. The general 
validity of the identity then follows from Proposition 12.181 □ 
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5.6 One-matrix rectangular symmetric Cayley identity 

The proof of the one-matrix rectangular symmetric Cayley identity is extremely 
similar to that of the two-matrix identity, but is slightly more complicated because 
it involves a perturbation of rank 4 rather than rank 2 [compare f l5.97p with fl5.72p ]. 
Luckily, the resulting 4x4 determinant turns out to be the square of a quantity 
involving only a 2 x 2 determinant [cf. fl5.98p / f l5.99p ]. Once again, we will need the 
full strength of Lemma [5. II to handle the all- minors case. 

Proof of Theorem 12.71 We begin once again by representing the differential 
operator as a Grassmann integral: exploiting Corollary IB. 6t we have 



det(aa^) = det 



0™ 


d ~ 




In , 



Here {1 < i < m) and rij,fij (1 < J < n) are Grassmann variables, and we use 

the same conventions as in the preceding subsection. For a general minor /, J C [m] 
with |/| = I J| = fc, we have (writing L = {m + 1, . . . ,m + n}) 



det[{dd')ij] = det 





d 




In 



IUL,JUL 



(5.87a) 
(5.87b) 

Applying the translation formula (15. ip to the whole set of variables {xij} produces 

det[{dd'^)jj]f{X) = e{I,J) y"l)^(7/>,^)P„(r/,r/)(j]7/;7/>)^^^/"V(^+^^^+^^^) 

(5.88) 

for an arbitrary polynomial /. We shall use this formula in the case f{X) = 
det{XX"^y where s is a positive integer. It is convenient to introduce the short- 
hand 

jS^trans ^ X + l/jT]^ + l/jf]^ (5.89) 

for the argument of /. 

Suppose now that X is a real m x n matrix of rank m that is sufficiently close 
to the matrix I^n defined in (15.640 . Then by Lemma IB. 191 we can find matrices 
P G GL{m) and Q G 0{n) such that X = PIQ [we drop the subscripts mn on I^n 
to lighten the notation]. We have 

det(XX'^) = det (P/gg^J^P^) = det(PP^) = det(P)2 (5.90) 

and 

j5^trans ^ X + l/jT]'^ + ^f]^ = P[T + p-^I^T]^ + lPf]^)Q^]Q . (5.91) 
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Let us now change variables from {ip, ip, r], r]) to ip', v'j v') = {P P Qv^ Qv)y 
with Jacobian (detP)-^ = det(XX'^)-^ In the new variables we have (dropping now 
the primes from the notation) 

and the translated determinant is given by 

^g^[^^trans)(^trans)T] ^ aet{XX^) det[(/ + t/jT]'' + ^f){P - T]^^ - #^)] , (5.93) 

SO that 

det[{dd"^)ij] det(XX^)" 

= e(/,J) det(XXT)^-iy"p^(^,V')2^„(r/,r/) (n(^^)(^^)),,^, x 

(5.94) 



e"^'"" det[(/ + ^r/^ + - vi^^ - #'^)]" 

Let us now split the vectors rj and f] as 

iVly ■ ■ ■ :Vn) = (Al, . . . , Am, Xl; • • • ; Xn-m) 
(r/1, . . . ,//„) = (Ai,...,Am,Xl,...,Xn-m) 



(5.95a) 
(5.95b) 



so that 



{I+^r]'+^f]^){P-r]i;'''-f]i;'') = Im+i^X'+tpX''-\^'''-\ij''+cipij'''-cipij''' (5.96) 

with c = A"'"A + x^X- This matrix has the form of a low-rank perturbation Im + 
ZlLi^a^Q' with vectors {ua}, {va} given by 



a 






1 




A 


+ ci/j 


2 




A 


— cip 


3 


A 






4 


A 




~^P 



(5.97) 



By Lemma IB. Ill we can write the needed determinant as the determinant of a 4 x 4 
matrix; after a few row and column manipulations we can write 

det[(/ + ijT]'^ + ^pf]^)iP - - #^)] 
/ 

= det-^ 



A 


x^x 
-x^x 




^^ip 


—ip'^ip 



AT 



[detA-(^T^)(x^x)]" 



(5.98a) 
(5.98b) 
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where 

1 + X-^i) AT*. 



and hence 

del A = 1 + X^iIj + X^iIj+{X^iIj){X^^) + {X'^iIj){iIj^X). (5.100) 
We therefore have 

det[{dd^)jj] det{XXy = 6(1, J) det{XX^y Jv^_^{x,x) e^"'' x 

J VmiiJ,i')V^{X,X) e^"^ (n(^V5) W))^^ [detA- -(5.101) 

Note in particular that e^^'*' and detA — {'4'^'ip){x^x) depend on ipjip, X,X only via 
scalar products. This allows us to apply Lemma [5A] to the integral over X,X] 
using also the fact that PP^ = XX^, we obtain 

det[{dd'^)ij] det(XX^)' = e(/, J) det(XX^)'-^ det[(XX^)/cjc] x 

j V^-UX, X) e^"^ j VUi^, i:) P„(A, A) e^"^ ( f[ V^.V^J [det A - (^^V^)(x^x)]" 



(5.102) 



This formula expresses det[{dd^)jj] det(XX'^)'' as the desired quantity det(XX'^)'^ ^ : 
e(/, J) det[{X X'^)jcjc] multiplied by the purely combinatorial factor 



psymrect^^^^^^^^) = j V^i^ , i>) V^X, X) V^.^X, X) X 

eAWx I pj [detA-{ij^ij){x^x)]-'' , (5.103) 

\a=k+l ) 

which we now proceed to calculate. 

First note that the factor H^fc+i 'i'a^a forces the Taylor expansion of [det A — 
(?/'■?/') (xx)]~^^ to contain no variables ^/'a, ^/'^ with A; + 1 < a < m. We can therefore 
drop the factor Hl^fc+i '^a^a-, forget about the variables t/'q, with + 1 < a < m, 
and consider -i/^, henceforth as vectors of length k. Let us also rename the vectors 
A, A as 

(Ai,...,Am) = (Ai, . . . , Afe,/ii, . . . ,/im-fc) (5.104a) 
(Ai,...,Am) = (Ai,...,Afc,yUi,...,/im-fc) (5.104b) 

Then the quantity fl5.103p is equivalent to 

P^y-'--*(s,m,n,fc) = jVk{tP,tP)Vk{X,X)V^^k{f^,fi)Vr,.m{x,x) x 
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where det A has the same expression as in (15.1001) but scalar products involving 
ijj,ip,X,X are understood as referring only to the first k variables. Integration over 
the variables /x, /2 is trivial and produces just a factor 1. So we are left with 

psymrect^^^ m, 71, k) = [ V,{ij, ^) Pfc(A, A) x) e'"'+^"^ [det A-{i,^ij){x^x)]-'- 



(5.106) 

Note that psy™i'ect(^g^ depends on n and m only via the combination n — m. 

The multinomial expansion of the integrand in (I5.106P (except for e-^ ^) is 



E 



-2s \ /t, + ...+t5\ ,T,(-iri ^ 



r > 
ti,...,t5 > 

(A'^A)'^(X^X)*'(A''V')*'+*'(A''V')*''"*1(A^V')(^''A)]**(V'''V')*', (5-107) 



The integration 



Vn-m{x.x)e^''^{x^xY = . (5.108) 

{n — m — i)\ 



is trivial. So we are left with two integrations over complex Grassmann vectors of 
length k. This is performed through a lemma that we shall prove at the end of this 
subsection: 

Lemma 5.5 For all integers a, a', b, b', c, d > 0, we have 

V^iij, ij) P„(A, A) (A^A)'^ (^^V^)"' (A^^)" {X'^Y {\'i,r {iJ'xY 

n ^ ,2 



Saa' Sbb' Sec' Sa+b+c,n { ^ ^] (a^-bld) ■ (5.109) 



Combining formula (15.1071) (after integration of x) the statement of the 
lemma, we are left with three independent summations; we choose the summation 
indices to be ti, and ts, which we relabel as h, I and j, respectively. The remaining 
indices are given by r = j, t2 = h, t4 = m— j — h — I. The resulting expression is 

p-«(.,™,n,..) = («y^(-f ) j: ("7) {'i') . 

j, h,l > 
j + h + l < k 

(5.110) 

Renaming —2s — k = a and n — m = b highlights the fact that there is no direct 
dependence of the summands on k. From Lemma [B.4I we have 



.hJ \ j J \ h J V k 

j,h,l>0 \ / \ / \ 

j + h + l < k 
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which finally gives 



-2s 
k 



-{2s + n — m — 1) 
k 



(klf 

k-l 

Y[{2s + j){2s + n - m - 1 + j). 

j=0 



(5.112a) 
(5.112b) 



This proves fl2.19p when X is a real m x n matrix of rank m lying in a sufficiently 
small neighborhood of Imn, and s is a positive integer; the general validity of the 
identity then follows from Proposition 12.181 □ 



Proof of Lemma 15.51 Let us rewrite fl5.109p by forming an exponential generating 
function: that is, we multiply both sides by «"(«')" /3*(/3')'' 'j'^i'j'Y /[a\a'\b\b'\c\c'\] and 
sum over a, a', b, b', c, c' > 0. So f l5.109p is equivalent to 



P„(^,V^)P„(A,A) e"^^+"''^ '^+^^^+'^'^^+^^^+^'^^ = (««' + + 77') 

(5.113) 

and this is the formula that we shall prove. Note first that the measure can be 
rewritten as 



Vn{ip, ip) Vn{X, A) = dipndipndXndXn " " " dipidtpidXidX^ 



(5.114) 



with no minus signs. So let us assemble ip, ip, A, A into a single Grassmann vector 6 of 
length 4n, with ip^ = ^ik, i^k = ^ik-i, X^ = 04k_2 and A^ = O^k-s] then the measure 
becomes T'4„(6'). In the exponential we have an antisymmetric bilinear form ^9Q9, 
in which Q is block- diagonal with n identical 4x4 blocks that we call Q': 



Q' 



/ a 7 /3 \ 

-a -i 

-7 -/?' a' 

\-/3 i -a' / 



(5.115) 



So the result of the integration is pf Q = (pf Q'Y- And indeed, pf Q' = aa'+(3(3'+'j'j', 
as was to be proven. □ 



5.7 One-matrix rectangular antisymmetric Cayley identity 

Let us now prove the one-matrix rectangular antisymmetric Cayley identity (The- 
orem [521) • Again, this proof is extremely similar to that of the two-matrix and the 
symmetric one-matrix rectangular identities. As always with pfaffians, we will deal 
only with principal minors, for which we will need the full strength of Lemma 15.31 
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Proof of Theorem 12.81 We begin once again by representing the differential 
operator as a Grassmann integral: exploiting Corollary IB.9t we have 



pf(aja^) 



02m 


d 


-6' 


J2n 



;-i)"pf 



(5.116a) 
(5.116b) 



where (1 < ^ < 2m) and rjj {1 < j < 2n) are "real" Grassmann variables, and 
we use the same conventions as in the preceding subsection. For a general even- 
dimensional principal minor I C [2m] with |/| = 2k, we similarly have (writing 
L = {m + 1, . . . , m + n}) 



pi[{dJd^)n] 



-l)'pf 



02m 


d 




J2n 



lULJUL 



(5.117a) 



= (-l)'e(/)yP2m(V^)2^2n(r/)(nV^),,e^''"^''"''^"''- (5-117b) 
Applying the translation formula (15. ip to the whole set of variables {xij} produces 

pi[idJd^)n]fiX) = y"l)2m(V^)2^2n(r])(n^),.e^''""'/(^ + V^^^) 

(5.118) 

for an arbitrary polynomial /. We shall use this formula in the case f{X) = 
pf{X JX"^)^ where s is a positive integer. It is convenient to introduce the short- 
hand 

^trans ^ X + ^7]"^ (5.119) 

for the argument of /. 

Suppose now that X is a real 2m x 2n matrix of rank 2m that is sufficiently close 
to the matrix l2m,2n defined in fl5.64p . Then by Lemma IB. 201 we can find matrices 
P G GL{2m) and Q G Sp{2n) such that X = Pl2m,2nQ- We recall the defining 
property of Sp{2n), namely that QJ2nQ^ = <^2n- We have 



so that 
and 



XJ2nX^ — Pl2m,2nQJ2nQ^l2m,2n-P^ ~ PJ2mP^ 

pf(XJ2„XT) = pf(PJ2mP^) = detP 



(5.120) 

(5.121) 
(5.122) 



^trans ^ X + ^JT]'^ = P[I + p-^t/jT]^ Q-^]Q . 

(we drop the subscripts on l2m,2n to lighten the notation). Let us now change variables 
from {ijj,ri) to {tp',ri') = (P'V, Q""^?/), with Jacobian (detP)"^ = pf(XJX'^)~^ In 
the new variables we have (dropping now the primes from the notation) 



X' 



P{I + ijr]^)Q 



(5.123) 
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and the translated pfaffian is given by 

= det(P) pf [(/ + ip7]^)J{P - r]^^)] 

= pf{XJX^)pf[{T+ij7]^)J{P -r]ij^)], (5.124) 

so that 

pf[{dJd^)n] pf{XJX^r 

eh^Jv pf [(J + ^ri^)J{P - ry^T)]s ^ (5^^25) 
Let us now spht the vector rj as 

(^/l, • • • ,^?2n) = (Al,...,A2m,Xl,---,X2{n-m)) , (5.126) 

so that 

{I + i,rj^)J2n{I^ -Vi^'^) = J2m-i'\'j2m- J2rrM'^ -Ciji^^ (5.127) 

with c = r]^J2nV = A""" J2mA + x'^'^2(n-m)X- The matrix ( I5.127P is manifestly antisym- 
metric, and has the form of a low-rank perturbation of J2m- Since pf J2m = 1 and the 
perturbation is purely Grassmannian, we can write 

pf (J - ipX^J^ - J\ii^ - ci)^^^) = det^/2(/ + JipX^J^ - Xip^ + cJip-^^) . (5.128) 

Here the argument of det is a matrix of the form l2m + Sa=i ^""^a ' with vectors {mq}, 
{v^} given by 



a 






1 




JX + cip 


2 


-A 





(5.129) 



By Lemma IB. Ill we can write the needed determinant as the determinant of a 2 x 2 
matrix, which, after a slight row manipulation, can be written as 



det-i 



1 — ip^X —x^Jx 



det(J + JipX^J^ - A^^ + cJipip'^) 

We therefore have 

pi[{d.Jd'')n] pi{XJXy = pf(XJXT)^-i(-l)'=e(/) j V2^n-m){x) e 



(5.130) 



-.x Jx 



X 



V2m{i^)V2^{X) e'^'"'' (n(^^)),. [(1 - ^^A)2 + (^^J^)(x^Jx)]" 



s/2 



(5.131) 
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In this expression we have both an ordinary scalar product (ip'^X) and symplectic 
scalar products {X^JX, ip'^Jip and x^Jx)- By a further change of variables A — )■ A' = 
— JA, we can reduce to symplectic products only (the Jacobian is 1). Dropping the 
primes, we have (also using J^JJ = J) 



pn{dJd'^)n] pf{X JX^y = pf{XJX^y-\~l)'e{I) Jv^^n^rrMe 



(5.132) 

Now the integral on the and A fields is of the form described in Lemma l5T3l Applying 
this lemma, and making use of fl5.120p . we have 

pf[{dJd^)n] pf(XJX'^)^ = e{I)pf{XJX^y~^pf[{XJX^)icj.] 

X (-1)'= j I?2(n-™)(X) e^^"^^ j V2m{X) 6^""^" 

^ (n^). . [{l-i^^JXf + {i,^Jij){x^Jx)V . (5.133) 

V-'--'- / {2k+l,...,2m} 

This formula expresses pf[{dJd'^)jj] pf(XJX'^)'^ as the desired quantity pf (XJX'^)*^^x 
e(/) pf [(X JX"'")/c/c] multiplied by the purely combinatorial factor 

P~(s,m,n,A;) = {-if j V^^n^raM) e^''" j V^^V^miX) e^"^^ 

^ (n^). Jl-^^JXf + i^-^J^yx^Jx)]-'", (5.134) 

V-'--'- / {2A:+l,...,2m} 

which we now proceed to calculate. 

Note, first of all, that the factor Y\^^2k+i'^a forces the Taylor expansion of [(1 — 
ijj'^JXy + {ijj^ J^lJ){x^ Jx)]~'^^'^ to contain no variables ipa with 2k + 1 < a < 2m, and 
thus also no variables Xa with 2k + 1 < a < 2m. The latter must all come from 
the trivial diagonal exponential ea^ and their integration can then be performed 
easily (it gives 1). We can therefore drop the factor na=2fc+i'^a' forget about the 
variables ipa and A^ with 2k + 1 < a < 2m, and consider ip and A henceforth as 
vectors of length 2k. We have 

P^^'-^^*(s,m,n,A;) = {-if Jv2in-m){x) e'^''"''' jv2kWV2kiX) x 

ei^"^^ [(1 - ij^JXy + i^P^J^PKx^Jx)]-'^' . (5.135) 



We see that p^srect^^^^ depends on n and m only via the combination n ~ m. 

In the expression (I5.135P we have symplectic scalar products of "real" Grass- 
mann fields. However, our summation lemmas (such as Lemma 15. 5p have heretofore 
been developed only for ordinary scalar products of "complex" Grassmann fields. 
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However, through a relabehng of the fields, we can rewrite fl5.135p in terms of or- 
dinary scalar products of "complex" fields of half the dimensions. More precisely, 
we relabel {ipi, . . . ,ip2k) {i'i,ipi, . . . ,ijjk,ipk), (Ai,...,A2fc) -> (Ai, Ai, . . . , A^, A^) 
and (xi, • • • , X2(n-m)) -> (Xi, Xi, • • • , Xn-m, Xn-m)- No slgus arlsc in the measure of 
integration, and we have the correspondences 



\^lj^Jip — > ijj'ilj (5.136a) 

iA^JA — ^ A^A (5.136b) 

^^JX — y tp^X + X^ij (5.136c) 

Ix^^X — > f'x (5.136d) 

which rewrites (15.1351) as 

P^^^^'^\s,m,n,k) = {-If jVk{iJ,i,)V,{X,X)V.^^^{x,x) 

X e^^^+*^^ [(1 - ^^A - X^ijf + A{^^i;){x^x)Y''^ ■ (5.137) 
We begin by getting rid of variables x ^"^^ Xi writing 

^Vr,_^{x, X) e^"^ [(1 - ^^A - X^i^f + 4(^T^)(x^x)]-^/' 

= 5^ ( ~M - ^^A - ATV')^'"'" ( a! . (5.138) 

a>0 V ^ V / 

So we have 

>asrect/'„ ^ „ l„\ / T\k I „T,\ ^ / \ 



P-^-^*(s,m,n,A;) = (-1)'= y Pfc(^, ^/.) P,(A, A) e-^ x 

^ r"^") (#T^)'^(1 - ^^A - AT^)-^-2a - ^"j a! . (5.139) 

a>0 V / \ / 

Next we expand fully the integrand in fl5.139p . yielding 

n — m\ ,f—s — 2a\ fh + h'\ 1 



a, a' > 
6, b' > 



b + b' J \ b I a'\ ^ 



Vk{^, ^) Vk{X, A) {^''mx''xr (^^A)^(V^TA)'' . (5.140) 
The fermionic integration here is a special case of Lemma 15. 5[ which thus gives 



a, a' > 
> 
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J2i-Ari-in-s - 2af-^ ("^ a\ (5.141) 



a=0 ^ / \ / 

where x- = x{x — 1) ■ ■ ■ {x — k + 1). Notice now that 

(-f )-(-s-2a)^^ = {-lTs{s+2) ■ ■ ■ {s+2k-2) (s+2a+l)(s+2a+3) ■ • ■ (s+2A;-l) . 

(5.142) 

The factors of the form s + 2j, namely s{s + 2) ■ ■ ■ (s + 2A; — 2), are independent of 
the summation variable a. What remains is 



2'^(s + 2a + l)(s + 2a + 3) ■ ■ ■ (s + 2fc - 1) 



a! 



2''k\J2 

a=0 

2^k\ ^ '~ 



' \k\ f n — m 



2 

, k — a 

a=0 

+ n — m + k 



k 

JJ(s + l + 2n-2m + 2j) (5.143) 



fc-i 



where the sum over a was a Chu-Vandermonde convolution (Lemma IB. 2p . Restoring 
the factors of the form s + 2j, we conclude that 



k-l 



P'^"'^\s,m,n,k) = Y[{s + 2j){s+l + 2n-2m + 2j) . (5.144) 

j=0 

This proves f l2.22p when X is a real 2m x 2n matrix of rank 2m lying in a sufficiently 
small neighborhood of l2m,2n, and s is a positive integer; the general validity of the 
identity then follows from Proposition 12.181 □ 



5.8 A lemma on Grassmann integrals of scalar products 

In the preceding proofs we have frequently had to evaluate Grassmann integrals 
over one or more sets of Grassmann variables, in which the integrand depends only 
on scalar products among those sets. A simple version of such a formula occurred 
already in (15. 9p . while more complicated versions arose in fl5.79p ff. and (I5.109p . [Also, 
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a version for real fermions arose in fl5.59p .] So far we have simply treated such integrals 
"by hand" . But since a much more complicated such integral will arise in the proof of 
the multi- matrix rectangular Cayley identity (Theorem I2.9p in the next subsection, it 
is worth stating now a general lemma that allows us to systematize such calculations. 
Indeed, we think that this lemma (Proposition 15.61 below) is of some independent 
interest. 

Let us start by rephrasing the simplest case (15. 9p in a suggestive way. Let n 
be a positive integer, and introduce Grassmann variables ipii^i (1 < ^ < '^). Let 
/(^) — YlT=o^kX^ ^ formal power series in one indeterminate. We then have by 
(EH) 

Jv4^,^)f{,lj^^) = n\an. (5.145) 
On the other hand, we also have trivially 

= nlttn (5.146) 

x=0 

and 

= nlttn. (5.147) 

x=0 

Hence 

(5.148) 

x=0 

Surprisingly enough, a similar formula exists for Grassmann integrals involving 
multiple sets of fermionic variables: 

Proposition 5.6 Leti andn he positive integers, and introduce Grassmann variables 
%pf,%pf (l<<y<i,l<i<n). Let f{X) be a formal power series in commuting 
indeterminates X = {xa/sYa^i^^i, and write d = {d/dxap). We then have 

y"p„(7^^^l)■■■p„(^^7A0/({rV}) = det{dr f{x)\^^, (5.149a) 

= /(9)(detX)"|^^o . (5.149b) 

Proof. We write 

f{x) = E n (5.150) 

where the sum runs over matrices = {uapYa p=i nonnegative integers. Note first 
that it suffices to prove (15.1491) for polynomials /, since all three expressions in f l5.149p 
have the property that only finitely many coefficients N contribute (namely, those 




fix) 
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matrices with all row and column sums equal to n). So it suffices to prove fl5.149p 
for all monomials := Yla i3=i ^I'jf ■ then it suffices to prove fl5.149p for the 
exponential generating function 

MX) = exptr(n^X) = E n ^"""'^"'f"' (5.151) 

N€^'x^ a,l3=l ^"l^' 

where fl = {u)a/3)ai3=i are indeterminates, since the value at the monomial can 
be obtained by extracting the coefficient [fi^] . Our goal is to prove that all the three 
expressions in (15.1491) . specialized to f{X) = 0q(X), equal (detf2)". 
For the Grassmann-integral expression we have 

V^{^\ ^1) ■ ■ ■ Vr,{tP', MKr'^^^}) 

= I V^{tlj\^') ■ • exp (j2 E V^r^a/3^f ) (5.152a) 

\i=l a,(5=l / 

= (detfi)". (5.152b) 

For the power series in derivative operators applied to a power of a determinant, a 
simple application of the translation formula (15. ip yields 



^{d){detxr\^^, = exp I ^ ) (detX)^ 



« 1 9Xal3 
\a,l3=l 



(5.153a) 



x=o 

= (det(X + ^]))"|;^^o (5.153b) 

= (detfi)". (5.153c) 

For the power of a determinant in derivative operators applied to a power series, it 
suffices to cite the proof (I5.153P and invoke transposition duality X 9 in the Weyl 
algebra; but for completeness let us give a direct proof. Use a fermionic representation 
of the differential operator 



^ %1 (5.154) 



det(9) = / Vi{r],r])exp j r]^ 
and apply again the translation formula (15. ip : 

det(9) [exp tr(^]^X)] = Jve{ri,f])expl ^ ^'^q^^'^] [exptr(fi^X)] (5.155a) 



J ^'£(r/,r/)exptr(^]'^(X + 7^7/^)) (5.155b) 
[exp tr(^]'^X)] JVi{r],f])expi^f]aUai3r]/3j (5.155c) 



(detl]) exptr(^]^X) . (5.155d) 
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Iterating n times, we get 

det(a)"0n(X)|^^o = (detfi)"exptr(fi^X)|^^Q = (detl^)'^. (5.156) 

□ 



5.9 Multi-matrix rectangular Cayley identity 

In this subsection we shall prove the multi-matrix rectangular Cayley identity 
(Theorem 12. 9p . which is the most difficult of the identities proven in this paper. The 
difficulty arises principally from the fact that the number of matrices appearing in 
the identity (which we call i) can be arbitrarily large, and we are required to provide 
a proof valid for all £. The proof nevertheless follows the basic pattern of the proofs 
of the other identities (notably the two- matrix rectangular identity), and divides 
naturally into two parts: 

(i) We represent the differential operator as a Grassmann integral, and after sev- 
eral manipulations we are able to express det[((9^^^ ■ ■ ■ d^^^)ij] det(X'^^) ■ ■ ■ X^^^Y 
as the desired quantity det(X(i) • ■ ■XW)"-ie(/, J) det[(X(i) ■ ■ •XW)/cjc] multi- 
plied by a purely combinatorial factor b{s) that is given as a Grassmann integral. 

e k 

(ii) We evaluate this Grassmann integral and prove that b{s) = Yl Yli^ ~^ — j)- 

a=lj=l 

The major new complications, as compared to the preceding proofs, come from the 
fact that our Grassmann integral involves i sets of fermionic fields, while in the 
preceding proofs we only had one or two, and the matrix arising from the application 
of the "low-rank perturbation lemma" (Lemma IB.lip is of size i x i, while in the 



preceding proofs we never had a matrix larger than 4x4. Despite these complications, 
step (i) follows closely the model established in the preceding proofs and is not much 
more difficult than them: the main novelty is that we need to use a variant of the low- 
rank perturbation lemma (Corollary IB.13P that is specially adapted to perturbations 
of the product form f l5.170p encountered here. The big trouble arises in step (ii): to 
evaluate the Grassmann integral involving scalar products among i sets of fermionic 
fields, we shall ffist rewrite it as a differential operator acting on a determinant (using 
Proposition 15. 6p and then exploit the fact that the matrix arising in this differential 
operator is upper Hessenberg (i.e. has zero entries below the ffist subdiagonal); the 
latter leads to some special combinatorial/algebraic computations (Theorem 15.71 and 
Corollary ES]). 

It is convenient to formally divide the statement of Theorem 12.91 into two parts, 
corresponding to steps (i) and (ii) above: 

Theorem 12.91 part (i) Fix integers £ > 1 and ni, . . . ,ni > and write n^+i = ni. 
Forl<a<i, let X(°) be an Ua X Ua+i matrix of indeterminates, and let d^"^ be the 
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corresponding matrix of partial derivatives. If I,J ^ [ni] with \I\ = \ J\ = k, then 



det[(5(i) ■ ■ ■ d^'^)ij] det(X(i) ■ ■ -XW)^ 

= bn,,...,n,,k{s) det(X(i) ■ ■ ■ e(/, J) det[(X(^) ■ ■ ■ xW),cj.] (5.157) 



where 



bni,...,ne;k{s) 



X exp 



a=2 



det{Ie + M)-' (5.158) 



and 



M, 



al3 



if a < 13 

_f^aT^a a = f3 + l 

otherwise 



(5.159) 



Theorem 12. 9L part (ii) We have 



a=l j=l 



(5.160) 



Let us now begin the proof of part (i) of Theorem 12.91 As before, we wiU need 
the full strength of Lemma [5. II to handle the all-minors case. 



Proof of Theorem 12.91 part (i). We can assume that > ni for 2 < a < i, 
smce otherwise det(X(i) ■ ■ ■ X^) is the zero polynomial. 

We begin, as usual, by representing the differential operator as a Grassmann 
integral, this time exploiting Lemma [B. 141 (rather than just Corollary IB.6|) . We have 



det(9(^)---a(')) 



Vr^,i^|J\^P')■■■n,^{^|J\^|J') 



X exp 



a=2 



.(5.161) 



Here ip"' , i/)"' (1 < a < i, 1 < i < ria) are Grassmann variables, and the subscripts 
on V serve to remind us of the length of each vector; shorthand notations for index 
summations are understood, e.g. ip"-^ d'^"'^ ip"''^^ = J2i=iJ2%V '^f'^j^^'^/'^^if ■ -^^^ ^ 
general minor J, J C [m] with |/| = \ J\ = k, the quantity det[{d^^^ ■ ■ ■ d^^'^)ij] has 

a representation like (I5.16ip but with an extra factor e(/, J) ( TT^^^'M 

V / i'=,j'= 

integrand. 
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It is convenient to make the change of variables — )■ as this makes the 

summand ijj^d^^^%jj^ analogous to the ijj°'~^d^"~^^ip"' arising for 2 < a < We shall 
exploit this structure by writing -0^"'"^ as a synonym for This change of variables 
introduces an overall factor (— = (—1)'^. 

Applying the translation formula f lS.ip to the whole set of variables {x[j^} produces 

det[(a«---9W),^]/(X«,...,xW) = e(/,J) 



X (-l)''exp 



a=2 



n fiX^lL • • • , ^£ns) (5.162) 



for an arbitrary polynomial /, where we have introduced the shorthand 

^Sns = X^-^-^P'^{r^'f (5.163) 

for the arguments of /. We shall use the formula fl5.162p in the case f(X^^\ ■ ■ ■ , X^^^) = 
det(XW---XW)^ where s is a positive integer. 

Suppose now that the X*^") are real matrices of rank mm(na,na+i) that are suf- 
ficiently close to the matrix luaUa+i defined in (15.641) . Then by Lemma [B.2 II we can 
find matrices Pa G GLijia) for 1 < a < £ + 1 such that X^") = Pq,/„„„„+i -Pq+i- We 
have X^^) ■ ■ ■ X*^^) = PiPl^^ since /„j„2-^„2„3 ■ ■ ■ In^m = In-^ as a consequence of the fact 
that ria > ni ioi 2 < a < I. Therefore 

det(X(^) ■ ■ ■ X(^)) = det(Pi) det(P£+i)-^ . (5.164) 

We also have 



^Sns ■ ■ -^Sns = n - K^r{r^^fP.^l\ P^ll (5.165) 

a=l 



where the product is taken from left (a = 1) to right (a = £). 

Let us now change variables from to defined by 



, PjV" for 2 < a < . 
r = <! - r (5.166) 



PZ-ii^^ for a = 1 



= p-i^" (5.167) 

The Jacobian is (det Pf ^)(det P^+i) = det(X(^) ■ ■ ■ X^)"^ [using flSAMj) ]. In the new 
variables we have (dropping now the primes from the notation) 

^Sns ■ ■ ■ ^£ns = Pi fn fc„n.+, " rir^'f]] , (5.168) 
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so again using fl5.164p we see that the translated determinant is given by 

detixilL ■ ■ ■ ^Sns) = det(X(i) ■ ■ ■ X^) (det M„,...,„J 
where the matrix depends only on the Grassmann variables: 



(5.169) 



ni,...,ne 



(5.170) 



Therefore 

det[(9«---9(^)),j] det(X«.--XW)^ 

= e(/,J) det(X«---XW)-i(-l)'= jVM\i,')...V^Xi^\i^') 



X 



a=2 



det(M„,, 



(5.171) 



To evaluate det(M„^_,,,_„^), we shall use Corollary IB.13t which is a variant of the 
low-rank perturbation lemma that is specially adapted to matrices of the form fl5.170p . 
Before doing so, it is convenient to split the vectors ip"" and as 

(^^...,CJ = (A°,...,A^^,Cr,...,CJ (5.172a) 

(<,...,<) = (A^...,A::^,Cr,...,CJ (5-172b) 

where ma := — rii (note in particular that ip^ = A^ and tp'^ = X^). We now apply 
Corollary IB. 131 with Xa = 4'°', Va = 4'°''^^ and e = — 1 to obtain 

det(M„,,...,„J = (detiV)-i (5.173) 

where the i x i matrix N is defined by 

i=i (5.174) 

5a/3-A"+^T^/3 a a>/3 

and m„ a := min m^. We now fall into the conditions for the application of 

a<7</3 ' 

Lemma |5TT] with = ■i/'^ = A^, = -j/;^ = A^ and = {A^, . . . , A^, A^, . . . , A^} (here the 
variables ( and ( just go for the ride), yielding 

det [(5(1) . . . ^W),^] det[X(i) ■ ■ ■ XW]^ 

= e(/, J) det[X(i) ■ • • X^'^Y^^ det[(X(i) ■ ■ ■ X^'^)jc^jc] (-1)'= 

X fv^.iWx')-.. (A^ AO v^, (c^ C') ■ ■ ■ (C, CO 



rti 

X I JJ A] A] I exp 

Kj=k+1 



^(A°^A" + C^C 



a=2 



det(X)- 



(5.175) 
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We have thus represented det[(c?^"^'' ■ ■ ■ d'^^'>)ij\ det(X*^-'^) ■ ■ ■ X'^^'^Y as the desired quan- 
tity det(X(i) ■ ■ ■XW)'^-ie(J, J) det[(X(i) ■ ■ ■XW)/cjc] multiphed by a purely combi- 



natorial factor h. 



'ni,...,n£;k 



's) that is given as a Grassmann integral: 



bn,,...,nAs) ■■= {-If jvr,,{\\\')...Vr,, (V , AO v^, (c^ c^) ■ ■ • v„,^ (CO CO 



ni 

X I J]^ AjA] I exp 

\j=k+l 

rni 



5^(A"^A" + C"V 



a=2 



det(X)"' . (5.176) 



In order to handle the factor Iljifc+i -^j-^]' convenient to further split the 
fields A and A as 



(A", . . . , A"J 
(A", . . . , A"J 



(A", . . . , A^, Xl , • • • , Xni-fc) 
(A", . . . , A^, Xl , • • • , Xni-fc) 



(5.177a) 
(5.177b) 



Notice now that the overall factor Hjifc+i ^j-^] ~ nr=r^ xlx] iii the integrand kills all 
monomials in the expansion of the rest of the integrand that contain any field or x^. 
Now, the factor (det X)~* in fl5.176p depends on the fields A, A, x, X, C? C o^ily through 
products of the forms {Af Af }i<^<„<£, {A^Uf }i<^<£, {x^xf }i</3<a<£, {XiXf }i</3<^ and 
{C"C^}2<a</3<£, while the exponential depends only on combinations A°A", x"x" aiid 
C°C" that are "charge- neutral" in each field separately. Therefore, the only monomials 
in the expansion of (det N)~^ that can contribute to the integral must also be charge- 
neutral in each field separately. But since no monomial containing any x^ can arise, 
it is impossible to make such a charge-neutral combination using any other x oi' X 
since all such terms are of the form x"X^ with [3 < a. In a similar way, the terms 
C°C^ with 2 < a < 13 < ^ cannot contribute. (The combinations C°C" do survive.) 
The integral 05.1761) will therefore be unchanged if we replace X by a new matrix N' 
in which all these "forbidden combinations" are set to zero: 



C^^C^ if a = 13 -I 

if a < /3 - 1 

5a/3-A°+iTA^ if a >/3 



(5.178) 



Note that the matrix N' is lower Hessenberg (i.e. has zero entries above the first 
superdiagonal) . Summarizing, we have 



-'ni,...,ni 



■AS 



V,{X\\')---Vu{X\\' 



X 2^„,-fc(xO XO ■ ■ • Vn,-k{x\ t) VmAC\ 



X 



exp 



i=l 



$^(A'^^A" + rV + CV^ 



a=2 



det{N' 



(5.179) 
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Since N' does not contain x or x, we can immediately perform the integrations 
over these fields, yielding 1. 

To make the indices in the matrix A^' look nicer, we perform the change of variables 
from A to A' defined by (A')" = A°^^ for 1 < a < i (and recalling that A^''^^ is a 
shorthand for A^); the variables A are left as is. The Jacobian is (— l)'^^^^^^ So, 
dropping primes, we have 



-'ni,...,n^ 



,fc(s) 



X exp 



5^(rTA"-^ + CV) 



a=2 



det{N" 



(5.180) 



where 



N" 







6. 



if a = /3 - 1 
if a < (3 - I 
A°TA/3 if a>/3 



(5.181) 



Finally, in order to put our scalar products of complex fermions in the standard 
forms C,'^C, and A"'"A, we anticommute all the bilinears (obtaining a minus sign); and in 
order to keep the indices in a natural notation, we replace A^" by its transpose. After 
a renaming A — )■ and ( —> rj, the result is fl5.158p / fl5.159p . where (A^")"^ = + M. 

This proves part (i) of Theorem 12.91 when the X^"'^ are real matrices of rank 
min(nQ,, Ua+i) lying in a sufficiently small neighborhood of InaUa+^i and s is a positive 
integer; the general validity of the identity then follows from Proposition 12.181 □ 



We now turn to the evaluation of the Grassmann integral fl5.158p / (15.1590 for 
bni,...,nt;k{.s) . The integrand depends on the fields only through the scalar products 
{'4'"'^i^^Ya^l3=i and {t]°'^ 'n°'Ya=2- "^^^ therefore apply Proposition l5.6l once to the en- 
tire set of variables ip, ip and separately to the variables r]", f/" for each a {2 < a < i). 
We therefore introduce indeterminates X = (2:^/3)0^^=1 and y = {1/0)0=2^ along with 
the corresponding differential operators d/dxa/s and d/dya- Using Proposition 15.61 in 
the form fl5.149bp . we obtain (after renaming y2, . . . ,ye as yi, . . . , ye-i) 



Jn\,...,ni 



'1-1 

a=l 



d 



dx, 



a+l,a 



_d_ 

dye 



£-1 



det(X)'=n^" 



a=l 



X=y=0 

(5.182) 
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where now M reads 

/ 1 + 



M 



8 


8 


8 


8 


dxii 


8x12 


8xi3 


8xii 


8 


1 + 

8X22 


8 


8 


8yi 


8X23 


8x21 





d_ 


1 + ^ 


8 


8y2 


8X31 



v 



8yi 



1 + ^ 



(5.183) 



We can now apply the translation formula (15. ip to the exponential of the differential 
operator: this transforms ntl^l/cT" i^^o Iltl^ll + Va)"^" and Xa/s into x'^f 



5, 



a,l3+l 



. Thus 



e-1 



hn,,...,nAs) = (-l)''(detM)-^ det(X')'=n(l + 



a=l 



XaH + 



(5.184) 



X=y=0 



Finally, we observe that M does not contain the differential operators d/dxap with 
a > /3, so in X' we can set those variables Xap to zero immediately. We thus have 
( I5.184P where now 

/ Xii Xi2 Xii ■■■ xu \ 

1 3^22 ■ ■ ■ X2I 

1 X33 ■ ■ ■ X3£ 



X' 



(5.185) 



\ ■■■ 1 Xu J 

Since the evaluation of (I5.184p will involve a recursion in £, it is convenient to in- 
troduce an infinite set of indeterminates {xQ,^}i<a</3<oo, along with the corresponding 
set of differential operators dap = d/dxap-, as well as another infinite set of indetermi- 
nates {aa}i<a<n~i (ouly finitely many of these will play any role at any given stage). 
Then define, for each i > 1, the quantities 



De{a) 



and 





/ l + dn 


di2 




die 








l + d22 


^23 


d2e 


det 







a2 1 


+ ^33 


d3e 




\ 







a 


e-i l + d( 








^ Xii X12 


Xl3 ■ ■ 


■ xie \ 










1 X22 


X23 ■ ■ 


■ X2e 






= det 


1 


X33 ■ ■ 


■ X3e 










I •■• 


1 


xee J 





\ 



(5.186) 



(5.187) 
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where we also set -Do(o) = 1 and Xq = 1. Note also that -D^(a) [resp. Xi] involves 
only those da/3 [resp. Xap] with a < (3 < i. 

Given a formal indeterminate s and a nonnegative integer k, our goal in the 
remainder of this subsection is to compute the expression 



Deia)-'X^\^^^ , (5.188) 



lx=o 



which abstracts the relevant features of (15.1841) [ clS 5 IS cl formal variable, the choice 
between s and —s has no special role and is made here just for convenience]. Since 
Di{a) is a polynomial in the quantities {da/B} and {tta} with constant term 1, Di{a)~'^ 
is here to be understood as the series 

Dei<^)-' = E / (5-189) 
h=o ^ ^ 

which, when applied to X^ as in (I5.188p . can be truncated to h < ki. We shall prove 
the following: 

Theorem 5.7 With the definitions above, we have 

D,(a)-'X}U__„ = *! (7) n E*! (7?,") ■ (5.190) 



a=l 6=0 

We remark that 

k 

k-h I h\ \ k 



b=0 

although we will not use this. 

What we shall actually need is a specific corollary of Theorem 15. 71 Let us introduce 
the variables {ya}i<a<oo and the associated derivatives da = d/dya, as well as the 
further indeterminates (or nonnegative integers) {mQ,}i<Q,<oo- 

Corollary 5.8 With the definitions above, we have 



De{-d)-'X^l[{l + ya] 



e-i k-i 

= (-i)''n n(^+^"+^) (5-192) 

X=y=0 a=0 i=0 



with the convention mo = 0. 

Given Corollary 15.81 the proof of Theorem 12.91 part (ii) is a triviality: 

Proof of Theorem I2.9L part (11), given Corollary 15 .Si It suffices to rec- 
ognize det[M] and det(X') as the operators D(^{—d) and Xi in Corollary 15.81 The 
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sign (—1)'^^ in (15.1841) combines with the one in (15.1921) . leaving exactly the pref actor 
claimed in Theorem 12.91 □ 

Let us next show how to deduce Corollary 15.81 from Theorem 15. 7t 



Proof of Corollary 15.81 given Theorem 15.71 We evaluate the left-hand side 
of (I5.192P by using (I5.190p with replaced by —da- Since 

'TTLr 



and everything is factorized over a, we obtain 



yc=0 V 



(5.193) 



a=l 



Then for each a we have 

k 



b=0 



-s-b 

k-b 



X=y=0 



TTLr 



— S 
k 



-1 k 



HE*! 

a=l fe=0 



-S-b 
k-b 



rrir 



b 

(5.194) 



k\Y,i-^)'-'( 

6=0 ^ 

{-irkif^f 

h—n ^ 



s + k-1 
k-b 



(-1) 



b ( 

b 



s + k-l 
k-b 



-l)''k\ 



6=0 

s + k + nia — 1 
k 

k-l 
i=0 



ma 

b 



(5.195) 



where we used the Chu-Vandermonde convolution (Lemma IB. 21) in going from the 
second to the third line. The prefactor k\ ("^'') gives an analogous contribution with 
mo = 0. □ 

Finally, we turn to the proof of Theorem 15. 7[ We shall need two main lemmas: one 
that essentially provides an inductive step, and another dealing with a special sum of 
multinomial coefficients. Henceforth we shall write Di as a shorthand for D^{a). 

We start with a pair of easy recursive formulae, obtained by expansion of the 
determinant on the last column: 



Lemma 5.9 



Q = l 



(the empty product aa - ■ ■ a^-i for a = I should be understood as 1) and 



Xe = ^(-1)^-"X, 



a—lXai ■ 



(5.196) 



(5.197) 



a=l 



70 



Related formulae for the determinant of a Hessenberg matrix can be found in |10H 
[102]. 

The induction lemma is the following: 

Lemma 5.10 Let m > 1 and ci, . . . ,Cm > be integers, and let t he an indetermi- 
nate. Then 



m m 
a=l a=l 



Xc,m=0 bl,...,bm>0 ' a=\ 



61 H Yhm= Cm 

m—1 m—l 

X 

a=l a=l 



Proof. Notice, first of all, that the factors and for a < m do not play any 
role, i.e. we can rewrite the left-hand side of f l5.198p as 



'm—l \ /m—l 



n (^^"^^*^^'"L..=o) n (5-199) 

,a=l / \a=l / 

and concentrate on the central factor alone. To compute I we expand 

D^~^^ and using Lemma 15.91 

^ \6iH Vhm)\ bi,...,bm J """^ 



6i,...,6™>0 

m 

X 

a=l 



II [(-1)— "(a„ ■ ■ ■ a^_,)D^_,d^^] ""^ (5.200) 

(s m 



bi,...,6™ > 

bl H +bm = Cm 



Since c?am^aml'^"™=n = ^b.b'h\, the two sets of summation indices, when combined 



inside D'S _ , must coincide, and we get 

6i,...,6m>0 ^ ^ 



61 H h bm = c„ 

m m 

X 



J] [(a, ■ ■ ■ a^_i)D„_i]'" {X^-if" . (5.202) 
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In the two final products, we can drop the factors D^q and Xq^ since Dq = Xq = 1. 
Reintroducing the missing factors from fl5.199p . we obtain f l5.198p . □ 

We will now apply Lemma [5 .101 for i "rounds", starting with the initial conditions 
m = i, t = —s — k, Ci = . . . = Q_i = and q = k. At round i we have m = i + l — i. 
Let us denote hj c\, . . . , c^+i_j the parameters {cj} immediately before entering round 
i, and let us denote by 6^, ... , the summation indices in round i. We therefore 

have the initial conditions 



0, 



k 



and the recursion 



for 1 < j < i — i 



[cf. f l5.198p ]. The summation indices b\, . . . , obey the constraint 



e+i-i 



Using f l5.203p and (15. 2 04 p we prove by induction that 



C 



e+i-i 



'+2-i 



for 2 < i < 



h=l 



It is convenient to arrange the summation indices {6*} 



i+j<e+i 



into a matrix 



B 



ft2 
"2 



"2 



(5.203) 
(5.204) 

(5.205) 



(5.206) 



(5.207) 



in which the first row sums to k and, for 2 < i < i, row i and column i + 2 — i have 
the same sum. Such matrices can be characterized as follows: 



Lemma 5.11 Fix integers i > 1 and k > 0. For a matrix B = (6*)j+j<^+i of 
nonnegative integers as in ( 5.207] ), the following conditions are equivalent: 

(a) The first row sums to k and, for 2 < i < i, row i and column i + 2 — i have the 
same sum. 



(b ) For 1 < h < i we have 



b] 



k. 



l<i<h 
1 < i < ^ + 1 - 
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(c) There exist nonnegative integers bj, b^_i, • • • , &2 completing the matrix B to one 
in which all the row and column sums are equal to k. (Such numbers are obvi- 
ously unique if they exist, and must lie in the interval [0,fc].j 

Note in particular that statement (b) with h = £ tells us that the first column sums 
to k. 

Proof, (a) =^ (b): By induction on h. By hypothesis the equality holds for h = 1. 
Then for h > 2 we have 

E E ^}=E^?-E b\^2-, = (5.208) 

l<i<h l<i<h-l ^<j<^-h l<i<h-l 

by hypothesis. 

(b) =^ (c): It is easily checked that the definition 

= ioT2<h<i (5.209) 

1 < j < h - 1 

i<i<^ + i- 'i 

does what is required. 

(c) =^ (a) is obvious. □ 

Remark. An analogous equivalence holds, with the same proof {mutatis mutan- 
dis), for matrices B of nonnegative real numbers where is a fixed nonnegative real 
number. □ 

The ^-fold application of Lemma [5. 101 with the initial conditions f l5.203p gives rise 
to a sum over matrices B satisfying the equivalent conditions (a)-(c) of Lemma [5. Ill 
To each such matrix there corresponds an £-tuple c = (c^, . . . , c^) of integers in the 
range [0,k], where c* = is the sum of row i (it also equals c^_|_i_j) and of 

course = k. It will be convenient to partition the sum over matrices B according 
to the vector c, so let us denote by B{c) be the set of matrices satisfying the given 
conditions with the row sums c: 

i+l-i 

B{c) = [b = (6;).+,<,+i G N^e^+i)/^: ^ 6} = for 1 < z < £ 

i=i 

j-i 

and = ^ for 2 < i < . (5.210) 

h=l 

The summand is then 

(n(':')"(-„..*,_,))(n-'*'') 
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Furthermore, it follows from Lemma 15.11( b) that 

i—a a 

J] J] 6} = k - c^+^-" . (5.212) 

Let us now show how to perform the sum over matrices B with a given vector c: 
Lemma 5.12 Let £ > 1, /c > and c^, . . . , > he integers, with = k. Then 

E ri(,: ) = ri(*)' 

Before proving Lemma I5.12[ let us show how it can be used to complete the proof 
of Theorem 15.71 

Proof of Theorem 15.71 given Lemma 15.121 Summing f l5.21ip over B e B{c) 
and using (I5.212p and fl5.213p along with t = ~s — k and = k, we obtain 

n(-r 1 - ( ■)) (n--'-^) ^ « (-;)n(-;-t) 

(5.214) 

where '5^ = k — c^~^^~". The sum over (f,...,c^ — or equivalently over c^, . . . , c^~^ — 
now factorizes and gives precisely f l5.190p . □ 



Proof of Lemma 15.121 The sum on the left-hand side of f l5.213p is nontrivial 
because the "row" and "column" constraints (15.2101) are entangled. We shall replace 
one of the two sets of constraints (say, the "column" one) by a generating function: 
that is, we introduce indeterminates C,j (1 < J < ^) and consider 



wo = E m-' 'XIL. \. ■ (5-215) 

fhH 7=1 1=1 ■ ■ • ' ^+1- 



The sum now factorizes over rows: we have 

-l-i 



-l-i ^ " 



n E ^0 ^^-^^^^^ 
i=i \ j=i / 

(6 + • . . + ^ifi^i + ... + ^i-ifi^i + ■■■+ ii-2f ■ ■ ■ it . (5.216c) 
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We must now extract the coefficient of the monomial 



e.^ur ■■■if (5.217) 

in <l>£_fc c(0- We ffist extract [^f] from fl5.216cp : here ^£ occurs only in the ffist factor, 
so we get (ij) (^1 + • • • + ii^iY"'^ times the remaining factors, i.e. 

Qt) (^1 + • • • + + ■■■ + ii-2f ■ ■ ■ it . (5.218) 

We can then extract [if_-^ in the same way, and so forth until the end, yielding 
right-hand side of (15.2131) . □ 

Let us conclude by apologizing for the combinatorial complexity involved in the 
proof of Theorem 15. 7[ The simplicity of the final formula (15.1901) . together with the 
simplicity of (15.2131) and the miraculous simplifications observed in its proof, suggest 
to us that there ought to exist a much simpler and shorter proof of Theorem 15. 7[ But 
we have thus far been unable to find one. 



i=2 



6 Proofs of diagonal-parametrized Cayley identi- 
ties 

In this section we prove the diagonal-parametrized Cayley identities (Theorems l2.10l 
and l2.11l) . We give two proofs of each result: the ffist deduces the diagonal-parametrized 
Cayley identity from the corresponding Cayley identity by a change of variables; the 
second is a direct proof using a Grassmann representation of the differential operator. 



6.1 Diagonal-parametrized ordinary Cayley identity 



Proof of Theorem I2.1U[ We change variables from (xjj)"^^^ to new variables 

(ti)"=i and {yij)i<i^j<n defined by 

U = Xii (6.1a) 

We also set ya = 1 for all i and define the matrices Tq, = diag(t"'), Ti_ct = diag(tj^^"') 
and Y = (yij), so that 

X = T^YT,_^ (6.2) 

and hence 

(detX)^ = (^n^lj (detr)^ (6.3) 
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A straightforward computation shows that the differential operators (vector fields) 
d/dxij can be rewritten in the new variables as 



d 



dx 



f~ 'f 



Wi 



if i = j 



if 2 ^ j 



(6.4) 



Let us denote by the matrix of differential operators whose elements are given 
by the right-hand side of (16 .4^ : please note that each element commutes with each 
other element not in the same row or column [that is, (Aq,)^^ commutes with 
whenever i ^ i' and j 7^ j']. 

Now (I6.3p . considered as a function of t = (tj), is of the form const x Hj^ij 
therefore, acting on (16.31) . each operator Ud/dti is equivalent to multiplication by s. 
It follows that the action of A^ on (16. 3p is identical to that of A^ defined by 



SjlJ 



s-aiJ2 Vik -pf— - (1 - tti) E 7¥- 



if i = j 



(6.5) 



V--ai7--(l-"j) d 



both commute with (A 



a,s)i'j' whenever i ^ i' and 



Furthermore, {^a)ij and {^a. s)ij 
3 ^ f- 

Now suppose we have a product (Aci)j^j^ ■ ■ ■ (Aq)^^^^ acting on (16. 3p . in which 
are all distinct and also are all distinct (this is the case that 

will arise when we take the determinant of any submatrix). Then the rightmost 
factor {Aa)iijf can be replaced by (Aa,s)j^j^; and the same can be done for the 
other factors by commuting them to the right, changing them from Aq, to Act,s, 
and commuting them back. It follows that, under the given condition on the indices, 
{A^)i,,, ■ ■ ■ {A^)i,,,{detXy = {A^,shn " " " (A„,s)i«,(det X)^ 

Translating back to the original variables (a;ij)"j=i, Aq, equals Da^s defined by 



d 

Oxij 



s-aiJ2 Xik iJ- - (1 - tti) E xii -fir. 
k^i /^i '^■^^■^ 



iii= j 
if i ^ J 



(6.6) 



and the Cayley formula (12. 2p tells us that 
det((L'„,,)/j) (detX)" = s(s + 1) ■■• (s + A; - 1) (det X)^-^ e(/, J) (det Xjcjc) . (6.7) 



On the other hand, -Da,/?,; 
diag(Xj^~^'), so that 



XfiDa,sXi-f3 where X^ = diag(x§) and Xi_^ 



det((D 



a,l3,s)lj) 



n 



n 



det{{D 



cx,s)lj) 



(6.8) 
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Combining (|621) and ([61]), we obtain (EJ?]). □ 



Now let us show an ab initio proof based on Grassmann representation of the 
differential operator det{{Da,s)ij)- Since Da^s contains not only terms d/dxij but 
also terms Xijd/dxij, the role played in Section [5] by the translation formula (15. ip 
will here be played by the dilation-translation formula flB.60|) : see Section [B.4I for 
discussion. 



Alternate proof of Theorem I2.10L For notational simplicity let us assume 
that the diagonal elements Xu are all equal to 1; the general case can be recovered by 
a simple scaling. 

Consider the matrix of differential operators D^^s defined by (16. 6p . In terms of 
the matrices E^^ defined by (13. 4p . we can write 



Dc.,s = si + J2 (^'' - x^j[a,E'' + (1 - aj)E^^] 



d 



dxij 



(6.9) 



As before, we introduce Grassmann variables i]i,rii {1 < i < n) and use the represen- 
tation 



det((D„,,),j) = e(/,J) J Vn{v,r])(^l[m)^^ 



X 



exp 



s^ViVi + "^[ViVj - Xij{air]ir]i + (1 - aj)r]jr]j)] 

i i^j 



d 



dx 



(6.10) 



Let us now apply this operator to a generic polynomial f{X): using the dilation- 
translation formula flB.60P on all the variables Xij {i ^ j), we obtain 



/(^|(1 - air]ir]i)xij{l - (1 - aj)VjVj) + ^i^JJ.^.j ■ (^-H) 

Note that the diagonal terms remain unchanged at their original value Xu = 1. Defin- 
ing the diagonal matrices M„ = diag(l — aif]i7]i) and Mi_q = diag[l — (1 — ai)f]ir]i], 
we see that the argument of / is 



X' = M^XM,_^ + f^ri" 



(6.12) 



(note that also the diagonal elements i = j come out right). We are interested in 
f{X) = det(X)'*, and we have 



detX' = (detM„)(detMi_„)det(X + M-^r7r7^Mf^ 
= (detM„)(detMi_c,)det(X + r/r/^ 



(6.13a) 
(6.13b) 
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where M^^rj = rj and rj'^M^]:^ = rp- by nilpotency. The factor (det Ma)*(det Mi_c»)'' 
exactly cancels the factor exp[s r^jr^j] in the integrand, and we are left with 

det((D«,,),j)det(X)^ = e(/,J) j P„(77, 77) ( J] r^T^)^^ det(X + r/r/T)^ (6.14) 

which coincides with f l5.13ap . Arguing exactly as in f l5.13p -f l5.15p . we therefore obtain 
det((D«,s)/j) (detX)" = s(s+l) ■ ■ ■ (s+fc-1) (det J) (det X/cjc) . (6.15) 

□ 



6.2 Diagonal-parametrized symmetric Cayley identity 

An analogous proof gives the symmetric analogue: 

Proof of Theorem 12 .111 We change variables from {xij)i<i<j<n to new variables 

(ti)"=i and {yij)i<i<j<n defined by 



ti 



Xi 



ii jj 



(6.16a) 
(6.16b) 



We also set yu = 1 for all i, use the synonymous yji = yij for j > i, and define the 
matrices T = diag(tji^) and Y = (yij), so that 



X 



TYT 



and hence 



(detX)^ = {UtAidetYY 



(6.17) 



(6.18) 



,i=l 



A straightforward computation shows that the differential operators (vector fields) 
d/dxij can be rewritten in the new variables as 



d 



t7 



dx 



d 



d 



+ 2+ 2 d 



(6.19) 



Let us denote by A the matrix of differential operators whose elements are given by 
the right-hand side of (16.191) : please note that each element commutes with each other 
element not in the same row or column [that is, Aij commutes with Aj/j/ whenever 
i ^ i' and j 7^ j']. 
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Now flG.lSp . considered as a function of t = (ti), is of the form const x Yli^i 'j 
therefore, acting on (16.181) . each operator tid/dti is equivalent to multiphcation by s. 
It follows that the action of A is identical to that of defined by 



(A 



SjlJ 



t7' 



ik 



d 



ik 



d 



(6.20) 



Furthermore, Ajj and {As)ij both commute with {As)i'j' whenever i ^ i' and j ^ j'. 

Now suppose we have a product Aj^jj ■ ■ ■ Aj^^^ acting on (16.181) . in which ii, . . . ,ie 
are all distinct and also ji, . . . ,je are all distinct. Then the rightmost factor Aj^j^ 
can be replaced by (As)j^j^; and the same can be done for the other factors by com- 
muting them to the right, changing them from A to A^, and commuting them back. 
It follows that, under the given condition on the indices, Aj^^— ■ ■ Aj^j^(det X)'^ = 
(A.)im---(As).«,(detX)^ 

Translating back to the original variables {xij)i<i<j<n, A^ equals D^^™ defined by 



sym\ 



d 



2 S ^ 

k^i 



ik 



d 



'5x 



ik 



2 "^li 



l^i 



d 



if i ^ j 



(6.21) 



and the Cayley formula (12.21) tells us that 



det((L' 



sym^ 



IJ. 



(detX 



sym\s 



+ 



s + 



k-l 



(detX^n'~'e(J, J) {AetXfJ^.) 

(6.22) 

On the other hand, -Dj™ = Xj^D^'^Xi-fj where = diag(a;§) and Xi_/3 = 
diag(x|~^'), so that 



det((D^y'" 



f3,s )IJ) 



1-/3, 



~i€l ^ ^jeJ 

Combining (l6:22|) and dOS]) . we obtain ([232]). □ 



(6.23) 



Now let us show the proof based on Grassmann representation of the differential 
operator: 

Alternate proof of Theorem 12.111 For notational simplicity let us assume 
that the diagonal elements Xu are all equal to 1; the general case can be recovered by 
a simple scaling. We are therefore using the matrix X^^^™ defined by 



sym\ 



1 for i = j 
for i < j 
for i > j 



X 



(6.24) 
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Consider the matrix of differential operators D^J^ defined by fl2.30p . In terms of 
the matrices E''^ defined by (13.41) . we can write 



i<j 



d 



dxij 



(6.25) 



We introduce Grassmann variables r]i,r]i {1 < i < n) and use the representation 
det((D|;n/.) = e(/,J) f 'Dn{v,v)[llm)^^^^ x 



exp 



d 



i<j 



dxi 



(6.26) 



Applying the dilation-translation formula (1B.60P on all the variables Xij {i < j), we 
obtain 



det((/^^^^)/j) /({xij}i<j) = e(/,J) J Vn{ri,ri) (^Y[rir]^ ^^exp s^rji^ 



X 



i<j 



(6.27) 



Note that the diagonal terms remain unchanged at their original value Xu = 1. Defin- 
ing the diagonal matrix M = diag(l — ^fjirji), we see that the argument of / is 

X«ym/ _ MX^^"'M + \{m^ -m^) (6.28) 

(note that also the diagonal elements i = j come out right). We are interested in 
j(^^sym) ^ det(X"y"^)^ and we have 



detX^y'"' = {del det[X + lM-\f]r]^ -r]f]^)M- 
= {det Mf det[X + - vf)] 



(6.29a) 
(6.29b) 



where the last equality again follows by nilpotency. The factor (detM)^'' exactly 
cancels the factor exp[s Xli^i^i] integrand, and we are left with 



det((I^^,, )ij 



^^■"^--^ det(X^n^ = e(/, J) J P„(r7,r^) (jl^^),.^. det[X^^^ mm^-m^)y 

(6.30) 



which coincides with f l5.4ip . The remainder of the proof is as in f l5.4ip -( l5.48p . □ 
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7 Proofs of Laplacian-parametrized Cayley identi- 
ties 



In this section we prove Theorems 12.121 and I2.14[ The proofs use a Grassmann 
representation of the differential operator, and are closely patterned after the proofs 
of the ordinary and symmetric Cayley identities in Sections 15.11 and 15. 3^ respectively. 

7.1 Laplacian-parametrized ordinary Cayley identity 

Let us begin by recalling the definitions of the matrices arising in Theorem 12.121 

Xij if « 7^ j 



if z = J 



and of course T = diag(tj). In order to maximize the correspondences with the proof 
in Section 15.11 it is convenient to prove instead the dual (and of course equivalent) 
result for co/wmn-Laplacian matrices: 

Xij if « 7^ J 



if ^ = j 



(note that (^^oi-Lap _ ^row-Lap jj-^ what follows we shall drop the superscripts "col- 
Lap" in order to lighten the notation, but it is important to remember the definitions 

Proof of Theorem I2.12[ We introduce Grassmann variables rji^fji {I < i < n) 
and use the representation 

det([/ + d)- det(9) = j Vniv, v) (e''''^'' - l)e^"''^'' , (7.3) 

It is convenient to introduce the Grassmann quantities 9 = X^^t]^ and 6 = ^^fje, 
so that f]'^Urj = Q^Q. Since 6^ = 6^ = 0, the first exponential in (17.31) is easily 
expanded, and we have 

det([/ + d)- detid) = J V^iv, v) (0^©) e''^^'' • (7-4) 
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Let us now apply (17.41) to det(T + X)'* where s is a positive integer, using the trans- 
lation formula (15. ip : by (17.2b ) we get Xij — )• Xij + fjirij for i ^ and by (17.2b ) this 
induces X — )■ X + fjrf- — dia.g{Qrji) . We therefore obtain 

[det(f/ + d)- det(9)] det(T + X)" 

P„(r/, n) (©""e) det(T + X + f^rj^ - dmg{Qr^i)y . (7.5) 



But now comes an amazing simplification: because of the prefactor 6"'"B and the 
nilpotency 0^ = 0, all terms in the expansion of the determinant arising from the 
term diag(9?7j) simply vanish, so we can drop diag(6?7j): 

[det(f/ + (9) - det(9)] det(T + X)' = j Vniv, f]) (6^6) det(T + X + r/r/^) . (7.6) 

Assuming that T + X is an invertible real or complex matrix, we can write this as 
[det(;7 + d)- det{d)] det{T + X)^ (7.7) 
= det(T + X)^ J Vn{ri, r/) {fUr]) det[I + (T + Xy'fjrj^Y ■ 

Let us now change variables from {ri,f]) to ij]' ,fj') = {rj, (T + X)^^ff): we pick up a 
Jacobian det(T + X))~^, and dropping primes we have 

[det(t/ + d)- det(a)] det(T + X)" (7.8) 
= det(T + Xy~^ j Vniv, V) {f{T + X^)U7]) det(J + m^Y ■ 

But since X is column-Laplacian, we have UX = and hence X'^U = 0, so that 
f]^iT + X'^)Ur] reduces to f]^TUr]lL_\ This expresses the left-hand side of the identity 
as the desired quantity det(T + X)'^"^ times a factor 

P{s, n,T) = J Vniv, n) ifTUr]) det(/ + m'^Y , (7.9) 

which we now proceed to calculate. The matrix / -|- fjr]^ is a rank-1 perturbation of 
the identity matrix; by Lemma IB. Ill we have 

det{I + f]r]^y = [l-fr])-' (7.10a) 

= ( /) (^V . (r.iob) 



For a row-Laplacian matrix we would have instead chosen to muhiply 7777"^ hy {T + X) ^ on 
the right, and then made the change of variables from (77, fj) to (77', ff) = {{T + X"^)^^ri, fj), picking 
up a Jacobian (det(r + X^))"^ = (det(r + X)y^ and obtaining a prefactor TfU{T + X^)ri; then 
UX'^ = because X is row-Laplacian. 
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Now 

f]^TUri = ^tif]irii + ^^tif]irij , (7.11) 

i i jjti 

but the terms fjirjj with i ^ j cannot contribute to an integral in which the rest of 
the integrand depends on f] and r] only through products fjirji] so in the integrand we 
can replace fj^TUr] with ff-Trj. Since 

Jv^{7],r])m{fvY = {n-l)\6e,n-i (7.12) 
for each index i, it follows that 

P{s,n,T) = {J2t^){-lr-'(~^^{n-lV. (7.13a) 

i ^ ^ 

= (^^ti)s(s + l)---(s + n-2) . (7.13b) 

i 

This proves fl2.36p when X is a real or complex matrix, t = (tj) are real or complex 
values such that T + X is invertible, and s is a positive integer; the general validity 
of the identity then follows from Proposition 12.181 □ 



7.2 Laplacian-parametrized symmetric Cayley identity 



Now we prove the corresponding result for symmetric Laplacian matrices. Let us 
recall the definitions: 



sym-Lap \ 



X 



X 



iii<j 
if i > j 
iii = j 



(7.14a) 



iff 



iSym-Lap\ 



d/dXij if i < j 

iij - - \ d/dxji ifi>j 
ifi=j 

Once again we drop the superscripts "sym-Lap" to lighten the notation. 



(7.14b) 



Proof of Theorem I2.14[ Again we introduce Grassmann variables rjijfji (1 < 
i <n) and apply the representation (17. 4p to det(T-l-X)* where s is a positive integer. 
Using the translation formula (15.1 p and defining and G as before, we obtain 

[det(f/ + d)- det(9)] det(T + Xf 

Vr^iv, v) (ee) det(T + X + r]ri^- 7]f - dmg{f]iQ) - diag(er7,))' • 

(7.15) 
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We again argue that because of the prefactor 6"'^6 and the nilpotencies 6^ = 6^ = 0, 
all terms in the expansion of the determinant arising from the terms diag(f/jO) and 
diag(6r7j) simply vanish, so we can drop these two terms: 

[det{U + d)- det{d)] det(r + X)' 

Vniv. n) (00) det(T + X + r/r/T - 7]fy . (7.16) 

Assuming that T + X is an invertible real or complex matrix, we can write this as 
[det(f/ + (9) - det(9)] det(T + Xf 

= det(T + Xy j Vn{r], f]) {fUr]) det[I + (T + X)-\r]r]^ - r]f)Y .(7.17) 

Let us now change variables from (r/, r/) to {ri',f]') = {rj, (T + X)''^f])] we pick up a 
Jacobian (det(T + X))~^, and dropping primes we have 

[det{U + d)- det{d)] det(T + xy 

= det(T + Xy~' J Vniv, n) {f{T + X^)U7]) det[/ + r/r/^ - (T + X)-'rif{T + XfY ■ 

(7.18) 

But since X is symmetric Laplacian (hence column-Laplacian), we have X'^U = 
0, and the prefactor reduces to ff-TUrj. Now we apply the Corollary IB. 121 to the 
determinant expression, and obtain 

[det(f/ + d)- dei{d)] det(r + Xy 

= det{T + Xy-' jVr,{v,v) {fTUiiyi-fii)-^' . (7.19) 

This expresses the left-hand side of the identity as the desired quantity det(T + X)'^~^ 
times a factor 

P(s, n, T) = j V^{r,, f]) {fTUr,) (1 - fr^)'^' , (7.20) 

which we now proceed to calculate. As in the row-Laplacian case, the terms in 
ff-{TU)ri of the form f/jZ/j with i ^ j cannot contribute to the integral, as the rest 
of the integrand depends on r] and fj only through products f/jr/j, so we can replace 
fj^{TU)rj with fj^Trj. Since 

Vn{v,f])m{fr]Y = {n-l)\5e,n-i (7.21) 



for each i, it follows that 



P(.,n,T) = {Y.t^){-lr-'(^\){n-lV. (7.22a) 
= (^^U^2s{2s + l)---{2s + n-2) . (7.22b) 
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This proves (12.401) when X is a symmetric real or complex matrix, t = (tj) are real or 
complex values such that T + X is invertible, and s is a positive integer; the general 
validity of the identity then follows from Proposition 12.181 □ 



8 Proofs of product-parametrized and border- 
parametrized rectangular Cayley identities 

In this section we prove the product-parametrized and border-parametrized rect- 
angular Cayley identities (Theorems 12.161 and 12.171) and then discuss the close rela- 
tionship between them. 



8.1 Product-parametrized rectangular Cayley identity 

Before beginning the proof of Theorem 12.161 let us observe that the quantity 
det{MA) appearing in the statement of the theorem [cf. (I2.44p for the definition of 
the matrix M in terms of X and B] has an alternate expression as follows: 

Lemma 8.1 Let A, B, X and M be as in Theorem \2. 16[ Then 

det(MA) = J2 ^{I,L){det{B^A)jL){det{XA)icLc) (8.1) 

L C [m] 
\L\ = k 

Proof. The definition (I2.44p of M can be rewritten as 

Mic^ = Xj.^ (8.2a) 

Mi, = (5^)j, (8.2b) 

Therefore 

{MA)ic, = {XA)i., (8.3a) 

{MA)j, = {B'^A)j, (8.3b) 

for any matrix A. We now apply multi-row Laplace expansion flA.8p with row set / 
and summation variable L; this yields (18. ip . □ 

We are now ready to prove Theorem I2.16[ In order to bring out the ideas behind 
the proof as clearly as possible, we will first fully develop the reasoning proving the 
"basic" identity (I2.42p — which is actually quite simple — and then describe the 
modifications needed to handle the all-minors case. 



85 



Proof of Theorem 12.161 Let us apply Corollary lB.6l to the determinant det{dB) 
and then introduce a Grassmann representation for the resulting block determinant: 
we obtain 



det{dB) = det 





d 


-B 


In 



Vmi^, tp) Vniv, V) exp [tp^dl] - Tj^Bi) + if 7]] 



(8.4) 

where ipii'ipi [1 < i < m) and rii,7]i [1 < i < n) are Grassmann variables. By the 
translation formula (15.11) . we have 

det{dB) f{X) = j VU^, ^) V^T], exp [f{-Bi, + r^)] f{X + ^r,^) (8.5) 

for an arbitrary polynomial /. We shall use this formula in the case f{X) = det(Xy4)^ 
where s is a positive integer, so that 

det{dB) det{X Ay = j P^(V^, iP) V^iv, v) exp [fi-BiP + r])] det {X + V'r/^)A ' . 

(8.6) 

It is convenient to introduce the shorthand 

^trans ^ X + t^T]^ . (8.7) 

Suppose now that XA is an invertible real or complex matrix. Then we have 

^trans^ = {X + ^ri^)A = (X A) [/„ + ( (XA) V^) (r/^ A)] . (8.8) 

Let us now change variables from {ip^ip^rj^ff) to , ip' , r]' , fj') = {ilj,{XA)~^ip,'r],f]), 
with Jacobian det(XA)~^ = (detXA)"^. In the new variables we have (dropping 
now the primes from the notation) 

^trans^ = {X A){Im + tpv'^ A) , (8.9) 

and the translated determinant is given by 

det(X*'^^^M) = det{XA)det{Im + tpV^A) , (8.10) 

so that 

det{dB) det(XA)" = det{XAy-^ jVmiiJ,^)Vn{v,f]) e-''''^^^'^^'^ det{I + ^r]'^Ay . 

(8.11) 

Applying Lemma [B. Ill to the rightmost determinant yields 

det{I + ipr]^A) = {1 -iP^A'^r])-\ (8.12) 
so that we are left with the Grassmann-integral expression 

det{dB) det(XA)^ = dei{X Ay-^ j Vm{ip, ^) V^iv, v) 6"^^^'^+^^'' (1 - ^^^t^)-^ _ 

(8.13) 
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We have therefore proven that det{dB) det{X Ay equals the desired quantity det(XA)' 
multiphed by a factor 



b{s,A,B) 



(8.14) 



that does not involve the variables X (but still involves the parameters A and B\ 
Now, in the expansion of 



fc=0 ^ ^ 



(8.15) 



only the term k = m survives the integration over the variables ip: so we can replace 
(1 - tfj^A^r])-' in the integrand of flHTTD by (-1)™ {~^) {ip^ A^r])"". Moreover, the 
same reasoning shows that we can replace [tp^ A^^t])"^ by ml explip'^A'^r]). We are 
therefore left with a combinatorial prefactor 



— s 
m 



s(s + l] 



(s + m — 1] 



multiplying the Grassmann integral 

■ 7- ^ 

V^{ip,ip)Vniv,v) exp ( _ 





A' 


-B 


In 



det 





A" 


-B 


In 



(8.16) 



(8.17) 



which equals det(y4^-B) by Corollary IB. 61 

This proves the "basic" identity fl2.42p whenever XA is an invertible real or com- 
plex matrix and s is a positive integer. Now, if A has rank < m, then both sides of 
f l2.42p are identically zero; while if A has rank m, then XA is invertible for a nonempty 
open set of matrices X. The general validity of the identity fl2.42p therefore follows 
from Proposition 12.181 

Now let us consider the modifications needed to prove the all-minors identity 
f l2.43p . For a while these modifications will run along the same lines as those used in 
the proof of the two-matrix rectangular Cayley identity (Section 15. 5p . Thus, a factor 

^{IiJ) \ n^V') gsts inserted into the Grassmann integral fl8.4p -( l8l6l) : after the 



So 



change of variables (and dropping of primes) it becomes e(J, J) y Y\i.-^Aip)ip 
we have, in place of equation (18. 130 . the modified expression 

dei[{dB)jj] dei{XAy = e{I,J) det{X Ay-' J VUi^,^)V4r],r]) (HiXA^)^'^ ^ 

fj^T] - ff^Bip] (1 - %lj^A^r])-' . 



X exp 



Once again we argue that the integration over variables ip allows to replace 



(8.1 



{l-^^A^7])-' 



k\ exp[^'^A^?7] 



L19) 
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since, in both cases, only the k-th term of the expansion survives. We therefore have 
det[{dB)ij] det{X Ay = e(/, J) det{X Ay-^s{s + 1) ■ ■ ■ (s + - 1) 

(8.20) 

Now we perform the integration over rj and fj using Wick's theorem for "complex" 
fermions in the "source" form [cf. flA.94p ]. yielding 

det[{dB)ij] det{XAy = e(/, J) det{X Ay-^s{s + 1) ■ ■ ■ {s + k - 1) 



X 



Vmi^,^) (^l[{M^X^)ij^^^^^ explij^A^Bi:] . (8.21) 



Next we perform the integration over ip and ip using Wick's theorem for "complex" 
fermions in the "correlation function" form (]A.97p . yielding'^ 



J2 L) (det Ij.l) (det(AT5)KcLc) {deiiA" X^)kic) 

\K\ = \L\=n~k 

= <K,r){det{A'B)K^j){dei{A'X^)Ki<^) 

\K\=n-k 

= J2 '^') {det{A^B)Kj) {det{A^X^)Kcj.) . (8.22) 

\K\=k 

We now use e(/s:^ J^) = e{J, K) and e(J, J)e{J,K) = e{I,K); it follows that e(/, J) 
times (K2^ equals (jHH]). □ 

Remark. The s = 1 special case of the all-minors identity fl2.43p has an easy 
elementary proof, which actually proves a stronger result. Note first that by the 
multilinearity of the determinant det(Xy4) in the variables {xij}, we have 

det(Myl) = ^idB),^j)j det(XA) . (8.23) 

Moreover, because permuting the indices ji, . . . ,jk amounts to permuting the rows of 
B^ and hence permuting a subset of the rows of M, we have 

det(MA) = sgn(a) (jlidBy^,^^^^ det(XA) (8.24) 



^° Here we have made in ()A.97P the substitutions A -s- A^B, B (/),7c^, C {A^X'^)i,i<= 
I ^ K, J ^ L. 
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for any permutation a E Sk- Summing this over a G Sk, we obtain 

k\ det(MA) = [det{dB)ij] det{XA) , (8.25) 

which is nothing other than the s = 1 case of f l2.43p . But it is amusing to note that 
f l8.24p holds for each a G Sk, not just when summed over a E Sk- 

8.2 Border-parametrized rectangular Cayley identity 



Proof of Theorem I2.17[ We introduce Grassmann variables {1 < i < m) 

and rii,fji {1 < i < n — m). We use the representation 



det(9) = J 'Dm{i',i')'Dn-m{v,v) 
By the translation formula f lS.ip . we have 



(8.26) 



det(9)/(X) = Vrni^,^)V 

n—m 

[f], 7]) exp 



f{X + 



(8.27) 



for an arbitrary polynomial /. We shall use this formula in the case f{X) = det(X)' 
where s is a positive integer, so that 



det{d) det{Xy 



n—m 

{T],r]) exp 



det 



X + 



0J\V 



It is convenient to introduce the shorthand 



f8.28) 



^^^1 trans ^ X + 



(8.29) 



for the argument of det. 

Let us now assume that X is an invertible real or complex matrix, and change 



variables from 



V 



to 



7] 



X 



V 



with Jacobian (detX) ^. Dropping primes 



from the new variables, we observe that the expression for the translated matrix can 
be written as 



(xy 



X 



I + 



so that 



det(X)* 



(det X) det 



I + 



(if) 



(8.30) 
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Applying Lemma IB.llI to the rightmost determinant yields 



det 



/ + 



(8.32) 



so that we are left with the Grassmann-integral expression 



det (5) det(X)^ 



X exp 



(8.33a) 



det(X) 



s-l 



X exp [n^BX^ip + f]^BA^7]] (1 - ^^V^)"" . (8.33b) 

As the integrand depends on the Grassmann variables only through combinations of 
the form ipiipj, fjifjj and fjiipj (i.e. there is no ipirjj), we can drop all the terms fjiipj, 
as these terms would certainly remain unpaired in the expansion. This removes all 
the dependence on X in the integrand, and proves that det (9) is a Bernstein-Sato 
operator for det(X). We are left with the determination of the prefactor b{s), which 
is given by 



bis) 



V^i^, ^) Vn-Uv, V) exp [n^BA^T]] (1 - ^j^^Y 



(8.34) 



Integration over rj and fj gives a factor det(-By4'^) = det(y4i?^), while the integration 
over ip and ip is identical to the one performed in the case of the ordinary Cayley 
identity [cf. f l5.7p -f l5TT0|) ] and gives s(s + l) ■ ■ ■ (s+m — 1). This proves f l2.46p whenever 
X is an invertible real or complex matrix, B is an arbitrary real or complex matrix, 
and s is a positive integer. Now, if A has rank < n — m, then both sides of fl2.46p 

are identically zero; while if A has rank n — m, then X = 



A 



is invertible for 



a nonempty open set of matrices X. The general validity of the identity therefore 
follows from Proposition 12.181 □ 



8.3 Relation between product-parametrized and border- 
parametrized identities 

Let us begin by recalling the product-parametrized Cayley identity f l2.42p and the 
border-parametrized Cayley identity fl2.46p . writing the matrices A and B occurring 
in them as A^^\ B'^'^^ in the former identity and A^^\ B^^^ in the latter: 

(^M : det(9EW) [det(XA(°))]^ = b{s) det(A(o)T5(o)) [det(XAW)]^-i ^g g^^ 
(M : detC--^) det[-^y = b{s) det(A«i?«^) detC-^Y ' (8.36) 
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where b{s) = s{s + 1) ■ ■ ■ (s + m — 1). Here X is an m x n matrix, while A^^^ and 
are n x m matrices, and A^^^ and B^^^ are (n — m) x n matrices. Note that A^^^ 
and B^^^ must have full rank m, otherwise (I8.35P is identically zero; likewise, A^^^ and 
B^^^ must have full rank n — m. 

We will construct the matrices A^^^ and A^-^^ out of a larger (ra x n) matrix A, as 
follows: Let A be an invertible n x n matrix, and define 

= y4j,jm] = first m columns of A (8.37a) 
A^^) = (A"^)[m]c^^ = last n - m rows of A'^ (8.37b) 

Likewise, let B be an invertible n x n matrix, and define B^^^ and -B*^^-* by the same 
procedure. We then have the following facts: 

Lemma 8.2 Let X and Y hemxn matrices, with m < n; let A and B be invertible 
nxn matrices; and define matrices A^^\ A'^^\ B^^\ B^^'> as above. Then: 

(a) det(XA(°)) = (detA) det 

(b) det(rs(°)) = (det 5) det 

^c; det(v4(°)T5(°)) = (detv4)(detE) det(A(i)5(i)T) 

Proof, (a) First we expand the left-hand side using the Cauchy-Binet identity (IA.3P : 
det(XA(°)) = ^ (detX,z.)(detAS) = (det X,z.) (det A^j^]) . (8.38) 

|L|=m \L\=m 

Next we expand the right-hand side using multi-row Laplace expansion (1A.8P with 
row set [m], followed by the Jacobi identity f lA.7p : 




= e(L)(detX,L)(detA^T„j.) 

\L\=m 

= (^^t ^-^) (^^t ^)"' (^^t ^L,[m]) . (8.39) 

\L\=m 

Comparing fl8.38p and fl8.39p proves (a); and (b) is of course identical, 
(c) First we expand the left-hand side using Cauchy-Binet: 

det(A(°)T5(°)) = ^(det(A(°)T).^)(deti?i°i) = (det A^,^) (det 5l,m) • 

\L\=m \L\=m 

(8.40) 
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Next we expand the right-hand side using Cauchy-Binet and then using the Jacobi 
identity twice: 

det(A«S«T) = Yl (detAS)(det(S«^)x.) 

\K\=n—m 

J2 idet{A-^)K,[mr){det{B~^)K,imY) 

\K\=n—m 
\K\=n—m 

= (det A)-^ (det B)-^ ^ (det (det . (8.41) 

|L|=m 

Comparing f l8.40p and fl8.4ip proves (c). □ 

Using Lemma 18. 2[ we see immediately the equivalence of (I8.35P and f l8.36p when- 
ever A^") and A*^^) are related by fl8.37p and likewise for and B^^\ 

On the other hand, given any n x m matrix A^^^ of rank m, it can be obviously 
be completed to yield a nonsingular matrix A (which is invertible at least when the 
matrix elements take values in a field). Likewise, an (n — m) x n matrix A^^^ of rank 
n — m can be completed to yield a nonsingular matrix A^^. So to each A^^^ there 
corresponds a nonempty set of matrices A^^\ and vice versa. 

Remark. It is not in general true that every pair arises from a matrix 

A. Consider, for instance, m = 1 and n = 2: an easy calculation shows that for 
arbitrary A we must have A'^^^A'^^^ = 0. 

9 Conjectures on minimality 

Let us recall that any pair Q{s, x, d/dx) and h{s) ^ satisfying 

Q{s,x,d/dx)P{xy = b{s)P{xy'^ (9.1) 

is called a Bernstein-Sato pair for the polynomial P{x). The minimal (with respect 
to factorization) monic polynomial b{s) for which there exists such a Q is called 
the Bernstein-Sato polynomial (or 6-function) of P. Our Cayley-type identities thus 
provide Bernstein-Sato pairs for certain polynomials arising from determinants. But 
are our polynomials b{s) minimal? 

For the ordinary Cayley identities (Theorems 12. lH2l9|) . it follows from the general 
theory of prehomogeneous vector spaces [52l[60j that the polynomials b{s) found here 
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are indeed minimal, i.e. that the correct 6-functions are 



b{s) 



s(s + l) 



sis + 



s(s + 2) 



[s + n -1) 

> + 2m-2) 
■ (s + m 



for an n X n matrix 
for an n X n symmetric matrix 
for a 2m x 2m antisymmetric matrix (pfaffian) 
-1) 



s{s + 1) ■ ■ ■ {s + m 



for a 2m x 2m antisymmetric matrix (determinant) 

l){s + n — m) ■ ■ ■ {s + n — 1) 

for a pair of m x n rectangular matrices 



for an m X n rectangular matrix (symmetric) 

s{s + 2) ■ ■ ■ {s + 2m - 2){s + 2n - 2m + 1) ■ ■ ■ {s + 2n ~ 1) 

for an m X n rectangular matrix (antisymmetric) 



£ ni-1 

n II i^ + ^a- rii+j] 

a=l j=0 



for i matrices of sizes Ua x Ua+i 



(9.2) 

Indeed, the polynomials P occurring in these identities all correspond to relative 
invariants of prehomogeneous vector spaces: the ordinary, symmetric and antisym- 
metric Cayley identities [f l2.ip . f l2.4p and fl2.7p ] correspond to cases (1), (2) and (3), 
respectively, in Kimura's pUl Appendix] table of the irreducible reduced prehomoge- 
neous vector spaces; the one-matrix rectangular symmetric and antisymmetric Cay- 
ley identities [f l2.18p and (12.210 ] correspond to cases (15) and (13) in the same ta- 
ble; while the two-matrix and multi-matrix rectangular Cayley identities [ (]2.16p and 
f l2.23p ] correspond to prehomogeneous vector spaces associated to equioriented quiv- 
ers of type A |100] . Whenever P is a relative invariant of a prehomogeneous vector 
space, the general theory [521 [60] allows the immediate identification of a suitable 
operator Q{d/dx) — namely, the dual of P itself — and provides a proof that the 
corresponding b{s) satisfies degb = degP and is indeed (up to a constant factor) the 
Bernstein-Sato polynomial of pF^ 

For the Laplacian-parametrized identities (Theorems 12 . 121 and 12 . 14p . we conjecture 
that the polynomials b{s) found here are also minimal, i.e. that the correct 6-functions 



See [521 Corollary 6.1.1 and Theorem 6.1.1] [501 Proposition 2.22] for the first two points, 
and [521 Theorem 6.3.2] for the third. We are grateful to Nero Budur for explaining to us the 
connection between our results and the theory of prehomogeneous vector spaces, and in particular 
for pointing out that this connection provides a general proof of minimality. 
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are 

{s{s + 1) ■ ■ ■ [s + n — 2) for a Laplacian-parametrized n x n matrix 
s{s + |) ■ ■ ■ (s + for a Laplacian-parametrized n x n symmetric matrix 

(9.3) 

Perhaps these identities can also be interpreted within the framework of prehomoge- 
neous vector spaces; or perhaps an alternate proof of minimality can be found. 

It is even conceivable that the following general fact about Bernstein-Sato poly- 
nomials is true: 

Conjecture 9.1 Let P{xi, . . . ,Xn) ^ be a homogeneous polynomial in n variables 
with coefficients in a field K of characteristic 0, and let b{s) be its Bernstein-Sato 
polynomial. Then degb > deg P. 

Simple examples show that we need not have deg b > deg P if P is not homogeneous: 
for instance, P{x) = 1 — x"^ has b{s) = s. Moreover, slightly more complicated 
examples show that one can have deg 6 > degP even when P is homogeneous: for 
instance, the Bernstein-Sato polynomial of P{xi,X2) = xiX2{xi + ax2) with a 7^ is 
(in our "shifted" notation) s^(s - |)(s + |) Corollary 4.14 and Remark 4.15] [HI 
5.4] @ However, no one that we have consulted seems to have any counterexample to 
Conjecture 19.11 

Let us remark that a necessary condition for a polynomial b{s) to be a Bernstein- 
Sato polynomial is that its roots should be rational numbers < 1: this is the content 
of a famous theorem of Kashiwara [56] [9l Chapter 6] [66t Proposition 2.11]@ Our 
polynomials (19. 2 p and (19. 3p satisfy this condition. 

For the diagonal-parametrized Cayley identities (Theorems 12.101 and 12. lip , a 
slightly more complicated situation arises. The polynomials b{s) arising from the 
basic case / = J = [n] of those theorems, namely 

s{s + 1) ■ ■ ■ {s + n — 1) for an n x n matrix 
b{s) = { _ (9.4) 

s{s + |) ■ • ■ (s + for an n X n symmetric matrix 



are definitely not minimal. Indeed, as remarked already in Section 12. 3[ a lower-order 
Bernstein-Sato pair can be obtained by taking I = J = [n]\{io} for any fixed io € [n]: 

s{s + 1) ■ ■ ■ {s + n — 2) for an n X n matrix 
s{s + |) ■ ■ ■ (s + for an n X n symmetric matrix 



■^^ More generally, this is the Bernstein-Sato polynomial for a homogeneous polynomial P{xi, X2) 
of degree 3 with generic (e.g. random) coefficients [57] . 

This is in our "shifted" notation (|1.3p . In the customary notation, the roots are rational 
numbers that arc strictly negative. See footnote [5] above. 
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We conjecture that these latter polynomials are indeed minimal, but we have no proof 
for general n. 

It is curious that these polynomials are the same as we get for the Laplacian- 
parametrized Cayley identities. Furthermore, also the Q operators corresponding to 
these (conjecturally minimal) polynomials h{s) of degree n — 1 are somewhat similar 
in one respect: namely, for the diagonal-parametrized case the Q operator is given by 
any principal minor of size n — 1 (i.e. I = J = [rt\ \ {iq}) of the relevant matrix Da^p^s 
of differential operators, while for the Laplacian-parametrized case the Q operator 
is a polynomial det(f/ + (9) — det(5) that is a sum over all minors (not necessarily 
principal) of size n — 1. We do not know whether this resemblance is indicative of 
any deeper connection between these identities. 

Note Added: After this paper appeared in preprint form, Nero Budur (private 
communication, July 2011) informed us that he has a proof of our minimality conjec- 
tures for the Laplacian-parametrized and diagonal-parametrized Cayley identities. 



A Grassmann algebra and Gaussian integration 

In this appendix we collect some needed information on Grassmann algebra (= 
exterior algebra) and Gaussian integration (both bosonic and fermionic). We begin 
by recalling the main properties of determinants, permanents, pfaffians and hafni- 
ans (Section lA.ll) . We then recall the well-known properties of "bosonic" Gaussian 
integration, i.e. Gaussian integration over R" or C" (Section IA.2|) . Next we define 



Grassmann algebra (Section IA.3[) and Grassmann-Berezin ("fermionic") integration 
(Section IA.4p . Finally, we explain the formulae for fermionic Gaussian integration 
(Section lA.Sp . which will play a central role in this paper. Our presentation in these 
latter three subsections is strongly indebted to Abdesselam [21 Section 2]; see also 
Zinn- Justin |119[ Chapter 1] for a treatment aimed at physicists. 



A.l Determinants, permanents, pfaffians and hafnians 

Notation: If A is an m x n matrix, then for subsets of indices / C [m] and 
J C \n\ we denote by Ajj the matrix A restricted to rows in / and columns in J, all 
kept in their original order. We also use the shorthand notation Ai,j = A^m] j when 
all the rows are kept, and Aj^, = Aj^n] when all the columns are kept. Finally, if A is 
invertible, we denote by A'"^ the matrix (A^^)"^ = {A^)~^. □ 



A. 1.1 Permanent and determinant 

Let i? be a commutative ring with identity; we shall consider matrices with entries 
in R. In particular, if A = {aij)^j^i is an n x n matrix with entries in R, we define 
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its permanent 

peiA = ^ ai,a(i) ■ ■ ■ (A.l) 

and its determinant 

detA= ^ sgn(a) ai^a(i) ■ ■ ■ an,a{n) ■ (A.2) 

Here the sums range over all permutations a of [n] = {l,...,n}, and sgn((T) = 
(^_X)#(even cycles of 0-) jg ^^iq sign of tlie permutatiou a. See |68j and [79] for basic 
information on permanents and determinants, respectively. 

In this paper we shall need only a few of the most elementary properties of deter- 
minants: 

Lemma A.l (Properties of the determinant) 

(a) detJ = 1. 

(h) det(Afi) = (detA)(det5). 

(c) (Cauchy— Binet formula) More generally, let A be an m x n matrix, and let 
B he an n X m matrix. Then 

dei{AB) = ^ (det A^/)(det E/^) . (A.3) 

/C [n] 
|7| = m 

(d) Let A he an n X n matrix, and define the adjugate matrix adj A hy 

(adjA),,- = {-iy+UetA{,y{,y (A.4) 

(note the transpose hetween the left-hand and right-hand sides). Then 

(adj A) A = A (adj A) = {det A) I . (A. 5) 

In particular, A is invertihle in the ring R"^^"- if and only if det A is invertihle 
in the ring R (when R is a field, this means simply that det A ^ 0); and in this 
case 

A^^ = (detA)-i(adjA) (A.6) 

(Cramer's rule). 

(e) (Jacobi's identity) More generally, if I, J C [n] with \I\ = \ J\ = k, then 

det{{A-^)ij) = {det A)-^e{I, J) {det Ajcjc) (A.7) 
where e{I, J) = (-l)E»e-r*+EieJ-?'. 
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(f) (Multi-row Laplace expansion) For any fixed set of rows I C [n] with 
\I\ = k, we have 

detA = e{I, J) {det Aij) {det Aicjc) . (A.8) 

J C [n] 
\J\ = k 



A. 1.2 Hafnian 



Let A 



define the hafnian [13 



^ij)f'j=i be a 2m x 2m symmetric matrix with entries in R. We then 



where the sum runs over all perfect matchings of the set [2m], i.e. all partitions of the 
set [2m] into m disjoint pairs. There are {2m — 1)!! = {2m)\/{2"^m\) terms in this 
sum. We have, for example, 



hf 



an ai2 

0-12 0,22 



hf 



/ 0,12 0,13 Oi4 \ 

012 022 O23 0-24 

013 O23 034 

y 0.14 024 0.34 0-44 y 



ai2 (A. 10a) 



«i20.34 + 013024 + 014023 (A. 10b) 



Note that the diagonal elements of A play no role in the hafnian. 

Equivalently, we can identify matchings with a subclass of permutations by writing 
each pair ij G M in the order i < j and then writing these pairs in increasing order 
of their first elements: we therefore have 



MA 



(2) ■ ■ ■ 0(j{2m-l)(7(2m) ; 



(A.ll) 



where the sum runs over all permutations cr of [2m] satisfying a{l) < cr(3) < . . . < 
a{2m - 1) and a{2k - 1) < a{2k) for A; = 1, . . . , m. 

Note that if we were to sum over all permutations, we would obtain each term 
in MA exactly 2™'m! times. Therefore, if the ring R contains the rationals, we can 
alternatively write 



MA 



2'"m! 



(2) 



■ O 



(T(2m~l)(7(2m) • 



(A.12) 
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A. 1.3 Pfaffian 



Finally, let A = iaij)tjL, be a 2m X 2m antisymmetric matrix (i.e. = —aji and 
On = 0) with entries in i?o We then define the pfaffian by 

piA = ^ e(M,Mo) n - (^-13) 

MGA<2m (i,j)eM 
j < j 

Here the sum runs once again over all perfect matchings M of the set [2m], and M 
is an (arbitrarily chosen) oriented version of M, i.e. for each unordered pair ij G M 
one chooses an ordering {i,j) of the two elements. The value of the summand in 
(IA.13P will be independent of the choice of M because e(M, Mq) will be odd under 



reorderings of pairs (see below), while A is antisymmetric. Here Mq is some fixed 
oriented perfect matching of [2m] (we call it the "reference matching"). The sign 
e(Mi,M2) is defined as follows: If Mi = {(21,22), (^3, "^4), • • • , {i2m-i,i2m)} and M2 = 
{(ji, J2), (j3, J4), • • • , (j2m-i, J2m)}, then e(Mi, M2) is the sign of the permutation that 
takes ii ■ ■ ■i2m into ji ■ ■ ■ j2m- (This is well-defined, i.e. independent of the order in 
which the ordered pairs of Mi and M2 are written, because interchanging two pairs 
is an even permutation.) This quantity is clearly odd under reorderings of pairs in 
Ml or M2, and has the following properties: 

(a) e(Mi,M2) = e(M2,Mi); 

(b) e(M,M) = +1; 

(c) e(Mi,M2)6(M2,M3) = e(Mi,M3); 

(d) e(Mi, M2) = —1 whenever Mi and M2 differ by reversal of the orientation of a 
single edge; 

(e) e(Mi, M2) = —1 whenever Mi and M2 differ by changing directed edges (a, b), (c, d) 
in Ml to (6,c), {d,a) in M2. 

Indeed, it is not hard to show that e(Mi, M2) is the unique map from pairs of oriented 
perfect matchings into {±1} that has these five properties^ It follows from (a)-(c) 



Such matrices are sometimes called alternating matrices, in order to emphasize that the con- 
dition an = is imposed. This latter condition is a consequence of = —aji whenever R is an 
integral domain of characteristic 7^ 2 (so that 2x = implies x = 0), but not in general otherwise. 
See e.g. [65l section XV. 9]. In this paper we use the term "antisymmetric" to denote 

Qii = 0. 

This also implies that e(Mi,M2) can be given an equivalent (more graph-theoretic) definition 
as follows: Form the union Mi U M2- Ignoring orientations, it is a disjoint union of even- length 
cycles. Looking now at the orientations, let us call a cycle even (resp. odd) if it has an even 
(resp. odd) number of edges pointing in each of the two directions around the cycle. We then have 

e(Ml,M2) = (-l)#(evencycles)_ 
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that the oriented perfect matchings fall into two classes (call them A and B) such that 
e(Mi, M2) equals +1 if Mi and M2 belong to the same class and —1 if they belong to 
different classes. The choice of reference matching Mq really amounts, therefore, to 
choosing one of the two equivalence classes of matchings as the reference class, and 
thereby fixing the sign of the pfaffian. 

The choice of Mq can be encoded in an antisymmetric matrix J defined by 

ri if(^,j)GMo 
^ otherwise 

and satisfying pf (J) = 1. The two most common conventions for the reference match- 
ing are 



Mo = {(1,2),(3,4), (2m-l,2m)}, 



and 



J 



1 




-1 






1 




-1 



V 



(A.15) 



Mo = {(1, m + 1), (2, m + 2), . . . , (m, 2m)} , J = 
In this paper we shall adopt the convention (]A.15|) . We thus have 

= ai2 



/, 



pf 



pf 

0-12 



ai2 
-ai2 



(A.16) 



(A.17a) 



^12034 — 0-13024 + 0-14023 (A. 17b) 



/ ai2 Oi3 ai4 \ 

— O12 0,23 O24 

— ai3 — O23 034 

y — oi4 — 024 ~034 y 

By identifying matchings M G A^2m with permutations a G iSg^ as was done for 
the hafnian, we can equivalently write 



pf A = Sgn(a-o) ^ Sgn(cr) 0,^(1)^(2) ■ • • a^(2m-lM2m) 



(A.18) 



where ctq G iSj^ is the permutation corresponding to the reference matching Mo. [For 
our choice flA.lSp . uo is the identity permutation.] If the ring R contains the rationals, 
we can alternatively write 



pf ^ = sgn(ao) Yl ^g^(^) 



)cr(2) ■ ■ ■ Oo-(2m-l)o-(2m) 



(A.19) 



0"G52t7 



Let us now recall the following basic properties of pfaffians: 
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Lemma A. 2 (Properties of the pfafRan) Let A be an antisymmetric 2m x 2m 
matrix with elements in a commutative ring R. Then: 

(a) pf J = 1. 

(b) (pf A)2 = detA. 

(c) pf(XAX"'") = (detX)(pf A) for any 2m x 2m matrix X. 

(d) (minor summation formula for pfafRans | |53[ , I54] ) More generally, we 
have 

pf(XAXT) = (detX,,)(pfA,,) (A.20) 

/ C [2m] 
|/| = 21 

for any 2i x 2m matrix X (£ < m). Here X^/ denotes the submatrix of X with 
columns I (and all its rows). 

(e) (Jacobi's identity for pfafRans) If A is invertible, then ]ii{A~'^) = (pf A)~^ 
and more generally 

pf{{A~^)H) = e(/)(pfA)-i(pfA,c,.) (A.21) 
for any I C [2m], where e(/) = 

See [SHlllQlllHlEllESlEHllHn] for further information on pfaffians. 

Remark. In this paper we wiU not in fact use the minor summation formula for 
pfaffians; but we will rederive it using Grassmann-Berezin integration. See Theo- 
rem lATTl] and the comments following it. 



A. 2 Bosonic Gaussian integration 

We shall use the following notation: li A = (aij) is an m x n matrix, and I = 
(«!, . . . ,ik) and J = {ji, . . . are sequences of indices (not necessarily distinct or 
ordered) in [m] and [n] , respectively, then we denote by Au the k x £ matrix defined 
by 

iAjj)^ls = ai^j^ . (A.22) 

This generalizes our notation Ajj for subsets I C [m] and J C [n], where a subset 
is identified with the sequence of its elements written in increasing order. We shall 
also use the corresponding notation for vectors: namely, if A = (Aj) is an n-vector 
and / = {ii, . . . ,ik) is a sequence of indices (not necessarily distinct or ordered) in 
[n], then we denote by A/ the fc- vector defined by (A/)q, = A,^. 
Let ip = {ipi)^=i be real variables; we shall write 

V^ = fl^ (A.23) 



100 



for Lebesgue measure on R" with a slightly unconventional normalization. Let A = 
{0'ij)7j=i be a real symmetric positive-definite nxn matrix. We then have the following 
fundamental facts about Gaussian integration on M": 

Theorem A. 3 (Wick's theorem for real bosons) Let A = (aij)"^^^ be a real 
symmetric positive- definite nxn matrix. Then: 

(a) For any vector c = (q)"^]^ in (or C"'), we have 

Jv(pexp{~lip^Aip + c^if) = (det exp(ic^A-^c) . (A.24) 

(b) For any sequence of indices I = [ii, . . . ,ir) in [n], we have 

if r is odd 

(detA)"^/^ hf((A"^)//) ifr is even 

(A.25) 



Jvf(pi^---(pi^exp[-lyD'^Aip) = I 



(c) More generally, for any real or complex r x n matrix C, we have 

if r is odd 

(det M{CA-^C^) ifr is even 

(A.26) 



Historical remarks. Physicists call these formulae "Wick's theorem" because 
Gian-Carlo Wick |115j proved the analogue of ( 1A.25I1 for the correlation functions 
of a free quantum field (see e.g. [97]). These formulae are called "bosonic" because 
the functional-integral formulation for bosonic quantum fields (see e.g. |119] ) leads to 
ordinary integrals over M" or C" (or infinite- dimensional generalizations thereof). By 
contrast, functional integrals for fermionic quantum fields lead to Grassmann-Berezin 
integrals, to be discussed in Sections IA.3l - IA.5l 

The formula flA.25p for the moments of a mean-zero Gaussian measure goes back 
at least to Isserlis j55j in 1918. We thank Malek Abdesselam for drawing our attention 
to this reference. □ 



Now let ip = (v5i)"=i be complex variables; we denote complex conjugation by 
and shall write 

J. J. 

1=1 

for Lebesgue measure on with a slightly unconventional normalization. Let A = 
{(^ij)7j=i be an n X n complex matrix (not necessarily symmetric or hermitian) whose 
hermitian part ^{A + A*) is positive-definite. We then have the following fundamental 
facts about Gaussian integration on C": 
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Theorem A. 4 (Wick's theorem for complex bosons) Let A = {aij)ij=i be an 
n X n complex matrix whose hermitian part \{A + A*) is positive- definite. Then: 

(a) For any vectors b = and c = (q)"^^ in C", we have 

j V{^,^) exp{-^^A^ + b^^ + (^^c) = (detA)-^ exp(FA-^c) . (A.28) 

(b) For any sequences of indices J = (zi, . . . , and J = {ji, ■ ■ ■ ,js) in [n], we have 

if r ^ s 

(detA)-i per((yl"^)/j) if r = s 

(A.29) 



(c) More generally, for any complex r x n matrix B and any complex n x s matrix 
C , we have 



I V{^,^) (fl{B^)c) (fliv^C)^ exp(-^TA^) 



if r ^ s 

{dQiA)-^ ^&i{BA-^C) ifr = s 
(A.30) 



Some final remarks. 1. In this article we shall use mainly the "source" versions 
of Wick's theorem, i.e. part (a) of Theorems IA.3I and IA.4I and the corresponding 
theorems for fermions. The "correlation function" versions, i.e. parts (b) and (c), will 
be used only in Sections H] and 18.11 and in the second proofs of Corollaries 15.21 and 15.41 
(Section [ESD- 

2. We have here presented bosonic Gaussian integration in an analytic context, i.e. 
integration on M" or C". A combinatorial abstraction of bosonic Gaussian integration 
can be found in [T]. 

A. 3 Grassmann algebra 

Let i? be a commutative ring. Every textbook on elementary abstract algebra 
defines the ring . . . , x„] of polynomials in commuting indeterminates Xi, . . . , x„ 

with coefiicients in i?, and studies its properties. Here we would like briefiy to do the 
same for the ring . . . , Xn] Grass of polynomials in anticommuting indeterminates 
Xi, ■ ■ ■ ,Xn- we call this ring the Grassmann algebra over R in generators xi, . . . , x„. 
(Of course, readers familiar with exterior algebra will recognize this as nothing other 
than the exterior algebra A(i?") built from the free i?-module of dimension racj) To 
lighten the notation, we shall henceforth omit the subscripts Grass? since it will always 
be clear by context whether we are referring to the Grassmann algebra or to the 
ordinary polynomial ring. 

Here is the precise definition: 



See e.g. [30l section 6.4] or [65', section XIX. 1] 



36 
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Definition A. 5 (Grassmann algebra) Let R be a commutative ring with iden- 
tity element, and let xi, ■ ■ ■ ,Xn be a collection of letters. The Grassmann algebra 
R[Xi, • • • , Xn] (or R[x] for short) is the quotient of the ring R{xi, ■ ■ ■ , Xn) of noncom- 
mutative polynomials in the letters X\i - ■ ■ iXn by the two-sided ideal generated by the 
expressions XiXj + XjXi A — ^ < J — '^y' ^''^^ xf A — ^ — can consider R[x] 

as a ring and also as an R- algebra. 

In other words, the generators Xi of ^[x] satisfy the anticommutation relations 

XiXj + XjXi = for all i,j G [n] (A.31) 

as well as the relations 

X- = for all ie[n]. (A.32) 

Please note that the anticommutation relation (]A.3ip for i = j states that 2%? = 0; 
but this need not imply xf = if the coefficient ring R does not contain an element 
|. For this reason we have explicitly adjoined the relations = 0. Of course, if the 
coefficient ring R contains an element | (e.g. if it contains the rationals), then this 
extra relation could be replaced by the cases i = j of XiXj + XjXi = 0- 
The first important property of R[x] is the following: 

Proposition A. 6 R[x] is a free R-module with basis given by the 2" monomials 
= Xh-- ■ Xip where I = {ii, . . . ,ip} C [n] with ii < . . . < ip. 

It follows that each element / G -R[x] can be written uniquely in the form 

/C[n] 

with // G R. The term f^ that contains no factors Xi is sometimes termed the body 
of /, and the rest X]/^0 fiX^ is sometimes termed the soul of /. 

Multiplication in the Grassmann algebra is of course i?-bilinear, i.e. 

(j2f^x')(j29JxA = E figjx'x' (A.34) 

^/Cfnl ^ ^JCH ^ /,JC[nl 



where 



, J a{I,J)x''" if/nj = 
XX = S A.35 

if / n J ^ 



Much of the elementary theory works also for coefficient rings without identity element. The 
main change is that polynomials and formal power series $ must have constant term cq in R (and 
not in Z) in order to be applied to elements / G + : see the paragraphs immediately after 
Proposition IA.9I But this means that we can consider the most important case, namely $ = exp, 
only when R has an identity element. For this reason it is convenient simply to make this assumption 
from the beginning. 
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and a{I, J) is the sign of the permutation that rearranges the sequence IJ into in- 
creasing order when / and J are each written in increasing order. 

We define the degree of a monomial in the obvious way, namely, deg(x''^) = 1-^1- 
The Grassmann algebra R[x] then possesses a natural N-grading 

n 

R[X] = ^R[x]p (A.36) 

where R[x]p is generated, as an i?-module, by the monomials of degree p. A coarser 
grading is the Z2-grading 

R[x] = R[xUcn®R[x]odd (A.37) 

where 

^Meven = ^R[x]p (A.38) 

p even 

^Modd = 0i?Mp (A.39) 

p odd 

Note that -R[x]even is a subalgebra of R[x] (but -R[x]odd is not). The parity operator 
P, defined by 

IQn] ^ /C[n] 

is an involutive automorphism of R\x\: it acts as the identity on i?[x]even and as 
minus the identity on -R[x]odd- A nonzero element / that belongs to either -R[x]even or 
-R[x]odd is said to be Z2-homogeneous , and its parity is p{f) =^ in the first case and 

p{f) =^ 1 in the second. It is easy to see that 2'2-homogeneous elements /, g satisfy 
the commutation/anticommutation relations 

fg = {-l)P(fM^)gf . (A.41) 

In other words, odd elements anticommute with other odd elements, while even ele- 
ments commute with homogeneous elements of both types. One easily deduces the 
following two consequences: 

Proposition A. 7 An even element f E R[x\even commutes with the entire Grass- 
mann algebra. (In particular, -R[x]even is a commutative ring.) 

Proposition A. 8 An odd element f G -R[x]odd is nilpotent of order 2, i.e. p = 0@ 
Proof of Proposition IA. 81 Write / = fiX^'^ then 

/odd 



104 



A further consequence of the relations (]A.31I) / (]A.32|) and the finiteness of the 
number n of generators is that every element of 



R[XU = 0^Mp (A.42) 

p=i 

(the set of elements with no term of degree 0) is nilpotent: 

Proposition A. 9 A "pure soul" element f G is nilpotent of order at most 

M = [n/2\ +2, I.e. = 0. 



Proof. Write / = fo±Ji with /o e i?[x]even n R[x]+ and /i e R[x]odd- Then /o 
and /i commute [by flA.4ip ] and = for n >2 [by Proposition IA.8j . so we have 
f = fk + kf^-^fi for all k. And it is easy to see, using f lATT]) / flX32|) . that = 
for k > n/2. □ 

If / = {ii, . . . , -ip} is a subset of [n], the Grassmann algebra R[xi\ '= R[Xii Xip] 
is naturally isomorphic to the subalgebra of R[x\ = R[xiy ■ ■ ■ :Xn] generated by 
{Xi}i€ii we shall identify these two algebras and use the same notation R[xi] for 
both. In particular, the degree-zero subalgebra R[x]o is identified with the coefficient 
ring R. 

If = XlfcLo ^fc^*^ is polynomial in a single indeterminate x with coefficients 
in either R or Z, we can of course apply it to any / G R[x] to obtain $(/) = 
Yl,k=o^kf^ £ R[x\- Moreover, if = Yl^=Q^kX^ a formal power series with 

coefficients in R or Z, we can apply it to any / G -R[x]+ because / is nilpotent and 
the sum is therefore finite. 

But more is true: suppose we have a formal power series = YlT=a^^k/ k^^x'^ 
with coefficients in i? or Z (note the factorial denominators!). If the coefficient ring 
R contains the rationals as a subring (as it is usually convenient to assume), then 
of course $(/) is well-defined for any / G -R[x] + - But we claim that this expression 
has an unambiguous meaning for / G -R[x]+ even if R does not contain the rationals. 
Indeed, let / = Yli^ts fi X^ a-^id consider the expansion of 




/'=|E^^') = E fi.---fi.x''---x''' . (A.43) 

7i,...,47^0 



I odd I, J odd 

I < J 



where I < J denotes that the distinct sets / and J are written in increasing lexicographic order. 
But {x^f = and x^X'' + X''x^ = by (|X3T|) / (|X32)) and (|X4T1) . □ 
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Whenever two or more of the sets Ii, . . . ,Ik are equal (or indeed have any elements 
in common), we have ' " 'X^'° = by (lA.3ip / (lA.32p . When, by contrast, the sets 
/i, . . . , /fc are all distinct, then there are k\ terms with the same coefficient fi^ ■ ■ ■ fi^ 
corresponding to the k\ different permutations of the factors ■ ■ ■ ) X^^i ^ind these 
terms are either all equal (if at most one of the sets la is odd) or else add to zero (if 
two or more of them are odd). It follows that 

at most one \Iol \ odd 

where /i -<...-< denotes that the distinct sets Ji, . . . , are written in increasing 
lexicographic order. So we can define 

n 

fc=o h -<...-< Ik 

at most one \Ia_\ odd 

even if the coefficient ring R does not contain the rationals. 

The most important case is $ = exp. Note that exp(/ + g) = exp(/) exp((yf) 
whenever f,g& R[x\+ ^ -^[xleven, but not in general otherwise. In this paper we will 
apply the exponential only to even elements. 

We also need to define one other type of composition. Suppose ffist that / G 
R{Xi, • • • , Xn) is a noncommutative polynomial in the letters xi, . . . , Xn with coef- 
ficients in R, and let ^i, . . . be elements of some ring R' (not necessarily com- 
mutative) that contains R within its center. Then the composition /(^i, . . . , ^n), 
obtained by substituting each Xi by the corresponding ^j, is a well-defined element of 
R'] furthermore, this composition satisfies the obvious laws 

{f+9m,---,u = /(ei,...,en) + ^7(ei,...,a) (A.46a) 
{fgm,---,^ = /(ei,...,en)^7(ei,...,a) (A.46b) 

Now suppose, instead, that / belongs to the Grassmann algebra R[xi, ■ ■ ■ ,Xn] = 
R{Xi, • • • , Xn) I {XiXj + XjXi^ xf} and that ^i, . . . are elements of R' that satisfy 

e^+Oe. = (A.47a) 

= (A.47b) 

for all Then the composition /(^i, . . . ,^„) is again a well-defined element of R', 
because the relations flA.471) guarantee that any representative of / in R{xi, ■ ■ ■ , Xn) 
will give, after substituting each Xi t>y ^i, the same element of R'; moreover, the 
laws flA.46p continue to hold. In particular, we can take R' to be another Grass- 
mann algebra over R (which may or may not contain some of the Xi &s generators, 
and which may or may not contain additional generators) and ^i, . . . to be arbi- 
trary odd elements of this Grassmann algebra; the relations (]A.47p hold by virtue of 
(lA.4ip and Proposition IA.8| respectively. We shall exploit this type of composition 
in Propositions lA.lOl and IA.12I below. 
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A. 4 Grassmann— Berezin (fermionic) integration 

Thus far we have simply been recalhng standard facts about exterior algebra. But 
now we go on to introduce a process called Grassmann-Berezin integration, which 
has become a standard tool of theoretical physicists over the last 40 years but is still 
surprisingly little known among mathematicians. As we shall see, the term "integra- 
tion" is a misnomer, because the construction is purely algebraic/combinatorial. But 
the term is nevertheless felicitous, because Grassmann-Berezin integration behaves 
in many ways analogously to ordinary integration over or C^, and this analogy is 
heuristically very fruitful. 

We start by defining, for each i G [n], the derivation di = (acting to the right) 
as the i?-linear map di. R[x\ ~^ R[x\ defined by the following action on monomials 
Xii ■ ■ ■ Xip with ii< ...Kip. 



dot I V ^) Xii ' ' ' Xia-lXia+1 ' ' ' Xip if i ia 

^^X^,■■■Xip = < ...... . . (A.48) 

lit ^ {ii, ...,ip\ 



(Of course, in the former instance there is a unique index a G [p] for which i = i^, 
so this definition is unambiguous.) It is then easy to see that flA.48|) holds also when 
the indices ii, . . . ,ip are not necessarily ordered, provided that they are all distinct. 
Clearly di is a map of degree —1, i.e. di: R[x]p R[x]p-u moreover, it takes values in 
i.e. the subalgebra generated by {Xj}j^i- Furthermore, the maps di satisfy 

df = (A.49) 

didj + djdi = (A.50) 
as well as the modified Leibniz rule 

diifg) = m)g + ij'fmg) (A.51) 

where P is the parity operator (]A.40|) . 

We now make a surprising definition: integration is the same as differentiation. 
That is, we define 



dx^ f = dj. (A.52) 

We always write the operator of integration to the left of the integrand, just as we 
do for the (completely equivalent) operator of differentiation. To lighten the notation 
we refrain from repeating the J sign in iterated integrals, so that J dxii ■ ■ ■ dxi^ f is 
a shorthand for J dxii ■ ■ ■ J dxi^ f and hence we have 

j dxh ■ ■ ■ dXik f = di,---dij . (A.53) 
Note by flA.50|) that changing the order of integration changes the sign: 

j dXi^w ■ ■ ■ dXi^ik) f = sgn(a) J dxi^ ■ ■ ■ dxt^ f (A.54) 
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for any permutation cr G iSfc. For instance, reversing the order of integrations gives 

J dx^, ■ ■ ■ rfXn / = (-l)^(^-^)/^ J rfXn " " " dx^, f • (A.55) 

In particular, when we integrate / = J2ic[n] //X^ with respect to all the genera- 
tors, we have 

JdXn ■■■ dxi f = f[n] ■ (A.56) 

That is, integration with respect to all the generators simply picks out the coefficient 
of the "top" monomial I = [n], provided that we write these integrations in reverse 
order^ We adopt the shorthand 

"Px = dxn ■■■ dxi , (A.57) 

and we sometimes write 'Dn{x) if we wish to stress the number of generators. 

An important special case arises when n is even, say n = 2m, and the generators 
Xi, . . . , Xn are divided into two sets ipi, . . . , ip^n and -01, ... , ijjm, where we think of 
each ipi as paired with its corresponding ipi. In this case we adopt the shorthand 
notation 

V{iI),iIj) = dipi dipi ■ ■ ■ dipm dipra (A.58) 

and we sometimes write if we wish to stress the number of pairs of gen- 

erators. Please note that since each pair dipi dipi is Grassmann-even, we can also 
write 

m 

= l[d^id^, (A.59) 

i=l 

where the terms in the product can be taken in any order. Note also that 

= (-l)™('"-i)/2i)^l)^ . (A.60) 

The notation is intended to be suggestive of complex conjugation, but we stress that 
it has nothing to do with complex numbers: it merely denotes the extra combinatorial 
structure on the index set [2m] that arises from the splitting of [2m] into two sets 
of cardinality m and the fixing of a bijection between the two sets@ In particular, 
the coefficient ring R is still completely arbitrary. The general case (x) and special 
case {ipiip) are nevertheless known in the physics literature as "real fermions" and 
"complex fermions", respectively. We shall retain this terminology but shall always 
put the adjectives "real" and "complex" in quotation marks in order to warn the 
reader that they are potentially misleading. 



Of course, since the integrations are performed from right to left, this actuaUy corresponds to 
performing Jdxi first and / dxn last. 

'^^ We can view this as dividing the 2m individuals into m males and m females and then pairing 
those individuals into m heterosexual couples. 
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Now let A = (Aj)rj=i be an n X n matrix, and define new Grassmann variables 
by 6 = Sj=i ^ijXj- We then have 

/„ n 71 

dXn--- dxi 6 ■ ■ ■ = dxn--- ^Xi Yl ^^hXji ■ ■ ■ X] ^"in^i" 
ii=i in=i 

= e(o- ) v4i^(i) ■ ■ ■ AncT(n) 

= detA. (A.61) 
This result can be reformulated as follows: 

Proposition A. 10 (Linear change of variables in Grassmann— Berezin integration) 

Let f e R[^i, ...,^n], and define F e R[xi, • • • , X»i] by 



F{Xl,---,Xn) = /(^l,---,^n) 



(A.62) 



[note that this substitution is well-defined by virtue of the discussion in the last para- 
graph of Section \A.3^ . Then 

Jdxn---dxiF = (detA) JdU---d^,f. (A.63) 

Please note that this is the reverse of the change-of- variables formula for ordinary 
multivariate integrals, i.e. in an ordinary integral the factor det(y4) would appear on 
the other side. 

We will also need: 

Proposition A. 11 (Fubini theorem for Grassmann— Berezin integration) 

Let I = {ii, . . . ,ip} with ii < ■ ■ ■ < ip be a subset of [n], and let P = {ji, . . . ,jn~p} 
with ii < ■ ■ ■ < jn-p- Then for any elements f G R[xi] and g e Rlxi'^] we have 

VxiVxi^ fg = i-if-'^^ (^j Vxi (1 Vxi^ g^ (A.64) 

where Vxi (resp. Vxi'^ ) is shorthand for dxip ■ ■ ■ dxh (resp. dxj„-p ■ ■ ■ dxji)- 

Proof. Expanding / and g in monomials, we see that the only terms contributing 
to the integrals on either side are the "top" monomials x^ and x^" in / and g, respec- 
tively; so we can assume without loss of generality that f = x^ and g = x^^ ■ Now 
use the fact that integration is the same as differentiation, and successively apply 
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the operators di for i G to the product fg using the Leibniz rule (1A.51|) . The 
differentiations hit only g, and we have P/ = (— 1)^/. It follows that 



/ Vxi^ fg = (1 Vxi^ g'j . (A.65) 

The result then follows by integrating both sides with Vxi- n 

Let us remark that the same formula flA.64p would hold for any choice of the 
orderings in defining Vxi and T>xi':, provided only that we use the same orderings 
on both sides of the equation. We have chosen to write Vxi = dxip ■ ■ ■ dxii for 
compatibility with our convention flA.57p that Vx = dxn ■ ■ ■ dxi- 

Finally, the following proposition shows that a Grassmann-Berezin integral over 
dxi is invariant under translation by an arbitrary odd element of the Grassmann alge- 
bra that does not involve the variable Xi- More generally, one can consider integration 
over a set / C [n] of generators: 

Proposition A. 12 (Invariance under translation) Let I = {ii,...,ip} C [n] 
and let Ci, ■ ■ ■ ^ -R[x/'=]odd satisfy = whenever j ^ /. Then 

j dXip ■ ■ ■ dxh fix + ^) = / ^^'p ' ' ' (AM) 

where + denotes the substitution defined at the end of Section \A.^A [Recall that 
the oddness of the C,i is required for this substitution to make sense.] 

Proof. The formula flA.66p can be rewritten as 

■■■d,,f{x + = drp---dnf{x)- (A.67) 

To prove this relation, it suffices to consider the cases in which f{x) = X"^ fo^^ some 
J C [n]. Now for any i G / and j G [n] we have 

9i{Xj+^j) = diXj = Sij (A.68) 

because G R[xi<:]- Using this relation together with the Leibiniz rule, we see that 
di{x + 0"^ equals the same object in which has been replaced by zero. Doing this 
successively for di^^, . . . ,di^, we set to zero for all i E I. But by hypothesis these 
are the only nonzero ^j. □ 
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A. 5 Fermionic Gaussian integration 

We are now ready to state the fundamental formulae ("Wick's theorem") for 
fermionic Gaussian integration, which are analogues of Theorems IA.3I and IA.4I for 
bosonic Gaussian integration. The main difference is that fermionic integration is a 
purely algebraic/combinatorial construction, so that it works over an arbitrary (com- 
mutative) coefficient ring R and does not require any positive-definiteness condition 
on the matrix A. 

We begin with the formula for a pure Gaussian integral, i.e. the integral of the 
exponential of a quadratic form. So let i? be a commutative ring with identity element, 
let xi, ■ ■ ■ ,Xn be the generators of a Grassmann algebra, and let A = (ajj) be an 
n X n antisymmetric matrix (i.e. aij = —aji and an = 0) with entries in i?o We use 
the notation 

Ix^^X = ^ XiaijXj (A.69) 

l<j<j;'<n 

and observe that the right-hand side makes sense even if the coefficient ring R does 
not contain an element |. We then have the following formula, which shows that a 
Gaussian fermionic integral equals a pfaffian: 

Proposition A. 13 (Gaussian integral for "real" fermions) Let A be an n x n 

antisymmetric matrix with coefficients in R. Then 



IT,, / IT,, ( r)f A if n is even 

dXn--- dxi e-^^ = Idxi--- dXn e"^^ = { ^ , , (A.70) 

L ij n is odd 

Proof. We expand the exponential using (IA.45|) with $ = exp and then integrate 
J dxn ■ ■ ■ dxi- When n is odd the integral vanishes, and when n is even (say, n = 2m) 
the only contribution comes from = m in (1A.45I) . yielding 



dXn ■■■ dXl 6 2^^^^^ 



dXn-- - dxi Yl ^'1^1 ■ ■ ■ "^-J™ ■ ■ ■ (A. 71a) 



h ^ ... ^ Im 

Ik = {ikJk} with jfe < jfc 

= ^ e(0") «'(j(1)(t(2) ■ ■ ■ «(T(2m-l)(T(2m) (A. 71b) 

a: a{2i - 1) < o-(2i) and 
o-(l) < ct(3) < ■ ■ ■ < a{2m - 1) 

= pfA. (A.71c) 

(Note that this holds whether or not the coefficient ring R contains the rational 
numbers.) 



See footnote |34] above concerning our use of the term "antisymmetric" . 
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For the variant in which A is replaced by —A and the order of integration is 
reversed, it suffices to observe that 

pf(-A) = (-l)"/2pf(^) (A.72) 
[this is an immediate consequence of the definition flA.18|) ] and use flA.SSP to get 

Jdxi-'-dXne-^''^^'' = (-l)"("-^)/2pf(_^) ^ (-l)"'/2pf(A) = pf(A) 

(A.73) 

because n is even. □ 

Remark. Equation ( 1A.70I) shows that there exist two equally natural conventions 
for Gaussian integrals with "real" fermions: either we write our quadratic forms as 
gax g^Yid our integrals as j dxn ■ ■■ dxi, or alternatively we write our quadratic 
forms as e~2^^^^ and our integrals as J dxi ■ ■ ■ dxn- In this paper we have adopted 
the ffist convention [cf. (lA.57p ]. 

In the special case in which n = 2m and the generators xi, ■ ■ ■ ,Xn are divided into 
two sets ipi, . . . , iprn and ipi, . . . , ipm, we have 

Proposition A. 14 (Gaussian integral for "complex" fermions) Let A be an 

m X m matrix with coefficients in R. Then 

#1 #1 ■ ■ ■ d^jjm d%l)m e^'^^^ = detA. (A.74) 

Please note that here A is an arbitrary matrix; no condition of symmetry or antisym- 
metry need be imposed on it. 

Proof. By the change of variables ip') = {Aip, we have from Proposition lA.lOl 



V{ij,ij)e^ = {detA) J V{tlj,tlj)e'^ (A.75a) 

„ m 

= (det A) J V{^, ^) J](l + ij.ij,) . (A.75b) 



1=1 



When we expand the product, only the term tpiipi ■ ■ ■ipmi^m. has a nonzero integral, 
and the integral of this term is 1. □ 

Here is an alternate proof, which treats "complex" fermions as a special case of 
"real" fermions and invokes Proposition IA.13I 

Alternate proof of Proposition IA.141 Let us write (xi, • • • , X2m) = 
{ipi,. . ., ipm, ^1, • • • , i'm)- Then by (1A.59P we have 

V{il;,ij) = i-l)"" dtpm d^m ■■■ dtfjidipi (A.76a) 
= (_i)™(™+i)/2 d^m--- d^i d^m-- - d^Ji (A.76b) 
= (-l)™("+^)/'rfx2m---rfXi (A.76c) 
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and 

V^^A^ = \x^Kx, (A.77) 
where K is the 2m x 2m matrix 

Then by Lemma [A. 2( c) we have 

viK = det(^° ^ ) •Pf(-/ ) = (-ir(detA) ■ (-ir{™-^)/% (A.79) 



so that 



□ 



jv{'4),^)e^'''^^ = (-l)'"(™+i)/2pf ^ _ ^g^^^ (A.80) 



Let us make one further observation concerning Gaussian integrals (both "real" 
and "complex"), which will be useful in proving certain aspects of Wick's theorem. 
If J = {ii, . . . ,ik) is an arbitrary sequence of indices in [n] — not necessarily distinct 
or ordered — then the "real" Gaussian integral (1A.70P can be generalized to 

1 T zi f pf ^11 if k is even 

dx^, ■■■rfXne^^^^"^^ = r (A.81) 

I if A; is odd 

To prove this, it suffices to observe, first of all, that if the indices ii, . . . ,ik are distinct, 
then this formula is merely flA.70p after a relabeling of indices: the point is that by our 



definition flA.22p . the variables arise in the same order in the integration measure and 
in the matrix Ajj. On the other hand, if the sequence ii, . . . ,1^ contains any repeated 
index, then the left-hand side vanishes because dxi d-Xi = by flA.49|) , while pf An = 
in such a case because the pfaffian changes sign under simultaneous permutations of 
rows and columns [this is a special case of Lemma lA.2( c)]. In a similar way, the 
"complex" Gaussian integral (lA.74p can be generalized to 

d^lJi, dtPj, ■ ■ ■ d^Pi, d^Pj, e^^^''^' = det Aji (A.82) 

for arbitrary sequences I = {ii, . . . ,ik) and J = (ji, . . . ,jk)- 

We are now ready to state the full Wick's theorem for fermions. Since the "source" 
version of Wick's theorem for fermions will involve a fermionic source A [cf. flA.83P ] 



this means that we will be working (at least when discussing this equation) in an 
extended Grassmann algebra R[x, 0] with generators Xi? • • • ? Xn and 9i, . . . ,9]\f for 
some > 1, and the sources Aj will belong to the odd part of the Grassmann 
stx^algebra R[9]. 

We assume without further ado that n is even, i.e. n = 2m. Recall also from 
(IA.22P the notation Ajj for arbitrary sequences of indices I and J. 
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Theorem A. 15 (Wick's theorem for "real" fermions) Let R be a commutative 
ring with identity element, and let A = {aij)fjLi be a 2m x 2m antisymmetric matrix 
with elements in R. Then: 

(a) If the matrix A is invertible, we have 

I Vx exp(ix^^X + A^x) = (pf ^) exp(iATA-iA) (A.83) 

whenever A = {Xi)f^i are odd elements of the Grassmann algebra that do not 
involve x, ^-C- Aj G -R[6']odd- 

(b) For any subset I = {ii, . . . ,ir} C [2m] with ii < . . . < ir, we have 

/( if r is odd 

'DxXn■■■X^reMlx'^x) = \ . (A.84) 

pt Ajcjc if r IS even 

(c) For any sequence of indices I = {ii, . . . ,ir) in [2m], if the matrix A is invertible 
we have 

if r is odd 

(pf y4) pf ((y4~'^)//) if r is even 

(d) More generally, for any r x 2m matrix C with entries in R, we have 

/( Q if r is odd 

Vx{Cx)i--- {Cx^e^vihx'Ax) = I ^ (deta/)e(/) pf(v4,e,c) if r is even 

[ U\=r 

(A.86) 

where Cj,/ denotes the submatrix of C with columns in I ; and if the matrix A 
is invertible. 



jl^XXh--- Xir ewikx^Ax) 



jvx{Cx)i---{Cx)rew{\x^Ax) = I 



if r is odd 

(j)f A) pf(CA-^C^) if r is even 

(A.87) 



Note that (1A.84P and (1A.85P are equivalent by virtue of Jacobi's identity for pfaf- 
fians [cf. flA.2ip ]. and that (IA.86P and flA.87p are equivalent by virtue of Jacobi's 
identity and the minor summation formula for pfaffians [cf. flA.20p ]. However, we 
shall give independent proofs of all these formulae; as a consequence, our argument 
provides ab initio proofs of these pfaffian identities by means of Grassmann-Berezin 
integration. 

Proof of Theorem IA.15[ To prove (a), we perform the translation x' = X + A^^X 
and use Proposition IA.12I We have 

Ix'^Ax' + A^x' = Ux^ + X^A-^)A{x + A^'X) + A^lx + A^'X) 

= \x^Ax + \X^A-^X (A.88) 
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since A"""" = —A^^ and = ~^^X- Therefore ( 1A.83P is an immediate consequence 
of Proposition IA.13I 

For (b)-(d), let us first remark that if r is odd, then the integral must vanish, 
because the expansion of the exponential can only give an even number of factors of 
X, so the "top" monomial in x cannot be generated (recall that the total number of 
generators is even, i.e. n = 2m). 

To prove (b), we compute 

= e(^) jl^Xi1^Xi-X^ exp{lxJcAi^icXic) 

= e{I) ) (^j Vxi x') (1 Vxi^ e^v{\xlA,.,.xi'^) 

= e(/) pf(A,c,c) . (A.89) 

The first equality holds because any factor with i G / arising from the expansion 
of the exponential would be annihilated by the prefactor x^'i the second equality is 
simply a reordering of the integration variables; the third equality is Fubini's theorem 
(Proposition lA. 1 : and the last equality is simply the evaluation of the two integrals 



(using Proposition IA.13I for the second) together with the fact that r is even. 

(c) We specialize result (a) to the case where the "sources" Aj are generators of 
an extended Grassmann algebra, then differentiate (or equivalently i nteg rate!) with 
respect to the sources Ajj, . . . , A,^, and finally replace all the A, by zercCj: 

hxXn■■■X^r^M\x^^x) = (Pf^) T^T 7^ exp(iATA-iA) 

J (J^ii (J^ir A=0 

= (pfA) jdX,,---d\^exp{'^Xj{A-')jjXj) 

= (pfA) {-ly^^-'y^ jvXj exp{lXj{A-')uXi) 
= (pfA)(-l)^(^-i)/2pf((A-i),,) 

= (pf A) (-l)-(-l)/2 (_l)^/2 pf((_A-l),,) 

= (pfA)pf((A-T),,) (A.90) 

Here the third equality involved reordering the integration variables from increasing 
to decreasing order; the fourth equality performed the Gaussian integral using (1A.81I) : 
the fifth equality used (1A.72I) : and the final equality used —A~'^ = A~^ and (—1)^ = 
1 (which holds since r is even). 



Setting A — can be interpreted as extracting the monomials that do not involve A, or 
alternatively as a special case of the substitution discussed at the end of Section lA.3l (since is odd). 
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(d) We have 

Jvx {Cx)i---{Cx)r exp{lx^Ax) 

= J2 ^1^1 ■ ■ ■ jT^XXn--- X^r exp(ix^^x) , (A.91) 

but the only nonvanishing contributions in the sum come when i all dis- 

tinct; so we can first require zi < . . . < v and then sum over permutations, yielding 

E E^iMi)---^".M sg^(^) jvxx'^MlX^^x) = 5^(deta/)e(/)pf(A,.,0 

\I\=r cr£Sr \1\=^ 

(A.92) 

by part (b). If A is invertible, we can write 

Jvx{Cx)i ■■■ {Cx)r expilx^Ax) = ^■■■^ j Vx exp{\x^ Ax + \^Cx) 



21 

d d 



A=0 



(pfA) exp{^\'^CA-'C^\)\^^^ 



(A.93) 



by using (a). Then by the same reasoning as in (1A.90I) we can see that this equals 
(pf A) pf(CA-TcT). □ 

Next we state Wick's theorem for "complex" fermions. Once again, when dis- 
cussing the "source" version of this theorem [cf. flA.94p ]. we will work in an extended 
Grassmann algebra R[ip,ip,9] in which the sources Aj, Aj belong to the odd part of 
the Grassmann subalgebra R[6]. 

Theorem A. 16 (Wick's theorem for "complex" fermions) Let R be a commu- 
tative ring with identity element, and let A = {aij)^j^-^^ be anxn matrix with elements 
in R. Then: 

(a) If the matrix A is invertible, we have 

Jv{^P,^P) exp{^P^A^ + Ft/; + t/^^A) = (det A) exp{-\^A-^X) (A.94) 

whenever X = (Aj)"^;^ andX = (Aj)"^| are odd elements of the Grassmann algebra 
that do not involve ip and ip, i.e. Aj, Aj G -R[6']odd- 

(b) For any subsets I = {ii, . . . ,ir} and J = {ji, . . . ,jr} of [n] having the same 
cardinality r, with ii < . . . < ir and ji < . . . < jr, we have 

[v{i:,^p) (fl^Pi^iJj^ expi^p^Atp) = e{I , J) {det Aj.j.) . (A.95) 
[If there is an unequal number of factors ip and ip, then the integral is zero.] 
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(c) For any sequences of indices I = (ii, . . . ,v) and J = (ji, . . . , jV) in [n] of the 
same length r, if the matrix A is invertible we have 

jvi^,^) (n^.>,,) exp(^TA7A) = (det A) det((A-T),j) . (A.96) 

[Again, if there is an unequal number of factors and ip, then the integral is 
zero.] 

(d) More generally, for any r x n matrix B and n x r matrix C with entries in R, 
we have 

[v{tlj,ij) (f[ {ijC)jBij)}j exp(^TA^) 

= <I, J) {det B,j) {det Ai.jc) {det C I,) , (A.97) 

\I\ = \J\=r 

and if the matrix A is invertible, 

fv{ij,^p)(f[{ijC)jBij)^ exp{i:^Aij) = {det A) det {BA-'C) . (A.98) 

Note that (IA.95P and flA.96P are equivalent by virtue of Jacobi's identity ( IA.7p . 
and that flA.97p and flA.98P are equivalent by virtue of Jacobi's identity together with 
the Cauchy-Binet identity flA.Sp . However, we shall give independent proofs of all 
these formulae; as a consequence, our argument provides an ab initio proof of Jacobi's 
identity by means of Grassmann-Berezin integration. 

The proof of Theorem IA.16I follows closely the pattern used in the proof of Theo- 
rem lATTSl but with slightly different combinatorics. 

Proof of Theorem IA.16[ To prove (a), we perform the translations tp' = ip—A^^X 
and ip' = ip — A^^A and use Proposition IA.12I We have 

ij''^A^p' + A^V' + V-'^A 

= - \^A-^)A{i, - A-^\) + \^{,p - A-^\) + (V-^ - \^A-^)\ 

= ij^AiP - X^A'^X . (A.99) 

Therefore (IA.94P is an immediate consequence of Proposition IA.14I 
To prove (b), we begin by observing as before that 

^^(^,^)(riV^.>j.) exp(V^T^V^) = !v{ij,ij)(T{i,,^ijAexv{i'lAi.jci;j.). 

(A.lOO) 
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We now proceed to reorder the integration measure. We have / = {ii, . . . ,ir} with 
ii < . . . < ir, and let us write P = . . . , in_r} with i[ < . . . < and hkewise 
for J. Then 'D{ijj,'ijj) is the product of factors dipidipi taken in arbitrary order; we 
choose the order to be IP (i.e., ii ■ ■ -iri'i ■ ■ -i'n-r)- '^o"^ leave the factors ip in 
place, but reorder the factors if) to be in order JJ'^ rather than IP: this produces a 
sign e(/, J). Using the notation 

V{ipj,iJi) = dipj^d^pi^ ■ ■ ■ dtpj^d^Ji^ (A. 101) 

and likewise for the complementary sets, we have proven that 

V{^P, ^P) = e(/, J) V{^j, ^j) V{^j., V'/O • (A.102) 

This in turn can be trivially rewritten as e(/, J) ^{ipjc^ipic) Vlipjji/ji) since the two 
factors V are Grassmann-even. We can now apply this measure to the integrand 
in flA.lOOp : integrating the monomial against Vi^pj^ipi) gives 1, and integrating the 
exponential against V^ipjc^ijjjc) gives detAjcjc by Proposition lA. 141 

(c) As in the case of "real" fermions, we specialize result (a) to the case where the 
"sources" Aj, Aj are generators of an extended Grassmann algebra, and then differen- 
tiate with respect to them: 



A=A=0 



= (detA) i-iy jvi\j,Xj)exp{-\^j{A-')jj\i) 
= (detA)(-l)'- detii- A-') ji) 

= (detA) det{{A~^)ij) . (A.103) 

Here the minus sign in the first equality comes from the fact that differentiation of 
the source term ip^X with respect to A, yields —ipi according to flA.Sip : the second 
equality says that differentiation is the same as integration; and the third equality 
uses flXMj) . 
(d) We have 



= E E (n^^-"^".i'^) M^'^) (n^^^^^^^ (A.104) 
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but the only nonvanishing contributions in the sum come when ii,. . . ,ir, and also 
ji, . . . ,jr, are all distinct; so we can first require ii < . . . < ir, and ji < . . . < jr and 
then sum over permutations, yielding 



^v(i),i ■ ■ ■ ^ia(r),r -^i^^ci) ■ ■ ■ ^^.j^^,.) 

\I\=r \J\=r o",TgiSr 



X 



sgn((T) sgn(r) / (11^*-^^- ) exp(^^A^) 



\I\=r \J\=r 



(A. 105) 



by part (b). If A is invertible, we can write 

f V{^,^P) (f[ (^C)^ {Bij)^ exp(^TA^) 

r'(^,^) exp(^^A^ + X^BiP + iP^CX) 



n 

r 

n 



9 9 



dXi dXi 

d_d_ 

dXi dXi 



A=A=0 



(detA) exp{-X^BA-^CX) 



A=A=0 



(A. 106) 



by using (a). Then by the same reasoning as in (1A.103P we see that this equals 
{detA)det{BA-^C). □ 



A. 6 Summary of bosonic and fermionic Gaussian integration 

The four types of Gaussian integration and their basic characteristics are summa- 
rized in Table [TJ Let us stress the following facts: 

• Fermionic integration is a purely algebraic/combinatorial operation: no analytic 
conditions on the matrix A (such as positive-definiteness) are needed, and in fact 
fermionic integration makes sense over an arbitrary (commutative) coefficient 
ring R. Bosonic integration, by contrast, is an analytic operation (at least as 
we have defined it here) and requires positive-definiteness of A or its hermitian 
part. 

• For complex bosons and "complex" fermions, no symmetry or antisymmetry 
conditions are imposed on the matrix A. By contrast, for real bosons and "real" 
fermions, the matrix A must be symmetric or antisymmetric, respectively. 
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Type of 
Variables 


Combinatorial 
conditions 


Analytic 
conditions 


Value of 
Gaussian int. 


Contractions in 
2fc-point correlation 


Real bosons 


A symmetric 


A pos.-def. 


(det.4)-i/2 


2k X 2k hafnian 


Complex bosons 


none 


A + A* pos.-def. 


(det A)-' 


k X k permanent 


"Real" fermions 


A antisymm. 


none 


pf^ 


2k X 2k pfaffian 


"Complex" fermions 


none 


none 


det^ 


k X k determinant 



Table 1: Summary of the four types of Gaussian integration. Here A is the matrix 
appearing in the quadratic form in the exponential. See text for details. 



B Some useful identities 

In this appendix we collect some auxiliary results that will be used at various 
places in this paper: identities for sums of products of binomial coefficients (Sec- 
tion [BT]) and for determinants and pfaffians (Section IB. 2p . lemmas on matrix factor- 
ization (Section IB.Sp . and an identity that we call the "dilation-translation formula" 
(Section El. 

B.l Binomial identities 

We collect here a few identities for sums of products of binomial coefficients that 
will be needed in the proof of the rectangular Cayley identities (Sections 15. 5l and 15. 6p . 
We use the standard convention [17] for the definition of binomial coefficients: 




[r — k + l)/k\ for integer k > 

forfc = (B.l) 
for integer k < 



where r is an indeterminate and k is always an integer. Multinomial coefficients 

'ai H hafc\ (ai H hOfc)! 



B.2) 



will, by contrast, be used only when all the are nonnegative integers. 

We now state some easy combinatorial lemmas involving binomial coefficients. 
They are all either contained in the textbook of Graham-Knuth-Patashnik ^47^ Chap- 
ter 5] or derived in the same fashion. 

Lemma B.l (parallel summation) 

" \Z^I form integer. (B.3) 

k<m V / \ / 
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Lemma B.2 (Chu— Vandermonde convolution) 



Efw\ f m \ fw + m\ 

Proof. For positive integers w and m, this is the number of ways of selecting p 
people out of a group of w women and m men. Since, for any fixed integer p, both 
sides are polynomials in w and m, the identity holds as a polynomial identity. □ 



Lemma B.3 



k,h,l>0 
k + h + I = m 



Proof. The left-hand side equals 




where the first and third equalities use the Chu- Vandermonde convolution. □ 



Lemma B.4 



E 



k,h,l>0 
k + h + I < m 



h } \ k } \ h } \ m 



Proof. Sum flR5l) and use flRSl) . □ 
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B.2 Determinant and pfafRan identities 



We collect here some identities for determinants and pfaffians that will be needed 
in Section El these concern the determinant or pfaffian of a partitioned matrix and 
the change of determinant under low-rank perturbation. These identities are well- 
known, but for completeness we will give compact proofs using Grassmann-Berezin 
integration. We will also give a (possibly new) "fermionic" analogue of the low-rank- 
perturbation formula, which will play a crucial role throughout Section O 

Let us begin with the formula for the determinant of a partitioned matrix, due to 
Schur [HHl Hilfssatz, pp. 216-217]: 

Proposition B.5 (Schur's formula for the determinant of a pcirtitioned matrix) 

Consider a partitioned matrix of the form 



M 



A 


B 


C 


D 



(B.8) 



where A, B, C, D are matrices of sizes mxm,mxn,nxm and n x n, respectively, 
with elements in a commutative ring with identity. 

(a) If A IS mvertible, then det M = {det A) det{D - CA-^B). 

(h) If D is invertible, then det M = {det D) det{A - BD-^C). 

The matrix D — CA~^B is called the Schur complement of A in the partitioned 
matrix M; see pH[3T|[76llll8j for reviews. One well-known proof of Schur's formula 
is based on the identity 



(B.9) 



Let us give a quick proof of Schur's formula using Grassmann-Berezin integration: 



A 


B 


C 


D 



A 




( 


A-^B 


C 


In J 


V 


D - CA-^B 



Proof of Proposition IB. 51 Let R be the commutative ring with identity in 
which the elements of A,B,C,D take values. We introduce Grassmann variables 
"^i^i^i (1 ^ ^ ^ ^) ciiid rij,f]j (1 < j < n) and work in the Grassmann algebra 
R[i), ijj, rj, f]]. We have 



det 



A 


B 


C 


D 



Vm{i',i^)'Dn{T],T]) CXp 



A 


B 


C 


D 



Vm{ip, ip) Vn{r], f]) exp [ip^Aip + ^j'^Bt] + fC^j + f Dt]] . 

(B.IO) 

If A is invertible, we can perform the integration over t/',?/^ using Wick's theorem for 
"complex" fermions (Theorem IA.16p . yielding 

(det A) jVniv,v) ewiv^Dri-f^^CA-^Bri) . (B.ll) 
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Then performing the integration over rj, i] yields (det A) det{D — CA ^B). This proves 
(a); and the proof of (b) is identical. □ 

Here are some important special cases of Proposition IB. 51 

Corollary B.6 Let U,V be mxn matrices with elements in a commutative ring with 
identity. Then 

detiUV^) = det ( S^. f ] , (B.12) 



where 0^ is the m x m zero matrix and In is the n x n identity matrix. 

Corollary B.7 (matrix determinant lemma) Let A be an invertible mx m ma- 
trix, let W be an invertible n x n matrix, and let U,V be m x n matrices, all with 
elements in a commutative ring with identity. Then 

det(A + UWV^) = (det A) (det W) det(W^"^ + V^A'^U) . (B.13) 

In particular, if we take W = In, then 

det (A + UV^) = (det A) det(/„ + V^A'^U) . (B.14) 

// in addition we take A = Im, then 

det{In, + UV^) = det{In + V^U). (B.15) 

Corollary IB. 71 is sometimes known as the "matrix determinant lemma" , and the 
special case flB.lSP is sometimes known as "Sylvester's theorem for determinants". 
When n ^ m, we can interpret ( lB.13p -( fB.15l) as formulae for the change of determi- 
nant under a low-rank perturbation: see Lemma TB . 1 1 1 below for an explicit statement. 

Analogues of Proposition IB. 51 and Corollary IB. 61 exist also for pfaffians: 

Proposition B.8 (PfafRan of a partitioned matrix) Consider a partitioned ma- 
trix of the form 



M 



A 


B 




D 



(B.16) 



where A, 5, D are matrices of sizes 2m x 2m, 2m x 2n and 2n x 2n, respectively, with 
elements in a commutative ring with identity, and A and D are antisymmetric. 

(a) If A IS invertible, then pf M = {pi A) pf(D + B^A'^B). 

(b) If D IS invertible, then pi M = {pi D)pi{A + BD'^ B"^) . 

Note that the matrices D + B^A^^B and A + BD^^B^ appearing here are just the 
usual Schur complements D — CA^^B and A — BD^^C specialized to C = —B^. 
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Corollary B.9 Let U be a 2m x 2n matrix with elements in a commutative ring with 
identity. Then 

/ n TT \ 

(B.17) 



pf(t/J2nf/^) = (-l)'"pf(-f/J2nf/^) = (-l)™pf 

where is the standard 2n x 2n symplectic form /IA.15\) . 



02 m 


u 




J2n 



Proof of Proposition IB.81 This time we introduce "real" Grassmann variables 
9i {l<i< 2m) and \i {I < i < 2n). We have 

pfM = jv2rn{0)V2n{\) e^v[\O^Ad + d^B\ + \X^D\\ . (B.18) 

If A is invertible, we can perform the integration over 6 using Wick's theorem for 
"real" fermions (Theorem IA.15I) . yielding 

(pf A) j V2n{X) exp (|A^DA + IX^B^A-^BX) . (B.19) 

Then performing the integration over A yields {pi A)pi{D + B'^A'^B). This proves 
(a); and the proof of (b) is identical. □ 



We next wish to prove an analogue of Corollary IB. 71 when the entries in the 
various matrices belong, not to a commutative ring, but to a Grassmann algebra. 
More precisely, the entries in A and W will be even elements of the Grassmann 
algebra, while the entries in U and V will be odd elements of the Grassmann algebra. 

Proposition B.IO (fermionic matrix determinant lemma) Let Q be a Grass- 
mann algebra over a commutative ring with identity; let A be an invertible m x m 
matrix and W an invertible n x n matrix, whose elements belong to Qeweni and let 
U,V be m X n matrices whose elements belong to QoAd- 

det(A + UWV^) = (det A) (det W)-^ det(l^-^ + V^A-^U)-^ . (B.20) 

In particular, if we take W = In, then 

det{A + UV^) = (det A) det(4 + V^A-^U)-^ . (B.21) 

If in addition we take A = Im, then 

det{Im + UV^) = det(4 + V^U)~^ . (B.22) 



^■^ Note that all the matrix elements in all these determinants belong to the commutative ring 
^cvon; therefore, these determinants are unambiguously defined. 
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Note that (]B.20p - (]B.22l) differ from (]B.13|) - (]B.15P by replacing some determinants 
witli tlieir inverse. 

Proof. Let us first observe that since A and W are invertible, we can rewrite f lB.20p 

as 

det{Im + A-^UWV"^) = detiln + WV^A-^uy\ (B.23) 
which is equivalent to 

det{Im + UV^) = det{In + V^U)~^ (B.24) 

under the (invertible) change of variables U = A'^U, V = VW"^. So it suffices to 
prove (Ib:22|) . 

Let us first prove (1B.22P when the ring R is M. We augment the Grassmann algebra 
by introducing generators f]j,f]j (1 < J < '^), and we also introduce bosonic variables 
(1 < < fn)- We wish to consider the mixed bosonic-fermionic integral 

j Vmiv. V) T^niv. n) exp [-^""^ + ^^Ut] + fv""^ + fr]] . (B.25) 

Please note that the quantities f/r], ff'V'^ and ff-rj all belong to the even "pure soul" 
part of the augmented Grassmann algebra, and in particular are nilpotent; therefore 
the integrand can be interpreted as 

oo ^ 

exp(-^Typ) - i^^Uv + fV^V + fvf (B.26) 

fc=0 

where the sum over k is in fact finite. The Gaussian integration over 99, (p can thus be 
interpreted as separate Gaussian integrations for the coefficients (which belong to M) 
of each monomial in the augmented Grassmann algebra; this makes perfect analytic 
sense. Let us now evaluate (IB.25P in two ways: Integrating first over 99, ip and then 
over 77, 77, we get 

j Vniv, V) exp [fri + f]'^V^Uri] = det(4 + V^U) . (B.27) 

On the other hand, integrating first over rj, fj and then over (f, (f, we get 

j V^{ip, if) exp [-^^yp - ip^UV^ip] = det(/^ + UV^)'' . (B.28) 

This proves (]B.22|) when the ring R is M. 

Let us finally give an abstract argument showing that (]B.22p holds for arbitrary 



commutative rings R. The point is that the matrix elements of U and V belong to 
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a Grassmann algebra over some finite set of generators xi, . . . , xn-, and fience can be 
written as 

Uij = '^^r,Kx'' (B.29a) 

Kodd 

"^^3 = E f^^r.Kx'' (B.29b) 

Kodd 



for some coefficients aij-x, Pij;K ^ R- Now, botli sides of (IB.22|) are of tlie form 



^ IlX^ wliere tlie coefficients 7^ are polynomials in {oiij-^K^ Pij;K} with integer 



L even 



coefficientso But we have just shown that these two polynomials coincide whenever 
{aij-^Ki Pij;K} are replaced by any set of specific values in M. Therefore, they must 
coincide as polynomials in the indeterminates {dij-^K, Pij;K}- But this implies that 
they are equal when {dij-K, Pij;K} are replaced by specific values in any commutative 
ring R. □ 

It is convenient, for our applications, to rephrase (IB.lSp and (1B.22P as explicit for- 
mulae for the change of determinant under a low-rank perturbation. So let Ui, . . . , m„ 
and Vi, . . . , w„ be vectors of length m whose entries are elements of a Grassmann al- 
gebra and are either all Grassmann-even or all Grassmann-odd. We call these cases 
e = +1 and e = — 1, respectively. We then have the following formula for the deter- 
minant of a rank-ra perturbation of the identity matrix: 

Lemma B.ll (low-rank perturbation lemma) Let ui, . . . ,Un and vi, . . . ,Vn be 

as above, and define the n x n matrix M by 

m 

Map = Va-Up = ^(t'a)i(M/3)i . (B.30) 



Then we have 



n 

det + J2 = det(/„ + MY (B.31) 



0=1 



where e = ±1 is as above. 



Here is one special case of the low-rank perturbation lemma that will be useful in 
Sections 15.31 and 17.21 in treating symmetric Cayley identities: 

Corollary B.12 Let rj = (771, . . . , rjm) and fj = (f/i, . . . , fjm) be Grassmann variables, 
and let A be an invertible m x m symmetric matrix whose elements are Grassmann- 
even (hence commute with everything). Then 

det(J^ + r/r/^- A-^r/r/^A) = (l-r/^r])-^. (B.32) 



The determinant on the right-hand side is an element of the Grassmann algebra whose "body" 
term is 1; therefore, it is invertible in the Grassmann algebra and the coefficients of its inverse are 
polynomials (with integer coefficients) in its own coefficients. 
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Proof. Applying Lemma [B.lll with ui = t], vi = rj, U2 = A ^rj and V2 = —A^t] 
—Afj gives 



det(J. + - A-'m^A) = det-^ ( ff^l"^^ j (B.33) 

Since A and A~^ are symmetric, the off-diagonal elements vanish, which gives the 
result. □ 

Another special case of the low-rank perturbation lemma will arise in Section 15.91 
here the m x m matrix + UV"^ occurring on the left-hand side on (]B.31[) will be 
written as a product of rectangular matrices, each of which is a rank-one perturbation 
of the corresponding rectangular pseudo-identity matrix. (The mxn pseudo-identity 
matrix has matrix elements {Imn)ij = ^ij-) Direct application of the low-rank 
perturbation lemma to such a product matrix yields a rather messy result, but after 
some row operations we can obtain a fairly neat alternative formula: 

Corollary B.13 Fix integers i > 1 and ni, . . . ,n£ > 1 with > ni for 1 < a < I, 
and write ra^+i = ui. Let Xi, . . . ,Xi and yi, . . . ,yi be vectors, where Xa is of length 
and ya is of length Ua+i, whose entries are elements of a Grassmann algebra and are 
either all commuting (e = +1) or else all anticommuting (t = —1)- Then we have 

det(l[{Tn^n^^,-x^yl)] = {detNY (B.34) 

where the product is read from left (a = 1) to right (a = i), and the £ x i matrix N 
is defined by 

_E yfx^. if a < /3 
N^p = I (B.35) 
Sai3 - E vT^i if a > /3 

i=l 

where ria^p = min n^. 

Proof. Note first that /„„„„+i/„^+i„^+2 ■ ■ ■ 4^-in/3 is an x matrix whose ij 
element is 1 if z = j < Ua^js and otherwise. In particular, /„iri2^n2n3 ' ' ' ^na^ma = 
IniHa- For a vector of length m > ni, define v as the vector restricted to the first 
ui components. So we can expand the matrix on the left-hand side of (1B.34P as 



a=l 



a=l 
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a + 2 



(B.36) 



a=l 



In this form, we are ready to apply Lemma IB.llI with vectors = ~Xa and 



a + 2 



a+1) 



{Intm+i - xeyj) for a = 1, . . . , This gives 



det(l[{Tn^r.a+,-x^yl)] = (detiVW)^ 



(B.37) 



with N^p = Sai3 + v^Ui3. Now observe that 

„,T -T 



I3>a 

"q+1,/3 „,a™/9 



(B.38) 



and call Ca[s = 2/J/n<,+in<,+2 " " " ^g.m^a;^ = T^iZV'" y?x'^ ■ Thus, defining the i x 
upper-triangular matrix C as 





if a 


<f3 




if a 


= P 





if a 


>(3 



we have 



(B.39) 



(B.40) 



Clearly det (7 = 1, so if we define = CA^(o) Yiave det iV = det N^°\ It is easy to 
see, using ( 1B.40p . that N is exactly the matrix given in (lB.35p . □ 



Please note that all the entries of the matrix N are polynomials of degree at 
most two in the variables x and y — unlike the matrix A^^"^ coming from the bare 
application of the low-rank perturbation lemma, which contains terms of degree as 
high as 2i. 



Corollary IB. 61 is in fact the case £ = 2 of a more general lemma that holds for 
i > 2, and will be needed in Section 15.91 

Lemma B.14 Fix integers i > 2 and ni,...,n£ > 1, and write ng^i = ni. Let 
Ui, . . . ,Ui be matrices with elements in a commutative ring with identity, Ua being of 
dimension x Ua+i ■ Define 



(B.41) 











■■■ 





\ 











-U2 









.,Ue) := 








Ins 
























V Ue 





... 




/ 
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where 0„ is the n x n zero matrix and In is the n x n identity matrix. Then 

det Me{Uu..., Ue) = det(f/i ■■■Ue) . (B.42) 

Proof. We prove this by induction on i. The case £ = 2 is aheady proven by 
Corollary IB.6I For £ > 2 we use the Grassmann representation of det(M^): 

detM,(f/i,...,f/,) = jVnM\i,^)---Vn,{il^',i;') 



X exp 



a=2 



[B.43) 



Now perform the integration over (tp^jip^): highlighting the factors in the integrand 
that involve these fields, we see that 



Vn,{ip',i>')exp[i^'^'-^'-'U,.,^' + ^%^'] = exp(^^-if/,_if/,V^i) . (B.44) 

Comparing this with ( ]B.43p . we see that 

detM,(f/i,...,?7,) = detMe.i{Uu...,Ue_2,Ue-iU,) , (B.45) 
which provides the required inductive step. □ 



Remark. Just as Corollary IB. 61 is a specialization of the more general Proposi- 
tion IB.5t so Lemma IB. 141 has a similar generalization, proven through an identical 
procedure (of which the details are left to the reader), namely: 

Lemma B.15 Fix integers i > 2 and ni,...,ne > 1, and write ra^+i = ni. Let 
Bi, . . . ,B£ be matrices with elements in a commutative ring with identity, Ba being 
of dimension Ua x ria+i. Let Ai,...,Ae be square matrices with elements in the 
same commutative ring, being of dimension n^x ria- Assume that A2, . . . ,Ai are 
invertible. Define 



M{Ai,...,Af,Bi, 



Bf) :-- 



1 ^1 




••• 










A2 


-B2 

























-Be- 


-1 


V B, 





... 


A, 


/ 



(B.46) 



Then 



detM(Ai,...,Af;5i,...,E^) = det(Ai + SiAg ■ ■ ■ ^7^5^) JJdetAj. (B.47) 

i=2 
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// also Ai is invertible, then we can obtain an expression with a form of cyclic sym- 
metry: 

i 

det M{Ai, ...,Ae]Bi,...,Bg) = det(/„, + A^^^iA^^^a ■ ■ ■ Aj'^Be) JJ det Aj . 

i=i 

(B.48) 

B.3 Matrix factorization lemmas 

In Sections 15.41 15.51 15. 6^ 15.71 and 15.91 we shall need some matrix factorization 
lemmas having the general form 

For any matrix X of the form there exists a matrix A of the form 

such that A) = [where $ denotes a specified collection of 

polynomial or rational functions]. 

or the multi-matrix generalization thereof: 

For any matrices X,Y, . . . of the form there exist matrices A,B,... 

of the form such that Y, . . . , A, B, . . .) = 0. 

The prototype for such matrix decomposition lemmas is the well-known Cholesky 
factorization [HI Theorem 4.2.5]: 

Lemma B.16 (Cholesky factorization) LetX be a real symmetric positive-definite 
n X n matrix. Then there exists a unique lower-triangular real matrix A with strictly 
positive diagonal entries, such that X = AA^ . 

A similar but less well-known result is the following factorization for antisymmetric 
matrices pfT2l[7T[ ITTilfTT7] : 

Lemma B.17 (AJA^ factorization of an antisymmetric matrix) Let X be a 

(real or complex) antisymmetric 2mx2m matrix. Then there exists a (real or complex, 
respectively) 2m x 2m matrix A such that X = AJA^ , where J is defined in iA.15\) . 
In particular, if X is nonsingular, then A G GL{2m). [The form of A can be further 
restricted in various ways, but we shall not need this.] 

In our applications we shall not need the uniqueness of A, but merely its existence. 
Nor shall we need any particular structure of A (e.g. triangularity) beyond lying in 
GL{n) or 0{n) or Sp{2n) as the case may be. Finally, and most importantly, we shall 
not need the existence of A for all matrices X of a given type, but only for those in 
some nonempty open set (for instance, a small neighborhood of the identity matrix). 
We shall therefore give easy existence proofs using the implicit function theorem. It is 
an interesting open question whether our decompositions actually extend to arbitrary 
matrices X in the given classes. 

More precisely, we shall need the following decomposition lemmas in addition to 
Lemma [B.17[ The matrices Imn are defined in (15.641) . 
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Lemma B.18 Let X andY he a (real or complex) mxn matrices (m < n) ofrankm 
that are sufficiently close to the matrix Imn- Then there exist matrices P,R& GL{m) 
and Q G GL{n) such that X = PImnQ O'^d Y = RImnQ~^ ■ 

Lemma B.19 Let X be a (real or complex) m x n matrix (m < n) of rank m that 
is sufficiently close to the matrix Imn- Then there exist matrices P G GL{m) and 
Q & 0{n) such that X = PlmnQ- 

Lemma B.20 Let X he a (real or complex) 2m x 2n matrix (m < n) of rank 2m 
that is sufficiently close to the matrix l2m,2n- Then there exist matrices P G GL{2m) 
and Q G Sp{2n) such that X = Pl2m,2nQ- 

Lemma B.21 Let £ > 1 and ni, . . . , n^+i > 1; and let {Xa}i<a<e be (real or complex) 
Ha X Ua+i matriccs of rank min(na, n^+i) that are sufficiently close to the matrix 
InaUc^+i- Then there exist matrices {Pa}i<a<e+i with Pa G GLijia) such that Xa = 

p T 

Lemmas IB.17l - [B!2n will be needed in Sections I5.4[ 15. 5[ I5.6[ 15.71 and 15. 9[ respectively. 
In addition, Cliolesky factorization could be used in Section 15.31 but we were able to 
avoid it; see the Remark at the end of that section. 

The proofs of these lemmas will all follow the same pattern. First we find an 
explicit pair Xq, Aq (or the multi-matrix generalization) with the needed properties. 
Then we linearize the functions $ in a neighborhood of {Xq,Aq), and we show that 
the tangent space for A at Aq is mapped onto the full tangent space for X at Xq. 
The required existence of A for X in a neighborhood of Xq then follows from the 
implicit function theorem^ For completeness we will also show how Lemmas IB. 161 
and lB.lTt which are known to hold globally, have simple proofs in this "infinitesimal" 
setting. Nearly all these proofs will be easy; only the last (Lemma IB.2ip turns out to 
be slightly tricky. 

Proof of Lemma [B. 161 for X near /. Linearizing X = AA^ in a neighborhood 
of (Xo, Ao) = (/, /) by writing X = I + X' and A = I + A' , we have X' = A' + {A'f 
(plus higher-order corrections that we always drop). Then an explicit solution is given 
by the lower-triangular matrix 



= <; X;, if z > J (B.49) 




We will use the implicit function theorem in the following form: Let U C and V C MP 
be open sets, and let f: U x V be a function {k > 1). Let uq G U and vq G V satisfy 

f(uo,vo) = 0, with {df/du){uQ,vo) nonsingular. Then there exist neighborhoods U' 3 uq and 
V 3 vo such that for all v G V' there exists a unique u G U' satisfying f{u,v) = 0; moreover, the 
map V u is . 
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The rest follows from the implicit function theorem@ □ 



Proof of Lemma [B. 171 for X near J. Linearizing X = AJA^ in a neighborhood 
of (Xo, Aq) = (J, /) by writing X = J + X' and A = I -A' J, we have X' = A'- (A')^. 
Then one explicit solution is given by the strictly lower-triangular matrix 



Another explicit solution is given by the antisymmetric matrix A' = X' /2. □ 

Proof of Lemma IB.181 Linearizing X = PImnQ and Y = RI^nQ'^ in a neigh- 
borhood of (Xo, Yo, ^0, ^0, Qo) = (^mn, ^mn, ^m, ^m, h) bj Writing X = I^^n + X' and 
so forth, we obtain X' = P'Imn + ImnQ' and Y' = R'Imn — ImniQ')^- In terms of 
the block decompositions Imn = {Im,Omx{n-m)), X' = {X[,X^), Y' = {YI,Y^) and 



We can choose Q[i and Q'22 arbitrarily; then the remaining unknowns P', R', Q'12, Q'21 
are uniquely determined. □ 

We shall actually prove the following generalization of Lemma IB.19I 

Lemma B.19' Let X he a (real or complex) m x n matrix (m < n) of rank m 
that is sufficiently close to the matrix Imn, o^nd let Y be a (real or complex) n x n 
symmetric matrix that is sufficiently close to the identity^ matrix In- Then there exist 
matrices P G GL{m) and Q G GL{n) such that X = PImnQ o,nd Y = Q^Q. 

For completeness let us make explicit how the implicit function theorem is used in this case; 
the analogous reasoning for the remaining lemmas can be supplied by the reader. 

If M = {mij)^j^-^ is an n X n matrix, let us write [M]lt = ("^ij)i<i<j<n to denote its lower- 
triangular part. We then use the implicit function theorem (see the preceding footnote) as follows: 
Let u = A = a generic lower-triangular matrix (aij)i<i<j<„, v = [X]lt, uq = vq = [/]lt, f{u,v) = 
[AA^ — XJlt- (Since AA^ — X is manifestly symmetric, it vanishes if and only if its lower-triangular 
part does.) Then {df /du){uo, vq) is the linear map A' i-^ A' + {A')'^; in other words, we have 




(B.50) 




X'2 



P' + Q'n 
Q12 

R' - iQ'nf 



(B.51a) 
(B.51b) 
(B.51c) 
(B.51d) 



duki 



{uo,vo) = (5(y)^(fci) (H- 4i) , 



which is a diagonal matrix with nonzero entries (namely, 1 and 2), hence nonsingular. 
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When Y = In this reduces to Lemma IB.19I 



Proof of Lemma [B. 191 . Linearizing X = PImnQ and Y = Q^Q in a neighborhood 
of (Xo, y^, Po,Qo) = {Imn, In, Im, In) by Writing X = Imn + X' and so forth, we obtain 
X' = P'Imn + ImnQ' and Y' = Q' + {Q')^ . In terms of the block decompositions = 

/ Y' Y' \ 

{Im, Omx{n-m)), X' = {X[, X^), Y' = I , l\rj, with and Y^^ symmetric, and 



x[ 


= ^' + Qn 


(B.52a) 






(B.52b) 




= Qn + iQuf 


(B.52c) 


Y' 

-'22 


— Q22 + ('522) 


(B.52d) 


Y' 

-'12 


= Q'l2 + iQ2lf 


(B.52e) 



We can choose arbitrarily the antisymmetric parts of Q'n and Q22 (e-g- by taking Q[-^ 
and Q22 lower-triangular, or alternatively by taking and Q22 symmetric); then 
the remaining unknowns are uniquely determined. □ 

Similarly, let us prove the following generalization of Lemma IB. 201 

Lemma B.20' Let X be a (real or complex) 2m x 2n matrix (m < n) of rank 2m 
that is sufficiently close to the matrix l2m,2n, CLnd let Y he a (real or complex) 2n x 2n 
antisymmetric matrix that is sufficiently close to J2n- Then there exist matrices P G 
GL{2m) and Q G GL{2n) such that X = Pl2m.2nQ ond Y = JQ. 

When Y = J2„ this reduces to Lemma IB. 201 

Proof of Lemma IB.20I . Linearizing X = Pl2m,2nQ and Y = Q^JQ in a 
neighborhood of (Xq, l^o, ^0, Qo) = {l2m,2n, -hn, hm, hn) by writing^X = /2m,2n + X' , 

Y = J2n + Y', P = hm + P' and Q = hn-'hnQ', we obtain X' = P' l2m,2n-l2ni,2nJ2nQ' 

and Y = Q' — {Q')^- In terms of the block decompositions = {Im,Omx{n-m)), 

f Y' Y' \ 

X' = (X(,X2), Y' = [ / / \T / ) F/]^ and F22 antisymmetric, and Q' = 

\~(^12) ^22/ 

/ Qii Qi2\ have 
\Q'2i Q22/ 



X[ = P'-J2mQ'n (B.53a) 

= -J2n.Q'i2 (B.53b) 

>^i'i = Q'n-iQ'nf (B.53c) 

>^22 = Q22-iQ22f (B.53d) 

^2 = Q'i2-iQ2if (B.53e) 
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We can choose arbitrarily the symmetric parts of Q'n and Q22', then the remaining 
unknowns are uniquely determined. □ 

In preparation for the proof of Lemma IB. 211 it is convenient to introduce some 
simple notation for decomposing rectangular matrices. Given an m x n matrix Y, 
we define the strictly lower-triangular m x m matrix C{Y) by taking the strictly 
lower-triangular part of Y and either deleting the last n — m columns (if m < n) 
or appending m — n columns of zeros (if m > n). Likewise, we define the upper- 
triangular n X n matrix U{Y) by taking the upper-triangular part of Y and either 
deleting the last m — n rows (if m > n) or appending n — m rows of zeros (if m < n). 
It follows immediately from these definitions that 

Y = C{Y)Trnn + ImnU{Y) . (B.54) 



Proof of Lemma IB.211 Linearizing Xa = Pa In^ua+i Pali i'^ ^ neighborhood of 
Xa = and Pa = In^ by writing Xa = /n„n,+i + X'^ and Pa = In^ + Pa, we 

obtain X'^ = -P^/ncn^+i — -^"c^nc+i-Pa+l- 

Let us decompose each P^ as a sum of its strictly lower-triangular part La and its 
upper-triangular part f/„. We therefore need to solve the equations 

X'a = {La + Ua)Inano,+i ^ -^nc.nc+i (-^a+1 + ) ) (B.55) 

where Xi, . . . ,Xi are considered as parameters and Li, . . . , L^+i, Ui, . . . , Ui+i are con- 
sidered as unknowns. But let us prove a bit more, namely that the matrices Ui and 
Li^i can be considered as parameters (i.e. can be chosen arbitrarily). Note first that 
the system (1B.55P can be solved by 

La = ^{X^ — UalriaUa+i + Inana+iLa+l) (B.56a) 
Ua+l = U{X'^-UaTn^n^+^+Tn^n^+^La+l) (B.56b) 

in the sense that any solution of flB.56p provides a solution of fIB.SSp . The equations 
flB.56p appear at first glance to be entangled, i.e. La depends on Ua and vice versa. 
But this is only an appearance, because the operators £ and U "see", respectively, 
only the strictly-lower-triangular and upper-triangular parts of the matrix on which 
they act. Therefore, the system (IB. 561) can be rewritten as 

La = C{X'^ + Tn^n^^,La+l) (B.57a) 
Ua+1 = UiX'^-Uat^n^^J (B.57b) 

But these latter equations can manifestly be solved sequentially for Li, . . . ,Li (given 
Li+i) and for f/2, . . . , f/^+i (given Ui). □ 
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Remark. It is instructive to count parameters and variables in Lemma [6.2 1[ The 
parameters are Xi, . . . ,Xi, and their number is Np = ^^^j^ '^a'^a+i- (The matrices 
Ui and Li^i, which can be chosen arbitrarily, do not count as extra parameters 
because they merely redefine the matrices Xi and X^, respectively.) The variables 
are Li, . . . ,Li and U2, ■ ■ ■ , U^+i, and their number is 

K = ^^i^^ +nl + ...+nU "^-^^^"^^^ + . (B.58) 

Therefore 

N^-Np = > 0. (B.59) 

The — Np extra variables were fixed by our choice of the operators C and U: we 
decided to append m — n columns of zeros to C when m > n, and n — m rows of zeros 
to U when m < n, but we could equally well have inserted an arbitrary strictly-lower- 
triangular matrix of size m — n into C and an arbitrary upper-triangular matrix of 
size n — m into U. This gives {m — n){m — n — 1) / 2 additional variables in both cases, 
which precisely accounts for A^^, — Np. □ 



B.4 Dilation-translation formula 

In Sections 16. II and | 6 .21 we will need the following well-known generalization of the 
translation formula flS.ip : 

Lemma B.21 (Dilation-translation formula) Let P{z) he a polynomial in a sin- 
gle indeterminate z, with coefficients in a commutative ring R containing the ratio- 
nals, and let a and b be indeterminates. Then 

exp + P{z) = P{e^z + ^-^«) (B.60) 

as an identity in the ring R[z,a\\^\\ of formal power series in b whose coefficients are 
polynomials in z and a. [Here the exponential is defined by its Taylor series, as are 
and {e^ - l)/b.J 

In particular, this identity can be evaluated at any nilpotent element b E R, as 
both sides then reduce to finite sums. 

Remark. If b is nilpotent of order 2 (i.e. 6^ = 0), then the formula simplifies 
further to 

exp (^{a + bz)^^ P{z) = p((l + b)z + (1 + |)a) (B.61) 

In our applications in Sections 16.11 and 16.21 we will have 6^ = and also 6a = 0, in 
which case the identity holds even when the ring R does not contain the rationals. 
□ 
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Proof of Lemma IB.211 When a = 0, the formula f lB.60|) states the well-known 
fact that the operator zd/dz generates dilations; it is easily checked by applying both 
sides to z"^. The general case is handled by the change of variables w = z + a/h. □ 



Remark. The formula (]B.60p is a special case of a more general formula for 
operators of the form exY>\^g{z) d/dz\: 

exp(^tgiz)-^^ P{z) = P{z{t-z)) (B.62) 

where 'z{t]z) is the solution of the differential equation dz{t;z)/dt = g(z(t;z)) with 
initial condition ^(O; z) = z. □ 
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